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Recently,  Evans^  h.,s  presented  a  dimensional  analysis  of  the  dependence 

of  strength  on  grain  size,  which  extended  some  earlier  conclusions  established 
2 

by  Singh  et  al.  In  this  analysis,  a  variation  in  local  fracture  toughness, 

Kgff,  with  crack  length,  a,  was  used  to  determine  the  fracture  criticality. 

Specifically,  the  fracture  toughness  was  chosen  to  vary  between  the  singTe 

crystal  and  polycrystal  values.  This  choice  permitted  a  range  of 
c  c 

Kp/Ks  to  be  identified  wherein  stable  crack  growth  preceded  fracture;  causing 
the  strength  within  this  range  to  become  independent  of  the  initial  crack 
length.  This  behavior  resulted  in  the  natural  emergence  of  a  reciprocal 
square  root  grain  size  law  for  the  fracture  strength. 

The  crack  extension  characteristics  were  illustrated  by  selecting  a 
particular  functional  relation  for  the  local  toughness  (fig.  1) 

Keff  *  <s  H  ♦  <X/(UX)]  (1) 

where  k  =  Kp/Kg  -  1,  x  =  a/dw,  d  is  the  grain  size  and  w  is  a  dimensionless 
constant.  This  relation  permitted  a  maximum  in  the  stress  to  be  established 
(fig.  2)  at  the  relative  crack  length  xc  given  by; 


This  stress  maximum  provided  the  final  relation  for  the  failure  strength"; 

af  =  ^  Kp/4/duT  (3) 

However,  this  choice  for  the  toughness  function  only  allowed  a  maximum 
stress  to  exist  for  <>8. 

X 

There  is  a  typographical  error  in  ref.  1  which  has  led  to  the  omission  of 
u)  in  eqns.  (14)  and  (15). 
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The  intent  >f  the  present  note  is  to  demonstrate  that  the  absence 

of  a  stress  maximum  (,  nd  hence,  of  stable  crack  extension)  for  <<3  is 

merely  a  consequence  of  the  particular  analytic  function  (eqn.  1)  chosen 

to  represent  Kg^.  Other  choices  yield  stable  crack  extension  for  much 

smaller  values  of  <,  in  accord  with  the  expectations  of  recent  experi- 

3  4 

mental  observations,  * 

A  particular  limitation  of  the  toughness  function  selected  by  Evans 
occurs  at  small  crack  lengths.  Although  the  function  reduces  to  !C  as 
a-*0,  it  does  not  approach  asymptotically  (only  asymptotic  behavior  at 
large  crack  lengths  was  invoked  in  the  choice  of  the  function).  Yet,  a 
range  of  crack  lengths,  for  which  the  crack  extension  resistance  is  given 
approximately  by  the  single  crystal  value,  K^,  is  expected  for  real 

systems.  It  can  be  readily  demonstrated  that  the  choice  of  a  reasonable 

toughness  function,  that  is  also  asymptotic  at  small  crack  lengths,  provides 

c  c 

the  requisite  stable  crack  extension  for  small  values  of  KJ/k;.  A  perti¬ 
nent  function  has  the  form; 

Kgff  =  [1  +  K(l-e'x2)]  (4) 

This  function  is  compared  with  eqn.  (1)  in  fig.  1.  Differentiation  yields 

the  critical  normalised  crack  length; 

X2  2 

[(<+!)/<]  e  c  -  1  =  4XC  (5) 

Inspection  of  eqn.  (5)  indicates  that  Xc  exists  for  all  positive 
values  of  k,  .A  regime  of  stable  crack  extension  will  thus  occur  for  any 
value  of  Kp/Kg,  as  illustrated  in  fig.  2.  It  is  noted  that  a  stable  growth 
region  is  encouraged  by  an  extended  range  of  uniform  crack  extension  resis¬ 
tance,  Kgf f  *  K^.  The  crack  extension  stress  diminishes  continuously 
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w  thin  the  single  crystal  range  and  hence,  permits  the  stress  to  increase 
when  the  crack  encounters  a  regime  of  increasing  crack  growth  resistance. 

However,  it  is  emphasized  that  a  region  of  increasing  crack  extension 
stress  does  not  ensure  that  stable  crack  extension  will  be  observed.  It 
is  also  required  that  the  initial  crack  length  be  large  enough  to  allow 
the  initial  extension  stress  to  be  less  than  the  stress  maximum,  as 
discussed  in  ref.  1.  An  extensive  single  crystal  region  of  crack  growth 
resistance  is  thus  conducive  to  a  crack  length  independent  fracture  stress 
only  if  the  pre-existent  crack  is  substantially  extended  into  the  single 
crystal  region. 

The  potential  for  stable  crack  growth  has  thus  been  demonstrated  for 
any  system.  Its  existence  in  a  particular  system  will  depend  upon  the 

functional  form  of  the  local  toughness  and  the  initial  crack  length.  A 
complete  description  of  strength  behaviors  thus  requires  an  experimental 
determination  of  the  local  toughness.  One  of  us  (A.  V.  Yirkar)  is  currently 
engaged  in  the  study  of  toughness  on  the  local  level. 


ACKNOWLEDGEMENTS:  The  above  work  was  supported  by  the  National  Science 
Foundation  under  Grant  No.  DMR-78-10016.and  the  Office  of  Naval  Research 


under  contract  No.  N0014-79-C-0159. 


-4- 


REFERENCES 


1.  A.  G.  Evans,  "A  Dimensional  Analysis  of  the  Grain  Size  Dependence  of 
Strength,"  J.  Am.  Ceram.  Soc.,  63^  [1-2],  115-116  (1980). 

2.  J.  P.  Singh,  A.  V.  Virkar,  D.  K.  Shetty  and  R.  S.  Gordon,  "Strength- 
Grain  Size  Relations  in  Polycrystalline  Ceramics,"  J.  Am.  Ceram.  Soc. 
62,  [3-4],  179-183  (1979). 

3.  Anil  V.  Virkar,  "Electrical -Mechanical  Analog  Applied  to  Crack  Growth 
Studies  in  Glass-Glass  Ceramic  Composites,"  J.  Am.  Ceram.  Soc.,  63_ 
[3-4],  219-223  (1980). 

4.  Anil  V.  Virkar,  "Simulated  Strength  -  Grain  Size  Study  Using  Glass- 
Glass  Ceramic  Composite  System,"  submitted  to  Journal  of  Materials 
Science,  (1979). 


CHAPTER  II 


MARTENSITIC  TRANSFORMATIONS  IN  ZIRCONIA: 
PARTICLE  SIZE  EFFECTS  AND  TOUGHENING 


by 

A.  G.  Evans,  N.  Burlingame,  M.  Drory  and  W.  M, 


Kriven 


Department  of  Materials  Science  and  Mineral  Engineering, 
University  of  California, 

Berkeley,  CA  94720 


MARTENSITIC  TRANSFORMATIONS  IN  ZIRCONIA: 


PARTICLE  SIZE  EFFECTS  AND  TOUGHENING 
by 

A.  G.  Evans,  N.  Burlingame,  M.  Drory  and  W.  M.  Kriven 

Materials  Science  and  Mineral  Engineering 
University  of  California 
Berkeley,  CA  94720 


Abstract 

The  incidence  of  martensitic  transformations  in  systems  containing 
ZrO^  particles  (or  precipitates)  depends  upon  the  size  of  the  particles 
The  origin  of  this  size  effect  is  demonstrated  to  reside  in  the  twinned 
or  variant  structure  of  the  transformed  phase,  by  virtue  of  a  surface 
area  related  strain  energy  term.  Predicted  critical  particle  sizes  cor 
relate  quite  well  with  observations  on  two  Zr02  systems.  The  particle 
size  effect  can  also  be  incoroorated  into  toughening  analyses,  in 
order  to  predict  both  octimum  toughening  conditions  and  toughening 
trends.  The  predictions  again  correlate  favorably  with  toughness  date, 
for  ZrO^  containing  materials. 
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I,  INTRODUCTION 

Martensitic  transformations  (involving  a  tetragonal  to  monoclinic 

crystal  structure  change)  have  been  observed  in  zirconia  precipitates'1  ’ 

4  5 

or  particles  ’  .  The  incidence  of  this  transformation  depends  upon  the 

size  of  the  particle  :  a  phenomena  that  has  not  yet  been  adequately 

explained.  The  initial  intent  of  the  present  paner  is  to  identify  the 

origin  of  the  size  effect  and  to  provide  a  basis  for  the  quantitative 

prediction  of  the  critical  transformation  condition.  These  results  cc,-:ain 

implications  for  the  increase  in  fracture  tounimess  that  can  be  achievcu  in 

the  presence  of  a  size  distribution  of  predominantly  tetragonal  ZrOo  par- 
5  7 

tides  ’  .  The  interpretation  of  toughening  effects  in  systems  containing 

zirconia  particles  thus  constitutes  a  second  theme  of  the  paper. 

Martensitic  transformations  in  zirconia  can  occur  in  accord  with 

either  of  the  two  lattice  invariant  deformation  processes,  twinning  or 
8  9 

slip  ’  .  A  relatively  large  shear  is  expected  to  accompany  transformation 
for  most  perm.issable  lattice  invariant  processes  (viz.  a  shear  strain  along 
the  habit  plane  of  d4';).  This  large  shear  strain  results  in  the  formation 
of  a  series  of  variants  (Fig.  1);  especially  when  the  transformation  occurs 
within  particles  embedded  in  a  matrix  (such  that  the  macroscopic  shape  oe- 
formation  is  subject  to  constraint).  Transformed  particles  thus  contain  a 
series  of  sheared  plates  in  which  alternate  plates  have  experienced  sheer 
deformations  of  opposite  sign  (Fig.  1).  A  significant  macroscopic  shear 

\ 

is  not  expected,  therefore,  whenever  there  are  either  a  large  number  of 
variants  or,  in  general,  for  an  even  number  of  variants. 

~A  size  effect  has  been  anticipated  for  situations  in  which  the  transfor~a- 
tion  yields  microcracks°  (notably,  single  phase  oolvcrystal 1 ine  Zr02).  ^o’:- 

ever,  the  transformation  of  isolated  Zr02  particles* i s’  not  usually  acccmar-.u 
by  microcracking1 .3.  An  alternate  explanation  must  ce  sought  for* tnese 
systems . 
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The  strain  energy  change  that  accompanies  the  martensitic  transfor¬ 
mation  should  be  strongly  influenced  by  the  variant  or  twin  structure  of 
the  martensite,  Evidently,  for  a  particle  of  specified  size,  the  higher 
the  variant  (twin)  density,  the  lower  is  the  strain  energy,  and  t.ne  larger 
the  interface  energy  (associated  with  the  parent/product,  or  twin,  inter¬ 
face),  The  energy  changes  attributed  to  the  twinned  structure  of  the  mar¬ 
tensite  are  considered  to  be  the  source  of  the  transformation  size  effect. 

An  analysis  of  the  strain  energy  of  spherical  twinned  particles  r,as 
been  presented  by  Kato  et  al^,  illustrating  the  diminution  in  strain 
energy  that  accompanies  an  increase  in  twin  density.  This  solution  is 
an  essential  source  of  the  quantitative  trends  in  strain  energy  pertinent 
to  the  analysis  of  size  effects.  This  solution  is  augmented  in  the  pre¬ 
sent  paper  by  a  two  dimensional  analysis  of  the  stress  and  strain  energy 
distributions  in  twinned  (or  variant)  martensites,  which  provides  tne 
additional  information  (regarding  strain  energy  localization)  needed  tc 
address  issues  related  to  size  effects.  Both  analyses  pertain  only  to 
systems  with  homogeneous  elastic  constants;  the  critical  size  predictions 
are  thus  restricted  to  such  systems  (a  close  approximation  is  a  system 
of  Zr02  precipitates  in  a  cubic  ZrO^  matrix).  However,  approximations 
pertinent  to  the  analysis  of  transformations  in  the  presence  of  elastic 

inhomogeneity  are  suggested.  The  critical  transformation  sizes  predicted 

1  3 

by  the  analysis  are  compared  with  results  for  ZrO^  ’  and  for  A1  ^O^/^rOg 
alloys'1^ , 

The  size  effect  also  has  a  significant  influence  on  the  toughening 
that  can  be  experienced  in  the  presence  of  an  array  of  tetragonal  ZrC^ 
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1  ?  10 

particles-  A  previous  toughening  analysis  “*  °  is  extended  to  include 
a  size  distribution  of  particles:  to  yield  predictions  of  the  optimum 
toughening  and  of  toughening  trends.  The  predictions  are  compared  with 
experimental  results^’7. 


2.  STRESS  AND  STRAIN  ENERGIES  ASSOCIATED  l.'ITH  TWINNED  PARTICLES 

The  stress  distributions  that  develop  in  martensite  plates  contain¬ 
ing  twins  or  variants  can  be  determined  straightforwardly  (in  the  absence 
of  elastic  modulus  inhomogeneity)  by  adopting  an  Eshelby  approach"*4  (Fig. 
2).  The  stresses  within  the  twinned  plates  are  determined  from  the  sum 
of  the  uniform  shear  stresses  needed  to  restore  the  unconstrained  twins 


to  their  original  shape  (Fig.  2)  and  the  non-uniform  stresses  that  arise 
from  applying  the  body  forces,  at  the  twin  boundaries  and  parent/product 
interfaces,  needed  to  establish  continuity  of  stress6, (Fig.  2).  The 
stresses  within  the  parent  phased  derive  exclusively  from  the  body  forces 
The  spatially  dominant  stresses  induced  by  twinning  are  the  shears  x  . 
The  shear  stresses  at  a  location  (y,z)  with  resneci;  to  the  twin  interface 
(Fig.  2c),  that  derive  from  a  line  body  force  of  magnitude  P,  at  a  dis¬ 
tance  x  from  the  interface  (Fig.  2),  are  given  by"*7; 


[(x+z)2+y2] 


(1-V)  +  2  (  1  +V  )  — -  y 
[(x+zr+y2] 


Normal  stresses  a  and  c  also  develop.  These  are  singular  at  the  twin/ 
xx  yy 

habit  intersection,  but  decay  very  rapidly  with  distance  (a  logarithmic  singu¬ 
larity)  *6,  The  normal  stresses  are  essential  to  considerations  of  micro- 
cracking"*6,  but  constitute  a  secondary  contribution  to  the  strain  energy. 

The  normal  stresses  are  thus  neglected  for  present  purposes. 
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where  v  is  the  Poisson  ratio.  Ti  .■  body  forces  needed  to  establish  stress 
continuity  are  given  by; 

P  =  -(uYT/2)dx 

where  u  is  the  shear  modulus  and  is  the  unconstrained  shear  strain 
associated  with  each  twin  (Fig.  2a).  The  stresses  around  a  centrally 
located  twin  can  be  deduced  by  integrating  the  body  forces  over  all  twin 
and  parent/product  interfaces.  The  stresses  within  the  parent  phase  for 
a  particle  with  length  L,  width  •'  and  twin  spacing  a  are  given,  for  >„  =0.2,  by; 


n  =  L/4d  1 

r  2~+y/d 

2-1  -y/d 

'xy  1 

PV  '  4~ 

E  1 

1 

•=2)ds  *  / 

1 

n  =  1 

2  —  1 +y/d 

2^-2 -y/d 

2n-l+y/d 

2-y/d 

z/d+l/d 

r 

-  / 

<V?>da 

-  /  (=r 

c2)d2  +  2  J 

2n-2+y/d 

2n-l -y/d 

z/d 

( 3 ) 

\  *  i 


where , 


A 

n 


(z/d)[0.8(z/d)2+3.252] 


’2 


(z/d-i/d)[0.8(z/d--/d)2+3.252] 

[z/d-t/d)2+S2]2” 


the  ;  terms  are  given  by; 


y[0.8y2+3.232] 

tiW? 


such  that  ■  is  ( 2 —  1  -w )  for  ,  (2--2-w)  for  :2 ,  (2--l^w)  *‘or  and  ,  2~-w) 

for  and  i  =  x/d. 
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The  strf  ,ses  within  the  twinned  particles  are  deduced  from  Eqn.  (3) 

by  inserting  the  appropriate  negative  values  for  z/d,  and  superimposing 

a  uniform  shear  stress  (of  unit  normalized  magnitude). 

The  resultant  stresses  are  plotted  in  Fig.  3  at  several  locations,  y/d. 

The  same  results  pertain  for  all  values  of  particle  length,  L  >10d,  and  for 

all  particle  widths,  l  >3d.  Two  important  characteristics  are  noted.  Firstly, 

the  zone  of  significant  stress  is  confined  to  a  region  close  to  the  plane  of 

termination  of  the  twins,  i.e.,  the  parent/oroduct  interface.  This  behavior 

arises  because  the  alternating  shears  associated  with  adjacent  twins  tend 

to  eliminate  long  range  stresses  (a  result  that  should  be  expected  from 

St.  Venant's  princiole).  Secondly,  the  stresses  are  a  unique  function  of 

the  normalized  distance  from  the  interface,  z/d;  independent  of  the  absolute 

twin  spacing,  the  twin  dimension,  and  the  particle  dimension,  (at  least  for 

t/d>3).  This  localization  of  the  stresses  and  the  scaling  with  d  had  pre- 

1 8 

viously  been  anticipated  by  Roitburd  and  Khachaturyan  .  The  stresses  de¬ 
duced  from  Fig.  3  can  be  approximately  exoressed  (for  0.5  y/d -0.1)  by  the 
analytic  function; 


=  F(z/d)  *  1.5(z/d)°'7e"2'z/d) 


,1. 


(4) 


For,  0.1>y/d<0,  the  stresses  decrease,  monotonically,  approaching  zero  at 
y/d=0.  A  different  result  will,  of  course,  pertain  for  twins  adjacent  to 
the  extremities  of  the  particle. 

Since  the  stresses  are  confined  to  a  zone  close  to  the  particle  inter¬ 
face,  it  might  be  anticipated  that  the  component  of  the  strain  energy 
associated  with  the  twinning  be  more  closely  related  to  the  surface  area 
of  the  particle  than  its  volume  (note  that  the  strain  energy  derived  from 
the  macroscopic  shape  change  is  solely  dependent  on  the  particle  volume  V), 


/ 
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especiall:  for  large  numbers  of  twins,  A  rough  estimate  of  this  strain 
energy  term  for  a  spherical  particle  of  radius  a  is  obtained  by  considering 
that  the  shear  stresses  in  each  twin  are  given  by  Eqn.  (4),  This  estimate 
of  the  strain  energy  change  is; 

F2(z/d)d  (z/d) 

-?/2d  ,,s 


Inserting  the  approximate  F(z/d)  suggested  by  Eqn.  (4),  the  strain  energy 
becomes, 


m  _ 


'('T> 


(4-a2)d 


/ 


where  n  is  the  number  of  twins  in  the  particle  (~=2a/d).  This  is  necessarily 
a  lower  bound  result  because  the  contribution  from  the  normal  stresses  has 
been  neglected.  Recognition  of  the  special  significance  of  the  surface  area 
provides  a  useful  perspective  for  interpreting  and  extending  the  numerical 
solution  of  Kato  et  al^.  For  a  large  number  of  tv;ins  it  might  be  anticipated 
that  the  normalized  energy,  )  ,  becomes  inversely  proportional  to 

the  number  of  twins,  as  exemplified  by  Eqn.  (6).  However,  in  order  to  permit 
the  existence  of  more  complex  behavior  at  limited  twin  densities,  a  more 
general  relation  for  the  strain  energy  can  be  obtained  by  adopting  the  func¬ 
tion; 
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A4>-i 


1 


i(yr)v 


x  ,  n 

2+ 


(7) 


where  X ^  and  X2  are  constants.  A  comparison  of  Eqn.  (7)  with  the  calcu¬ 
lations  of  Kato  et  al"*0  (Fig.  4)  indicates  a  good  correlation,  with 
X-j  ^0.64  and  X2^2.4.  The  result  for  large  ", 


Ap 


0.64 


u(yt)‘ 


(8) 


is  in  modest  agreement  with  the  estimate  obtained  above,  by  focussing 
on  the  stresses  in  the  vicinity  of  the  particle  interface. 

3.  TRANSFORMATION  THERMODYNAMICS 

The  changes  in  thermodynamic  potential  that  accompany  the  formation 
of  a  twinned  martensite  are  the  mechanical  potential,  Ac™  (strain  energy 
and  interaction  energy),  the  surface  work,  Aos,  and  the  chemical  potential, 
AFc_  The  latter  is  independent  of  the  extent  of  twinning,  and  for  spherical 
particles  of  radius,  a,  the  chemical  free  enemy  change  per  particle  is 
simply; 

AFC  =  -(4/3ba3AFo  (9) 


where  AFQ  is  the  Helmholtz  free  energy  difference  between  unit  volume  of 
the  two  phases.  The  surface  work  reflects  chances  that  occur  at  the  mar¬ 
tensite  interface  as  well  as  the  formation  of  twin  (variant)  boundaries. 
If  the  twin  spacing  is  d,  the  total  surface  energy  change  Ab$  for  a 
spherical  particle  is; 


Abe 


=  4,a2r>2’a2  (n'' >  rt 


00) 
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where  7.  is  the  interface  energy  and  I\  is  the  twin  boundary  energy. 
1 

.Evidently,  for  a  large  number  of  twins,  d<<a,  Eqn.  (10)  reduces  to; 


The  mechanical  energy  change  can  be  considered  tc  consist  cf  two 

m 

components:  that  due  to  the  macroscopic  shaos  change  :.om  and  that 

associated  with  the  twinning  Aq-p  The  macrosconic  component  under  cordi- 

tions  of  constant  applied  stress  P.^  includes  a  strain  energy  and  an 

1  3 

interaction  energy(with  the  applied  stress)  and  is  given  by  ; 

<  *  ra3|“v  [o.14£iV-pA]  +  [o.21Eedef -pTet]l  (12 

where  AV  is  the  volume  stra in ,  e|j  is  the  deviator ic  component  of  the 
average  particle  transformation  strain  and  o  ?.f,i  pp  are,  respectively 
the  dilational  and  deviatoric  components  of  the  anolied  stress..  The 
twin  component  is  obtained  directly  from  Eon .  (7)  fS; 

AOf  =  (V3)-a3EYT2C0. 13/(1. 2+a/d)]  (13) 


This  component  does  not  involve  a  significant  interaction  energy  with  the 
applied  stress,  because  of  the  alternating  character  of  the  twinning  shear 
strains.  It  is  not  generally  permissable,  of  course,  to  sum  strain  ere-gy 
terms,  unless  they  derive  from  different  constituents  of  the  applied  stra:-. 
Hence,  the  twin  component  (which  emanates  frr:  t.se  shear  stresses)  car  02 
combined  directly  with  the  dilational  constituent  of  the  macroscopic  rote”t 
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t>u t  should  not  be  combined  with  the  deviatoric  const'  tuent.  However,  as 
noted  above,  typical  variant  structures  of  transformed  ZrO^  particles  in¬ 
dicate  a  small  macroscopic  shear  deformation.  It  would  thus  appear  per- 
missable  to  neglect  the  deviatoric  part  of  the  macroscopic  strain  energy, 
for  present  purposes,  thereby  averting  this  difficulty.  Additionally,  it 
is  noted  that  the  spatial  locations  of  the  stresses  created  by  twinninq 
(localized  near  the  interface)  are  quite  differenct  from  those  associated 
with  the  macroscopic  shape  change  (uniformly  distributed  throughout  the 
particle),  Any  error  incurred  if  a  simple  summation  of  the  mechanical 
energies  were  invoked  should  thus  be  quite  small. 

The  total  Dotential  change  associated  with  transformation  can  now  be 
expressed  as; 

— ^--y-  =  -aF  +rt/d  +  3rH/a  +  0.14E  :v2-paav 

(4/3)ra3  . 

+  0. 1  3E'<t2/(1  ,2+a/d)  (14) 

This  basic  result  can  be  used  to  examine  the  existence  of  a  driving  force 
for  the  transformation  and  thereby,  to  estimate  size  effects,  toughening, 
etc.  The  four  variables  that  influence  transformation  in  a  given  material 
are  the  chemical  free  energy  AFQ  (which  is  dependent  on  temperature  and 
solute  content),  the  particle  size,  the  twin  spacing  and  the  applied  stress. 
The  only  variable  that  is  not  readily  specified  a  priori  is  the  twin  spacing. 
The  twin  structure  presumably  evolves  in  the  nucleative  stage  of  the  trans¬ 
formation  and  hence,  a  minimization  of  A:  vit:o  respect  to  d  (at  constant  a 

A 

and  P  )  does  not  necessarily  provide  a  useful  prediction  of  the  expected 
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twin  spacing.  Such  an  analysis  (3A:/3d=0)  would  suggest  that 


1.9Cta) 


1/2 


YTE1/2[1-4.6:^/Y7(Ea),/t] 


1/2 


,1/2- 


(15a) 


or,  for  a  >>  d. 


d 


(15b) 


Experimental  results  obtained  by  Kato  et  al^  for  iron  particles  in  a 

Cu-Fe  alloy  are  not  consistent  with  this  prediction-  rather,  d  would 

appear  to  be  approximately  independent  of  the  particle  size.  This 

experimental  result  is  preferred  for  all  subseouent  developments  in  this 

paper  because  preliminary  observations  suggest  a  similar  benaviour  for 
1119 

ZrC^  ’  .  Refinements  that  include  a  size  dependence  of  d  could  be  incorporates 

later,  if  it  is  substantiated  that  a  significant  size  effect  exists  in  toe 
systems  of  interest.  However,  it  is  noted  here  that  the  twin  spacing 
should  be  independent  of  the  chemical  free  energy  change. 


A  two-dimensional  nucleation  analysis  involving  a  simple  twin  or  variant 
has  been  reported  by  Roitburd  and  Kurdjumovc®.  At  the  critical  nucleation 
condition  a  twin  spacing/length  relation  identical  to  Eqn.  (15),  excect 
for  the  numerical  coefficient,  was  derived.  Then,  bv  assuming  that  the 
twin  (variant)  extended  fully  across  the  particle  before  broadening  to 
the  final  dimension,  a  twin  spacing  proportional  to  the  twin  length  was 
predicted.  This  prediction  is  even  less  consistent  with  the  experimental 
results.  Further  studies  of  this  phenomenon  are  clearly  demanded. 
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4,  PARTICLE  SIZE  EFFECTS 

Particle  size  effects  on  transformation  emerge  by  considering  the 

condition  wherein  the  total  free  energy  chanae  becomes  incrementally  negative. 

This  condition  will  yield  a  minimum  possible  value  for  the  critical  particle 

size,  ac .  ,  at  which  transformation  will  initiate.  This  estimate  of  the 
min 

critical  size  is  most  pertinent  when  the  nucleation  barrier  is  relatively 
small,  i.e.tne  free  energy  at  successive  staaes  of  partial  transformation  is 
only  marginally  in  excess  of  the  free  energy  of  the  fully  transformed  p?r- 
t i c  1  e o  Justification  for  the  neglect  of  an  '^reciable  nucleation  barriey' 
can  only  be  effectively  estaolished  (for  a  specific  system)  by  a  comprehen¬ 
sive  comparison  of  the  predicted  and  observed  critical  transformation  dimen¬ 
sion  and  its  dependence  on  composition,  temperature  and  elastic  stiffness. 

(It  is  noted  that  this  comparison  is  most  readily  achieved  with  ceramic 
systems,  because  they  exhibit  negligible  relaxation  of  the  elastic  strains 
by  dislocation  activity  near  the  particle).  Some  favorable  preliminary 
comparisons  for  ZrC^  will  be  presented  in  this  section. 

Setting  the  total  thermodynamic  potential  It  in  Eqn.  (14)  to  zero,  the 
following  relation  for  the  lower  bound  critical  size  obtains; 


0 . 1 3E»  Td 

l,2d+ac,n 

min 


;FQ+PAiV  -  P.14EAV2 


If  the  twin  spacing  is  considered  to  be  insensitive  to  particle  size  ,  as 


Critical  transformation  dimensions  considerably  in  excess  of  the  values 
based  on  the  free  energy  of  the  fully  transformed  oarticle  would  suggest 
a  size  effect  based  on  statistical  arguments,  through  a  spatial  distribu¬ 
tion  of  nucleation  sites. 

*f*T 

It  is  noted  that  the  present  analysis  would  yield  a  size  effect  (attributabl 
to  strain  energy  localisation  due  to  variant  formation)  provided  that  6  *  an 
with  n<l .  It  is  not  imperative  that  d  be  independent  of  a, 
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the  limited  experimental  observations  suggest,  then  the  critical  particle 
size  becomes; 


3I\  +  0.l3EyT2d 
+  d(AFQ  +  pAAV  - 


0.1 4c  AV  ~ ) 


-  1.2 


(17) 


2-2  -  ] 

When  the  interface  energies  are  relatively  snail  (r./Evy d< 1 0  ,rt/dt,F  <1:3  ), 

the  relation  for  the  critical  size  reduces  to; 


amin 

d 


_c 

2 


0.13Eyt2 

J.F0+PA^V -0 . 1 4E_  V2 


-  1.2 


(18) 


which,  in  the  absence  of  an  applied  stress,  becomes: 

0.27(yT/AV)2 

p  =  - 1 -  -  2-4 

c  ? 

IF  /ElV  -0.K 
o 


(19) 


A  specific  number  of  twins  (or  variants)  is  thus  nrsdicted  at  the  critical 

4- 

condition. ' 

The  preceding  analysis  pertains  to  systems -v/ith  homogeneous  elastic  and 
thermal  expansion  properties.  Modifications  that  incorporate  elastic  inhomo¬ 
geneity  and  thermal  expansion  mismatch  effects  are  of  considerable  importance. 
The  effect  of  modulus  inhomogeneity  is  difficult  to  anticipate  because  only 
the  component  the  strain  energy  associated  with  the  macroscopic  shape 
change  can  be  readily  analyzed.  However,  since  the  stresses  associated  v.'itii 


tSince  it  is  not  yet  known  whether  the  variant  snacing  can  be  changed  for  a 
given  material,  it  should  not  be  assumed  that  a  critical  number  of  variants 
could  be  achieved  for  a  variety  of  actual  particle  radii. 


* 


the  matrix  and  particle  <E>  may  be  used  to  afford  an  approximate  measure 
of  the  strain  energy  change,  viz; 


Afry  ~  (4/3)  ira3<E>Yy  [0.1 3/(1. 2  +'a/d)] 


The  mechanical  energy  change  associated  with  the  particle  dilation,  for 
v  =  0.2,  is13 


£Um  =  ( 4/ 3  )-rra3  Ep6AV 


3. 6(1+3) 


where  S  is  the  ratio  of  the  elastic  modulus  of  the  matrix  E  to  that  of 

m 

the  particle,  E  .  Modifying  the  preceding  analysis  in  accord  with 
these  energy  terms  yields  a  new  value  for  the  variant  density  at  the 
critical  size; 
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The  above  predictions  can  be  compared  with  experimental  results  ob- 
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tained  for  partially  stabilized  ZrC^  ’  ,  and  for  AiyyZrC^  .  Transforma¬ 
tion  in  partially  stabilized  ZrC^  results  in  the  formation  of  variants 

1  3 

with  about  four  variants  within  each  particle.  *  The  pertinent  transfor- 

12 

mation  parameters  are  :  AFQ^230MPa  (for  transformation  close  to  room 
temperature  in  partially  stabilized  ZrQ^,  which  is  the  critical  condition 
of  present  interest),  EM70GPa,  4VM1.Q58.  The  shear  in  the  habit  plane 
Yj  is  ^-0 .14  '  and  (for  an  even  number  of  variants)  there  is  no  macro¬ 
scopic  shear  deformation  of  the  particle.  The  measured  variant  spacing^ 

is  Ml  x  10~^m.  Hence,  for  typical  upper  bound  values  of  the  interface 

-2  7 

energies,  r.'vr.'vlJm  ,  the  magnitudes  of  the  quantities,  r./EyT  d  and 

t  1  1  1 

-  -3  -2 

ydAFo  are  Aj2  x  10  and  ^2  x  10  respectively.  The  interface  energy 

2 

terms  in  Eqn.  (17)  are  thus  of  negligible  magnitude.  Noting  that  AFQ/EAV 
is  Ml. 4,  substitution  of  the  remaining  transformation  parameters  into 
Eqn.  (17)  indicates  that  the  critical  number  of  variants  for  transformation 
in  the  osence  of  an  external  stress  is  ^4.  The  prediction  is  thus 

I  “7 

consistent  with  the  observed  variant  structure  of  transformed  particles  . 

It  is  also  noted  that  an  applied  stress  will  only  significantly  reduce  the 

critical  variant  density  (i.e.,  to  M?)  if  the  hydrostatic  component  is  in 

excess  of  ^3GPa.  Such  large  stresses  can  only  be  encountered  in  the 

21 

vicinity  of  a  crack  tip  . 

For  the  AlgOyZrOg  system,  the  modulus  ratio  S  is  ^2.5  and  hence,  the 
number  of  variants  at  the  critical  size,  derived  from  Eqn.  (23),  is 


°-25Vt2 

ifo-0.2Ep4V2 


2.4 


(24) 


J6- 


The  (  fective  volume  strain  in  the  ZK^/A^O-j  svsten  is  smaller  than  in 
a  Zr02  matirial  because  of  the  thermal  expansion  mismatch.  La,  such  tnat 

£V  ^  ~  L\l-3LaAT 

where  iT  is  the  cooling  range.  The  incorporation  of  the  effective  volume 
change. is  equivalent  to  recognizing  that  the  system  is  subject  to  a 
strain  energy  prior  to  transformation  that  derives  rrom  the  expansion 
mismatch  and  hence,  that  the  change  in  strain  enemy  uDon  transformation 
is  reduced  (when  La  is  positive).  Insertion  the  effective  volume  change 
and  the  other  transformation  oarameters  into  Eqn.  {?.-),  the  number  of 
variants  at  the  critical  size  is  ".6.  Observations  (Fig.  5)  indicate 
6-10  variants  in  critical -sized  particles. 

Particle  size  effects  predicted  by  the  thermodynamic  requirements 
for  complete  transformation  are  thus  consistent  with  observations  of 
martensitic  transformations  in  two  ZK^  systems.  A  basis  for  further  study,  em 
phasizing  effects  of  temperature  and  solute  content  on  the  transformation,  is  t 
established.  Also,  a  more  comprehensive  analysis  of  toughening  than  has 
previously  been  possible  can  now  be  attempted,  as  discussed  in  the  follow¬ 
ing  section. 

5.  TRANSFORMATION  TOUGHENING 
5.1  Analysis 

The  hydrostatic  component  of  the  stress  needed  to  induce  transforma¬ 
tion  (neglecting  interface  energy  terms)  is  obtained  from  Eqn.  (14)  as; 


-17- 


A  ?  0.27EyT 

p  lV  =  0.14E-/-IF  +  — - - 

o 


(26) 


2.4+n 


22 

Noting  that  the  hydrostatic  stress  near  a  crack  tip  is  ; 


pA  .  MlivI  cos  ( 6/2 ) 
3/2  nr 


(27  ) 


where  K  is  the  stress  intensity  factor  and  (r,r)  are  the  coordinates  with 
respect  to  the  crack  tip,  the  distance  rc  from  the  crack  tip  at  which 
a  single  particle  will  transform  can  be  deduced,  for  ■.  =  0.2,  as; 


2~rr  =  8.8  < 

c 


,.v  (2.4+n)(2.4+-  }  I 

c-  cos(e/2)  * 


Ey, 


(nc-n) 


(2 


n  ' 


This  prediction  makes  the  implicit  assumption  that  each  particle  ex¬ 
periences  a  sufficient  shear  stress  that  the  transformation  may  nucleate 
when  the  thermodynamic  conditions  are  satisfied.  This  assumption  re- 
auires  further  investigation. 

The  tougheni'in  7-  associated  with  the  transforming  particles  may 
now  De  determined,  using^  ; 


?T  •  27P*rTVf/(4/3)ra3+:o 


(29) 


where  is  the  volume  fraction  of  particles,  Tq  is  the  intrinsic  tough¬ 
ness  of  the  matrix,  id*  is  the  free  energy  change  in  the  absence  of  an 
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aDplied  stress  and  r-r  is  the  transformation  zone  width,  given  by  , 


( 
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.(Vf)rc 


(3D) 


where  r  is  the  zone  dimension  ;hat  coincides  with  the  maximum  zone  width 

(8%tt/3)  .  The  quantity  u  is  a  proportional ity  constant  that  allows 

for  the  relaxation  of  the  stress  field  by  the  prior  transformation  of 

particles  closer  to  the  crack  tip  and  for  a  partial  reverse  transformation 

of  particles  near  the  periphery  of  the  transformation  zone.  A  much  more 

extensive  analysis  than  that  conducted  herein  would  be  needed  to  evaluate 

u-  However,  for  present  purposes,  ^  is  estimated  bv  comparing  predicted 

21 

with  observed  transformation  zone  dimensions  .  The  transformation  zone 

size  for  a  partially  stabilized  ZrO^  with  two  complete  transformation 

variants  within  each  particle  and  a  fracture  toughness3of  4.9MPa»m  is 

predicted  from  Eqn.(23)  to  be  1.5um.  This  compares  with  a -measured 

21 

zone  size  of  n,0.6  urn  :  suggesting  that  the  relaxation  parameter  ^  is 
n-0.4.  This  magnitude  is  assumed  herein  to  pertain  for  all  conditions. 
However,  it  is  recognized  that  there  could  be  a  sicinificant  influence 
of  the  particle  volume  fraction  on  u:  an  effect  that  is  not  included 
in  the  present  analysis.  Substituting  r^  from  Eons.  (23)  and  (30) 

(the  maximum  zone  width)  into  Eqn.  (29),  the  tounhening  becomes 


‘T 


6.6 

r r 


K2(AV/yt)2 


(2.4+n)(2.4+nc) 


:Vn) 


(31) 


Inserting  the  requirement  that  K  at  the  criticality  is,  *E ?j,  we  obtain; 


r 

i 

0 


(32  ) 
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where 


-2Vf(,V/,T)2  (2-4")12-4tr'c>  ... 


<Vn) 


As  noted  in  the  previous  section,  lV/y-  is  0.41.  The  toughening  ex¬ 
pected  for  ZKU  is  thus  expressed  (wi th  ujy.0.4)  by, 


0.1?.  V, 


(2.4+n)  (2.4+r,c) 


For  the  Al203/Zr0?  system,  the  modulus  mismatch  (3=2. 13)  causes  ;  to  increase, 
by  a  factor  of  M  .2. 

The  toughening  anticipated  by  Eqn.  (32)  pertains  when  the  particle 
sizes  are  uniform.  Significant  deviations  can  be  expected  when  a  size  dis¬ 
tribution  of  particles  exists.  If  the  particle  size  distribution  is  charac¬ 
terized  by  the  frequency  function,  q>(a)da,  then  the  toughening  given  by 
Eqn.  (30)  becomes; 


2K2(aV/yT)2 


(2.4+nc) 

<a3> 


•/ 


(2.4d+a! 

(a'-a) 


c(a)da 


(33) 


Assuming  an  extreme  value  distribution  for  $(a)da  (since  the  larger 
particles  are  of  principal  interest23), 

at 

<$(a)da  =  k  exp 


where  aQ  is  the  scale  parameter  and  k  is  the  shaoe  parameter,  the  touchen- 
ing  becomes; 


I 
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:  T 


=  2r 


(A V/7T)  Vf(2.4>nc) 


<U)> 


Ka0/ac 


<) 


(3: ) 


where  I(a,/aQ,k)  is  a  function  obtained  by  integration.  Similarly,  the  rela¬ 
tive  toughening  is  expressed  through  the  parameter. 


r  =  2<u)>(aV/yt)2  Vf(2.4+nc)  I(a0/ac,k)  >33 ) 

Inserting  the  known  values  for  the  transformation  strains  and  the 
critical  variant  density  for  ZK^  particles  in  a  ZrC?  matrix  into  Eqn.  (35), 
the  predicted  toughening  is  plotted  as  a  function  of  particle  size,  in  Fig.  6, 
for  several  values  of  the  shape  parameter  k.  It  is  observed  that  a  peak 

A 

toughening,  Fy ,  occurs  for  each  k,  coincident  with  a  oarticular  value  of 
the  ratio  of  the  scale  parameter  of  the  distribution  to  the  critical  size 

A 

a. 

5.2  Comparison  With  Experiment 

The  predictions  of  the  above  model  can  be  comoared  with  toughness 
data  obtained  for  PSZ^and  for  Al^O^/Z^^.  For  this  purpose,  the  solu¬ 
tions  pertinent  to  a  size  distribution  of  particles  are  used.  Emphasis 

is  placed  on  solutions  in  the  range,  5<k<10,  as  representati ve  scale 

24 

parameters  for  typical  microstructural  entities  (discreteness  of  the 
transformation  that  may  occur  in  the  presence  of  a  small  number  of 
variants  per  particle  is  neglected).  The  comparisons  with  available 
data  can  only  be  superficial,  because  of  the  ancilliary  information  needed  to 


insert  into  the  models  (particularly  the  particle  sizes)  is  not  available. 
More  complete  comparisons  will  be  presented  in  a  subsequent  publication^. 
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For  the  PS7  material  of  Porter  and  Heuer\  the  primary  variable 
that  influences  the  toughening  is  the  particle  size,  which  is  consistenly 
varied  by  extending  the  annealing  period.  Changes  in  volume  fraction  are 
also  occurring  during  the  precipitation  process,  but  these  changes  are 
assumed  to  be  secondary,  within  the  range  of  the  optimum  toughness. 
Particle  coarsening  following  precipitation  often  occurs  as  a  function 
of  time,  t,  in  accord  with  a  t  dependence  (the  interface  control 
result).  If  this  behaviour  is  assumed  to  pertain  in  the  Porter  and 
Heuer  experiments^,  a  comparison  between  the  toughening  measurements 
and  the  predictions  of  the  present  model  are  shown  in  Fig.  7  for  k  =  6 
(the  Schoenlein  result  is  used  for  the  optimum  toughness,  since  this 
result  excludes  any  contribution  from  surface  compressive  stresses). 

The  relatively  good  correlation  is  encouraging. 

7 

The  results  of  Lange  for  the  A^O^/ZrO^  system  can  also  be  examined 
in  the  context  of  the  present  model.  The  principal  variable  in  Lange's 
experiments  is  the  volume  fraction.  However,  there  may  also  be  a 

systematic  increase  in  particle  size  as  the  volume  fraction  is  increased. 
If  this  latter  effect  is  neglected,  the  experimental  results  can  be 
compared  with  the  model  predictions  for  volume  fractions,  "<  0.3 
(i.e.,  isolated  particles):  another  model  would  be  required  to 
establish  the  behaviour  expected  for  a  continuous  ZrO^  phase.  Compari¬ 
sons  in  this  range,  illustrated  in  Fig.  8,  indicate  that  the  observed 
trends  are  adequately  predicted  by  the  model.  It  is  reemphasized,  how¬ 
ever  that  the  comparison  is  superficial  in  the  absence  of  size  dis¬ 
tribution  information. 


-22- 


6.  CONCLUSION 

The  energy  changes  that  accompany  the  complete  martensitic  trans¬ 
formation  of  an  isolated  particle  have  been  examined.  It  has  been 
demonstrated  that  one  of  the  dominant  energy  terms  depends  primarily 
upon  the  surface  area  of  the  particle,  whereas  the  other  principal 
terms  depend  upon  the  particle  volume.  A  particle  size  dependence 
of  the  transformation  thus  naturally  emerges  from  energy  change  consi¬ 
derations.  The  energy  term  that  depends  upon  the  particle  surface 
area  is  the  strain  energy  that  derives  from  the  twinned  or  variant 
structure  of  the  martensite.  The  areal  character  of  this  dependence 
arises  because  the  twinning  stresses  are  localized  near  the  particle 
interface. 

Lower  bound  estimates  of  the  critical  particle  size  for  stress 
free  transformation  have  been  obtained  from  the  energy  changes  asso¬ 
ciated  with  complete  transformation.  The  analysis  predicts  a  specific 
number  of  twins  or  variants  within  particles  at  the  critical  size. 
Predictions  performed  for  two  modes  of  transformation  in  Zr02  par¬ 
ticles  provide  a  consistent  description  of  the  observed  size 
effects. 

The  size  effect  has  been  incorporated  into  analyses  of  the  toughen¬ 
ing  induced  by  crack  tip  transformations.  The  toughening  that  can  be 
realized  has  been  shown  to  depend  on  the  particle  size  distribution. 
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as  well  as  the  chemical  tree  enerqy  change  and  the  oertinent  transfor¬ 
mation  strains.  The  choice  of  a  typical  size  distribution  permits  the 
observed  toughening  trends  for  partially  stabilized  Zrf^  and  for  Al^O^/ 
ZrOg  alloys  to  be  correlated  with  the  predictions  of  i:he  model. 
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FIGURE  CAPTIONS 


Fig.  1. 

Fig.  2. 

Fig.  3. 

Fig.  4. 

Fig.  5. 

Fig,  6. 

Fig,  7. 

Fig.  8. 


A  schematic  indicating  the  constrained  and  unconstrained  shapes 
of  monoclinic  ZrC^  particles  containing  variants. 

A  schematic  indicating  the  Eshelby  method  of  calculation  applied 
to  a  twinned  (or  variant)  martensite. 

The  normalized  shear  stress  plotted  as  a  function  of  distance 
from  the  interface. 

A  plot  of  the  normalized  strain  energy  as  a  function  of  the  number 
of  twins  (Eqn,  7)  compared  with  the  discrete  results  of  Kato  et 


Transmission  electron  micrograph  of  transform, ed  ZrOg  particles 
in  an  AI^O-^  ^  matrix. 

The  effects  of  a  size  distribution  on  the  toughness  of  Zr02  pre¬ 
dicted  by  the  analysis. 

A  comparison  of  the  toughness  data  obtained  by  Porter  and  Heuer1 
and  by  Schoenlein-*  with  the  behavior  predicted  by  the  model. 

A  comparison  of  the  toughness  data  for  Al^O^/ZrC^  obtained  by 
Lange^,  with  the  predictions  of  the  model. 
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ABSTRACT 

A  thermal  shock  test  has  been  designed  which  permits  the  thermal 
fracture  resistance  and  the  mechanical  strength  of  brittle  materials  to 
be  quantitatively  correlated.  Thermal  shock  results  for  two  materials, 
A1 2O3  and  SiC,  have  been  accurately  predicted  from  biaxial  strength 
measurements  and  a  transient  thermal  stress  analysis  (performed  using  a 
finite  element  method).  General  implications  for  the  prediction  of 
thermal  shock  resistance,  with  special  reference  to  ceramic  components, 
are  discussed. 


1 .  INTRODUCTION 


Thermal  transients  are  a  ubiquitous  source  of  fracture  in  ceramic 
components.  A  thorough  comprehension  of  the  variables  that  dictate 
thermal  fracture  are  thus  essential  if  the  reliable  performance  of  cer¬ 
amics  in  fracture  critical  applications  is  to  be  achieved.  A  logical 
approach  for  evaluating  the  susceptibility  of  a  ceramic  component  to 
failure  during  a  thermal  transient  is  based  on  a  parity  between  thermal 
shock  resistance  and  mechanical  strength.*  However,  this  approach 
has  not  been  uniquely  validated;  in  fact,  certain  observations  on 

ceramics  appear  to  be  superficially  inconsistent  with  the  approach^  K 

(2  3) 

Two  comparative  studies  have  been  conducted  *  '  that  afford  some  con- 

fidence  in  the  pertinence  of  the  mechanical  strength  for  thermal  shock 

(2) 

prediction.  Manson  and  Smith'  presented  statistical  relations  between 
mechanical  strength  (measured  in  flexure)  and  thermal  shock  resistance, 
and  demonstrated  that  the  statistical  shape  parameters  deduced  from 
both  thermal  shock  and  flexural  strength  tests  (performed  on  steatite) 
were  similar.  However,  they  did  not  attempt  a  prediction  of  the  absolute 
thermal  failure  condition  from  their  mechanical  strength  data.  More  re¬ 
cently,  the  thermal  failure  of  a  precracked  polymer^  has  been  predicted 
from  independent  measurements  of  the  crack  dimensions  and  the  critical 
stress  intensity  factor,  by  employing  a  numerical  thermal  stress  analysis. 
In  this  study,  an  approximate  choice  of  the  heat  transfer  coefficient 
was  used  to  provide  reasonable  predictions  of  the  failure  condition. 

The  intent  of  the  present  paper  is  to  provide  a  fully  quantitative  pre¬ 
diction  of  thermal  failure.  This  is  achieved  by  developing  a  suitable 

thermal  shock  test  and  by  devising  a  method  for  the  accurate  calibration 
_  • 

The  mechanical  strength  refers  to  the  stress  level  when  the  stress 
intensity  factor  at  the  fracture  initiating  flaw  attains  the  critical 
value,  «T£. 


-2- 


of  the  heat  transfer  coefficient;  a  method  which  involves  the  conduct  of 
failure  tests  on  a  model  brittle  material  (e.g.  A1 203 ) .  The  thermal 
failure  of  another  ceramic  material  (SiC)  is  then  predicted  from  mechanical 
strength  data  by  adopting  the  previously  calibrated  heat  transfer  coeffi¬ 
cient. 

A  quantitative,  but  simple  thermal  shock  test  for  evaluating  and  com¬ 
paring  ceramic  materials  (especially  those  suitable  for  advanced  applica¬ 
tion  such  as  gas  turbine  engines,  heat  exchangers  and  solar  collectors) 

is  not  presently  available.  A  procedure  commonly  adopted  for  thermal 

(4) 

shock  evaluation  is  based  upon  a  water  quench  test'  the  results  of  which 

can  be  incorporated  into  established  theories  of  fracture  initiation  and 
(5) 

crack  propagation'  .  Under  the  high  heat  transfer  conditions  that  prevail 
during  a  water  quench,  heat  transfer  and  size  effects  have,  in  most  previous 
studies,  been  assumed  to  be  of  negligible  importance.  The  peak  thermal 

A, 

stress,  o,  has  then  been  calculated  by  applying  the  simple  equation; 

:  -  ^  ■  (i) 

where  E  is  Young's  modulus,  a  is  the  coefficient  of  thermal  expansion, 

AT  Is  the  temperature  differential  across  the  specimen  and  v  is  Poisson's 
ratio.  By  setting  the  peak  tension  equal  to  the  equivalent  fracture  stress 
of  the  material,  an  estimate  of  the  maximum  temperature  differential,  aTc  . 
that  can  be  sustained  by  a  material  prior  to  severe  strength  degradation 
is  then  provided.  However,  the  water  quench  test  has  several  problems 
which  limit  its  utility  as  a  quantitative  test  for  evaluating  thermal  shock 
resistance. 

Firstly,  it  has  recently  been  shown  that  the  above  relation  is  only 


independent  of  sample  size  when  the  samples  are  extremely  large^.  Thus, 

ATc  shows  a  strong  size  dependence  in  many  experiments,  rendering  material 
ranking  somewhat  questionable.  Secondly,  the  heat  transfer  rates  encountered 
in  typical  ceramic  applications  (cited  above)  are  appreciably  smaller  than 
those  enforced  by  a  water  quench.  The  magnitudes  of  the  thermal  stresses 
expected  in  these  applications  are  thus  smaller  than  anticipated  by  Eq.(l), 
and  described  by  the  general  result: 

where  h  is  the  heat  transfer  coefficient,  k  is  the  thermal  conduc¬ 
tivity,  c  is  the  specific  heat,  o  is  the  density  and  r  is  a  speci¬ 
men  dimension;  the  quantity  hr/k  is  known  as  the  Biot  modulus, 
kt/ccr2  is  the  Fourier  number  and  F  and  G  are  functions^ .  Hence, 
many  more  material  parameters  influence  typical  thermal  shock  failures 
than  the  elementary  water  quench  test  analysis  admits.  Thirdly,  complex 
variations  in  the  heat  transfer  coefficient  with  test  conditions  present 
appreciable  interpretation  difficulties.  For  example,  in  one  study  per¬ 
formed  on  a  metal  wire  in  water,  h  was  found  to  vary  by  >104  between 
20  and  400°C,  due  to  boiling  at  the  metal-water  interface^ .  Although 
the  heat  transfer  variability  can  be  essentially  eliminated  by  quenching 
into  oil,  edge  and  corner  effects  cause  futher  difficulties/9^ 

A  thermal  shock  resistance  testing  and  evaluation  procedure  which 
alleviates  the  difficulties  associated  with  variable  and  very  high  heat 


transfer  rates  and  edge  effects  is  proposed  in  this  paper.  The  rationale 
and  the  experimental  procedure  are  firstly  presented  and  discussed.  A 
numerical  stress  analysis,  essential  to  the  detailed  interpretation  of 
the  experimental  results,  is  then  described.  Finally,  some  important 
implications  for  fracture  prediction  under  thermal  transients  are  dis¬ 
cussed. 
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2.  EXPERIMENTAL 

2.1  Rationale 

The  mechanical  fracture  of  a  ceramic  exhibits  appreciable  statis¬ 
tical  variability.  The  statistical  nature  of  the  fracture  also  results 
in  a  fracture  criticality  that  depends  upon  the  stressed  volume  and  the 
stress  state.  These  statistically  derived  effects  pose  considerable 
limitations  upon  the  direct  comparison  of  mechanical  fracture  and  ther¬ 
mal  shock.  This  difficulty  is  alleviated  if  the  comparisons  are  con¬ 
ducted  on  samples  with  controlled  pre-cracks,  which  extend  to  failure 
with  minimal  statistical  variability.  The  principal  experimental 
results  are  those  obtained  on  pre-cracked  specimens.  The  pre-cracks 
are  introduced  at  the  locations  of  peak  tension,  to  provide  a  basis  for 
subsequent  comparisons  of  normally  prepared  surfaces  at  the  location  of 
maximum  fracture  probability. 

The  experiments  are  designed  to  provide  a  comparison  of  the  mech¬ 
anical  strength  and  the  critical  temperature  for  thermal  fracture  ini¬ 
tiation,  at  a  prescribed  location  and  under  equivalent  stress  states 
(equi-biaxial  tension).  This  is  achieved  by  employing  disc  samples, 
in  which  the  edge  stresses  (and  hence,  undesirable  edge  initiated 
failure)  can  be  essentially  eliminated.  The  mechanical  tests  are  con¬ 
ducted  in  a  biaxial  flexure  mode;  while  the  thermal  tests  utilize  a 
disc  at  uniform  initial  temperature,  cooled  by  a  fluid  jet  impinging 
upon  the  disc  center. 


2.2  Procedures 

Two  materials  were  selected  for  the  thermal  shock  study:  a  99% 
dense  sintered  a-SiC  (~  7  ym  grain  size)  and  a  fully  dense  slip  cast 
A1 2O3  (~  15  ym  grain  size).  The  samples,  5  cm  diameter  discs,  0.25  cm 
thick,  were  ground  and  polished  on  one  surface.  Subsequently,  the  cen¬ 
ters  of  the  polished  faces  were  precracked  using  a  Knoop  indenter  with 
loads  of  15  to  33N.  Residual  stresses  caused  by  the  indentation  pro¬ 
cess  were  removed  by  fine  grinding. 

Each  sample  (both  precracked  and  non-precracked)  was  individually 
tested  in  the  thermal  stress  apparatus  shown  in  Fig.  1.  The  sample  was 
mounted  horizontally  on  fibrous  insulation  (with  two  thermocouples 
placed  against  the  surface  of  the  specimen)  and  heated  in  a  MoSi2  resis¬ 
tance  tube  furnace.  The  sample  was  allowed  to  equilibrate,  at  which 
time  is  was  subjected  to  a  rapid  temperature  change,  through  the  use  of 
high  velocity  air.  The  air  was  channeled  onto  the  disc  center,  using  a 
0.32  cm  diameter  silica  tube,  at  a  velocity  of  -  100  ms”1.  After  the 
quench,  the  sample  was  examined  for  crack  extension.  If  the  crack  did 
not  extend,  the  temperature  differential  between  the  sample  and  the  air 
jet  was  incrementally  enhanced  until  crack  extension  was  detected. 

This  critical  temperature  differential,  ATc>  was  recorded. 

Fractographic  analysis  was  performed  to  insure  that  fracture  originated 
at  the  precracks  in  the  indented  samples.  Typical  specimens  failed  by 
thermal  shock  are  shown  in  Fig.  2. 

Identical  samples  were  tested  in  biaxial  flexure  at  room  temperature 
using  the  apparatus  desioned  by  Uachtnan  et  al.^111.  A  5  cm  disc  was 


suoported  on  three  equally  SDaced  balls  (concentric  with  the  load)  and 
loaded  with  a  flat  piston  at  the  rate  of  "^MPas”1. 


3.  RESULTS 


The  variations  in  the  critical  imposed  temperature  differential 
with  indent  load  obtained  for  Al^O^  are  shown  in  Fig.  3.  The  results 
were  obtained  by  employing  successive  10°C  temperature  increments,  prior 
to  failure.  Fracture  in  the  precracked  specimens  occurred  within  a 
relatively  narrow  AT  range  for  each  precrack  size.  The  critical  tem¬ 
perature  differential  increased  as  the  precrack  size  diminished,  ranging 
between  480°  and  580°C.  The  specimens  with  as-machined  surfaces  failed 
over  a  wider  range  of  temperature:  a  phenomenon  related  to  the  flaw 
size  distribution.  The  approximate  failure  times  (obtained  from  the 
time  when  the  temperature  at  the  lower  thermocouple  registered  a  rapid 
decrease)  were  ~  5  s. 

Results  for  sintered  a-SiC  were  more  difficult  to  obtain,  because 
this  material  is  subject  to  time  dependent  strengthening  (by  oxidation 
or  by  surface  diffusion)  within  the  temperature  range  required  to  induce 
thermal  fracture,  >  1100°C.  It  was  elected,  therefore,  to  obtain  upper 
and  lower  bound  values  for  aTc  by  performing  only  a  single  thermal 
shock  test  on  each  precracked  specimen  (in  contrast  to  the  sequential 
testing  employed  with  Al^O^).  Firstly,  a  preliminary  upper  bound  aTc 
was  established  for  each  specific  precrack  size,  by  inducing  thermal 
failure.  Subsequent  tests  on  specimens  with  the  same  precrack  size 
were  then  performed  at  sequentially  lower  temperatures.  The  lowest  tem¬ 
perature  at  which  failure  occurred  was  designated  the  upper  bound, 
while  the  highest  survival  temperature  was  denoted  the  lower  bound.  The 
results  obtained  using  this  procedure  are  plotted  in  Fig.  3.  It  is 
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noted  that  samples  without  precracks  could  not  be  thermally  shockec 

within  the  temperature  limitations  of  the  present  system  (M500°C). 

The  biaxial  flexure  data  for  both  materials  are  plotted  as  a  func- 

(12 ) 

tion  of  the  indentation  load  in  Fig.  4.  The  theoretical  slope  of 

-1/3  is  superimposed  to  emphasize  the  trend.  It  is  noted  that  SiC 

exhibits  lower  indentation  strengths  than  A^O^;  a  trend  *hat  probably 

(12) 

reflects  the  fracture  toughness  characteristics  of  the  two  materials 
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4.  STRESS  ANALYSIS 

The  thermal  stresses  that  develop  within  the  disc  following  impinge¬ 
ment  of  the  air  jet  were  determined  using  axisymmetric  finite  element 
schemes.  The  temperature  distributions  were  established  by  direct 
application  of  the  program  Do/13/  The  temperatures  were  then  used  to 
ascertain  the  thermal  stress,  by  an  adaptation  of  the  program  SOLSAp/^ 

The  calculations  were  conducted  by  allowing  heat  transfer  through 
the  upper  surface  of  the  disc,  to  the  air  stream.  Radiation  from  the 
lower  surface  of  the  disc  into  the  cooled  regions  was  also  permitted. 

The  heat  transfer  coefficient  was  assumed  to  be  uniform  over  the  diam¬ 
eter  of  the  air  jet,  with  a  magnitude  proportional  to  that  expected  for 
a  fluid  stream  flowing  against  a  thin  pi  ate  Z1 5  ^ 


where  k^  and  are  the  thermal  conductivity  and  kinematic  viscosity 
of  the  fluid  at  the  film  temperature,  d  is  the  jet  diameter,  v*,  is 
the  fluid  velocity  and  P  is  the  Prandtl  number;  X  is  a  proportion¬ 
ality  constant,  to  be  determined  by  a  calibration  procedure,  described 
below. 

Heat  transfer  was  also  considered  to  occur  as  a  consequence  of  the 
jetting  of  the  fluid  over  the  surface.  The  heat  transfer  coefficient 
for  boundary  layer  flow  over  a  flat  surface  was  selected, 


and 
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all  owed  to  diminish  inversely  with  distance  from  the  jet,  to  account  for 
flow  attenuation.  However,  preliminary  calculations  indicated  that  rea¬ 
sonable  choices  for  this  heat  flow  exerted  a  relatively  minor  influence 
on  the  thermal  stress  developed  at  the  disc  center  (the  location  of 
present  interest). 

Four  axi symmetric  element  groups  were  employed  in  the  study:  an 
interior  two  dimensional  element  (group  1),  a  surface  convection  element 
with  constant  heat  transfer  (group  2),  a  surface  convection  element  with 
h  inversely  proportional  to  radial  distance  (group  3)  and  a  uniform 
radiation  element  (group  4).  The  assignment  of  these  elements  is  illus¬ 
trated  in  Fig.  5.  Four  different  meshes  were  selected  (Fig.  6)  in  order 
to  ascertain  the  influence  of  the  mesh  definition  and  to  establish  con¬ 
vergence  conditions.  It  was  firstly  established  that  consistent  results 
emerged  provided  that  the  mesh  distribution  was  relatively  uniform  both 
in  the  central  region  (of  high  heat  transfer)  and  across  the  adjacent 
discontinuity  in  h.  Thereafter,  it  was  determined  that  a  grid  with 
216  elements  (grid  (c)  in  Fig.  6)  was  the  preferred  choice  for  subsequent 
analysis.  This  grid  yielded  temperature  distributions  within  IS  and 
stresses  within  6 %  of  those  calculated  with  the  finest  grid  (559  elements) 
an  accuracy  deemed  sufficient  for  present  purposes. 

Temperature  and  stress  calculations  were  conducted  by  admitting 
temperature  dependent  material  properties:  thermal  conductivity,  speci¬ 
fic  heat,  thermal  expansion  coefficient  (Fig.  7)  and  elastic  modulus 
(~400  GPa  for  both  materials).  Typical  temperature  distributions  within 
the  test  specimen,  determined  both  across  the  surface  and  through  the 


thickness,  are  plotted  in  Fig.  8.  It  is  noted  that  there  is  an  appre¬ 
ciable  axial  temperature  gradient.  Results  for  a  thin  disc  can  not, 
therefore,  be  expected  to  apply. 

The  time  dependence  of  the  maximum  stress,  at  the  center  surface 
location,  typically  exhibits  the  form  depicted  in  Fig.  9.  It  is 
observed  that  the  peak  stress  develops  at  a  time  approximately  coinci¬ 
dent  with  the  instant  when  the  temperature  at  the  disc  periphery  begins 
to  diminish. 

The  procedure  adopted  for  the  calibration  of  the  heat  transfer 
coefficient  involved  calculations  of  the  influence  of  h  on  the  peak 
stress  for  A^Oj  at  a  specific  temperature  differential  575°C,  pertinent 

to  a  15N  indentation  precrack.  The  value  of  the  heat  transfer  coeffi- 

★ 

cient,  h  ,  that  provided  a  peak  tensile  stress  coincident  with  the 

biaxial  fracture  strength,  for  precracks  of  equivalent  size  (290  MPa), 

was  then  considered  to  be  the  correct  value  of  h  for  that  film 

temperature.  A  unique  X  value  for  the  test  system  was  then  deduced 
★ 

by  inserting  h  into  Eq.  (3).  The  calibrated  value  of  X  was  employed 
for  all  subsequent  stress  calculations,  yielding  values  of  h  that  vary  with 
temperature,  through  the  temperature  dependence  of  k^  and  (for 
the  fluid  medium). 

★ 

Preliminary  credence  in  the  deduced  value  of  h  was  established 
by  comparing  the  measured  failure  times  with  estimates  derived  by 
analysis.  For  this  purpose,  it  is  noted  that,  since  10°C  temperature 


increments  were  used  to  determine  the  incidence  of  fracture,  the  experi¬ 
mental  estimate  of  ATC  could  exceed  the  actual  value  by  <  10°C. 

Stress  levels  developed  at  a  10°C  temperature  separation  were  thus  cal¬ 
culated,  and  a  lower  bound  on  the  failure  time  established  from  the  time 
when  the  stress  at  the  upper  temperature  attained  the  median  failure 
stress  (Fig.  9).  The  time  of  12  s  determined  using  this  procedure  is 
of  the  sane  order  as  the  measured  failure  times. 

A  series  of  calculations  were  conducted  for  aT  values  within  the 
vicinity  of  the  values  measured  for  each  of  the  two  materials.  These 
calculations  permitted  relationships  between  AT  and  the  peak  stress 

A 

a  to  be  deduced  for  each  material  within  the  pertinent  AT  range,  as 
plotted  in  Fig.  10.  Imposing  the  median  biaxial  failure  stress  for  each 
precrack  size  (Fig.  4)  then  allowed  the  critical  AT  for  each  thermal 
shock  test  to  be  predicted. 

The  trends  in  thermal  fracture  predicted  in  this  manner  (Fig.  3) 
are  remarkably  consistent  with  the  measured  behavior.  Of  particular 
merit  are  the  accurate  predictions  obtained  for  SiC,  which  fractures  in 
a  very  different  AT  range  than  the  calibration  material  (Al^),  by 
virtue  both  of  appreciable  differences  in  k  and  a  and  of  significant 
variations  in  Kc  and  h.  It  is  also  noted  that  superior  predictions  are 
obtained  by  excluding  radiation  from  the  lower  surface  of  the  specimen. 
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5.  IMPLICATIONS 

The  close  prediction  of  thermal  shock  failure  from  (surface  crack 
dominated)  biaxial  strength  measurements  strongly  supports  the  notion 
that  thermal  fracture  can  be  uniquely  related  to  the  stress  field  pro¬ 
duced  by  the  temperature  transient.  Many  thermal  shock  failures  origi¬ 
nate  from  surface  cracks  (because  the  thermal  stress  gradient  encourages 
failure  from  near-surface  located  defects).  The  present  validation  of 
the  thermal  failure  process  for  this  important  defect  type  should  thus 
be  of  substantial  merit  for  the  prediction  of  thermal  shock  failure  in 
ceramic  components. 

The  calibration  of  the  heat  transfer  conditions  for  the  present 
test  arrangement  permits  the  apparatus  to  be  used  for  the  quantitative 
determination  of  thermal  shock  resistance  for  a  range  of  ceramics. 

The  comparison  can  be  conducted  on  surfaces  prepared  in  a  manner  analo¬ 
gous  to  that  used  for  actual  components,  thereby  averting  the  extraneous 
influence  of  specimen  edges.  (Also,  the  relative  temperature  invariance 
of  the  heat  transfer  coefficient  permits  a  direct  ranking  of  thermal 
shock  resistance.)  Such  quantitative  thermal  shock  tests  would  allow 
both  the  elucidation  of  the  statistical  character  of  fracture  and  would 
permit  the  magnitude  of  thermal  (e.g.,  interrupted  heat  flow)  sources 
of  stress  intensification  pertinent  to  fracture  from  fabrication  defects 
(voids,  inclusions,  etc.)  to  be  evaluated. 

For  intrinsically  thermal  shock  resistant  materials  (SiC,  Si-jN^), 
more  severe  heat  transfer  conditions  than  those  presently  employed 
would  need  to  be  devised,  in  order  to  induce  fracture  within  the  temperature 
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capabilities  of  the  test  system.  This  might  be  achieved  by  increasing 
the  specimen  dimensions  or,  more  expediently,  by  increasing  the  conductivity 
and  decreasing  the  viscosity  of  the  fluid  (see  Eq.  (3)).  An  appropriate 
choice  of  fluid  might  be  helium. 

Finally,  it  is  observed  that  the  heat  transfer  calibration  pro¬ 
cedure  adopted  in  the  present  study  could  provide  an  unequivocal  method 
for  heat  tr.-cfer  calibration  in  actual  systems.  Specifically,  a  small 
precrack  could  be  placed  at  the  location  of  maximum  stress.  Then,  the 
time  at  which  fracture  initiates  from  the  precrack  could  be  measured, 
during  a  test  run.  Thereafter,  a  combination  of  the  failure  time  with 
the  magnitude  of  the  failure  stress  at  the  crack  would  permit  the  heat 
transfer  coefficient  to  be  determined  by  comparison  with  a  series  of 
thermal  stress  calculations. 
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FIGURE  CAPTIONS 


Fig.  1  A  schematic  of  the  thermal  stress  test  apparatus. 

Fig.  2  Thermally  shocked  discs 

(a)  a  substantially  supercritical  thermal  shock  that 
activates  many  surface  cracks 

(b)  a  slightly  supercritical  thermal  shock  activating 
a  central  precrack. 

To  enhance  the  cracks  the  samples  were  immersed  in  a  dye 
penetrant  and  illuminated  under  ultra-violet  light. 

Fig.  3  The  influence  of  indentation  load  upon  the  critical  tempera¬ 
ture  for  thermal  fracture  of  A1 2O3  and  SiC.  Also  shown  are 
the  critical  temperatures  predicted  by  the  stress  analysis 
from  the  biaxial  strength  results. 

Fig.  4  The  influence  of  indentation  load  on  the  biaxial  flexure 
strength  of  A1 2O3  and  SiC. 

Fig.  5  The  assignment  of  the  four  types  of  element  used  in  the 
finite  element  analysis. 

Fig.  6  The  four  meshes  employed  in  the  finite  element  calculations. 
Mesh  c  was  employed  for  the  majority  of  calculations. 

Fig.  7  Temperature  dependent  material  properties  of  AI2O3  and  SiC 

(a)  thermal  conductivity 

(b)  specific  heat 

(c)  thermal  expansion 

Fig.  8  The  axial  and  radial  temperature  distributions  in  an  A1 2O3 
sample  subject  to  a  temperature  differential  of  575°C, 
obtained  at  two  times,  8s  and  32s. 

Fig.  9  The  variation  of  the  maximum  tensile  stress  with  time  for 
A1 2O3  subject  to  several  values  of  IT.  Also  shown  is  the 
construction  for  estimating  the  lower  bound  failure  time. 

Fig.  10  The  variation  of  the  peak  tensile  stress  with  AT  for. 

(a)  AT 2O3 ,  (b)  SiC.  Essentially  identical  results  obtain 
for  A1 2O3  both  with  and  without  radiation;  but  a  significant 
disparity  exists  for  SiC,  as  illustrated. 
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ABSTRACT 

An  approximate  analysis  is  presented  of  the  toughening  induced 
by  the  incorporation  of  second  phase  particles  subject  to  microcracking. 
The  toughening  is  demonstrated  to  become  appreciable  for  a  narrow  size 
distribution  of  particles  of  appropriate  size.  The  magnitude  of  the 
toughening  is  determined  by  the  mismatch  strain  (due  to  thermal  contrac¬ 
tion  incompatability)  and  the  microfracture  resistance  of  the  particle/ 
matrix  interface.  Implications  for  designing  optimally  tough  ceramics 
are  presented. 
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1.  INTRODUCTION 

The  incidence  of  microcracking  to  form  a  process  zone  around 

macrocrack  tips  is  now  a  well-substantiated  phenomenon,  for  certain 

12  3 

classes  of  brittle  solids.  *  *  A  comprehensive  analysis  of  this  problem 
has  been  elusive,  however,  because  several  complex,  interactive  effects 
accompany  the  microcracking  process.  Firstly,  the  microcracking 

4 

increases  the  compliance  of  the  process  zone  and  thereby,  modifies  the 

crack  tip  stress  field  from  the  initial  linear  elastic  field  (Fig.  1). 

c 

Some  simple  approximations  have  been  used,  but  these  are  not  fundamen¬ 
tally  satisfying.  Secondly,  the  microcrack  density  will  tend  to  decrease 
with  distance  from  the  crack  tip,®  resulting  in  a  diffuse  process  zone 
(Fig.  1),  and  a  tendency  for  interaction  with  the  boundaries  of  test 
samples.  Thirdly,  the  criterion  for  coalescence  of  microcracks  is  un¬ 
certain.  Coalescence  with  the  primary  crack  is  especially  significant 
because  this  event  determines  the  crack  extension  condition®  and  hence, 
dictates  the  toughness  of  the  material. 

A  detailed  study  of  the  microcrack  problem  will  inevitably  in¬ 
volve  a  computer  simulation.  A  first  rational  attempt  at  simulation  has 
recently  been  reported  by  Hoagland  and  Embury,7  using  an  image  force 

O 

(based  on  a  Green's  function  derived  by  Hirth  et  a!)  to  account  for 
compliance  effects.  Much  additional  study,  based  on  this  approach,  will 
be  needed  to  develop  a  comprehension  of  effects  of  microcracking  on 
crack  extension  resistance  curves  and  on  toughness.  This  evaluation 
would  be  greatly  facilitated  by  the  availability  of  complementary  approx¬ 
imate  analytic  solutions  for  specific  microcrack  problems.  Such  solutions 
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would  also  provide  invaluable  guidance  for  research  studies  on  the  micro- 

structural  design  of  high  toughness  ceramics.  One  such  solution  is 

described  in  the  present  paper. 

It  must  be  recognized  at  the  outset  that  stable  microcracking 

in  brittle  materials  generally  results  from  large  localized  residual 

stresses  that  develop  due  either  to  thermal  contraction  anisotropy  or  to 

9  10- 

a  phase  transformation.  *  The  tensile  residual  stresses  superimpose 
on  the  applied  crack  tip  stresses  to  initiate  microcracking.  A 
knowledge  of  the  residual  stresses  is  thus  an  important  constituent 
of  microcrack  analyses. 

A  microcrack  system  of  appreciable  practical  significance  (that 
is  also  amenable  to  approximate  analytic  treatment)  is  a  two-phase  sys¬ 
tem  containing  second  phase  particles  with  a  larger  thermal  contraction 

1 1 

coefficient  than  the  host.  Such  particles  tend  to  microcrack  circum- 
12 

ferentially  and  thus,  are  intrinsically  stable  (the  microcracks  being 

confined  to  the  immediate  vicinity  of  the  particles).  This  microcrack 

problem  can  be  addressed  using  a  simple  thermodynamic  approach,  anal- 

13 

ogous  to  that  recently  described  for  martensitic  toughening,  if  appeal 

12 

is  made  to  two  significant  observations.  Firstly,  Davidge  and  Green 

noted  that,  for  this  type  of  particle,  the  onset  of  microcracking  can  be 

adequately  described  by  the  condition  that  the  total  elastic  energy 

exceeds  the  energy  needed  to  create  the  crack  surface.  Some  rationale 

14 

for  this  result  has  recently  been  presented  by  Ito,  et  al. 

Secondly,  the  computer  simulations  of  Hoagland  and  Embury^  have  indi¬ 
cated  that  the  extent  of  the  microcrack  zone  is  not  appreciably 
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different  from  that  predicted  by  neglecting  compliance  effects. 

Hence,  an  approximate  zone  size  can  be  estimated  from  the  linear 
elastic  field  equations.  These  observations  permit  the  zone  size 
problem  to  be  treated  by  evaluating  the  distance  from  the  crack  tip, 
r  ,  at  which  a  single  particle  satisfies  the  requirement  that  the  change 
in  thermodynamic  potential  of  the  system  (before  and  after  microcracking), 
becomes  incrementally  negative.  The  toughening  can  then  be  deduced  from 
the  process  zone  size. 

The  sequence  of  analysis  in  the  present  paper  consists  of  a  pre¬ 
sentation  of  the  governing  thermodynamic  relations,  followed  by  an  esti¬ 
mation  of  the  process  zone  size.  Thereafter,  the  toughening  increment 
is  evaluated  and  the  implications  for  microstructural  design  are  dis¬ 
cussed. 

2.  THE  GOVERNING  THERMODYNAMIC  RELATIONS 
2.1 .  The  Elastic  Energy 

The  increase  in  strain  energy  A4>0  of  a  system  containing  a 
single  particle  subject  to  a  dilational  transformation  strain  eT  is;^ 

*  -  (4)  (,) 

where  Vp  is  the  particle  volume  and  p*  is  the  stress  within  the  particle. 
The  transformation  strain  pertinent  to  the  present  problem  is  the  thermal 
expansion  mismatch  strain;  a  strain  that  is  essentially  dissipated  fol¬ 
lowing  separation  of  the  particle  from  the  matrix.  Hence; - 


eT  =  -3ActAT 
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where  Aa  is  the  thermal  expansion  mismatch,  assumed  to  be  isotropic 
for  both  matrix  and  inclusion  for  this  analysis,  and  AT  is  the  cooling 
range.  For  an  ellipsoidal  inclusion,  p1  is  constant,13  andeqn.  (1)  becomes: 


A4>0  = 


T  T 

Zve  p  abc 


where  a,  b  and  c  are  the  semi-axes  of  the  ellipsoid.  The  stress  within 
the  ellipsoid  is; 

I  2E  (e^/3) 

P1  -  -  - - S -  (3) 

0+vm)+26(l-2vp) 

where  E^  is  Young's  modulus  for  the  matrix  and  8  is  the  ratio  of  moduli 
between  the  matrix  and  the  particle  (^Ep).  Combining  eqns.  (2)  and  (3) 
and  simplifying  for  the  case  of  a  spherical  inclusion,  radius  b,  gives; 


A$o  * 


(4iT/3)b3(eT)  Em 
3[(Hvm)+26(l-2v  )]" 


which  for  vm  %  Vp  »  0.2  reduces  to; 


A<*>0  as 


(4m/3)b3(eT)  En 
3.60+5) 


The  Increment  inelastic  energy  (that  derives  from  particle  di¬ 
lation}  is  modified  in  the  presence  of  a  uniform  applied  stress  to,1 

a.  _  /4tt\.3,T,2  f  pA  0,3Em  1 


/4tt  \  3,  T,2  n"  0.3E 

**•  =  Wb  (e  1  £r, +  rrar 


HU 
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where  is  the  dilational  component  of  the  aoplied  stress.  Note  that 
an  applied  tension  (as  pertains  to  a  crack  tin  field)  augments  the  elastic 
energy  in  the  presence  of  a  particle  with  a  larger  expansion  coefficient 
than  the  matrix,  and  thus  encourages  separation.  However  an  additional 
elastic  energy  term  is  obtained  in  the  presence  of  an  applied  stress;  a 
term  that  derives  from  the  interaction  of  the. stress  field  with  the  elastic 
in homogeneity,  as  reflected  in  the  elastic  modulus  mismatch  between  the 
particle  and  matrix.  The  pertinent  interaction  energy  is1^; 


KrV 

m 


(4um  *  V 

*4um  +  V 


nA  A 
P  P 


15(u  -u  )(l-v)  p^.  pj. 

+  n  P  _ L3 _ 2J _  (7) 

2um^2up^4"5v^+um^7‘5v^  J 

•  /\ 

where  k  is  the  bulk  modulus,  u  the  shear  modulus  and  p  ^  is  the 
deviatoric  component  of  the  applied  stress.  For  a  Poisson's  ratio  of  0.2, 
eqn.  (7) reduces  to; 


^ 2  =  (~r)  b 


.)  (S-D 

fo.4  PV  +  2.4  'oA  'p*  1 

Em(9*l) 

L  ij  uj 

(8) 


When  a  particle  separates  from  the  matrix,  the  comoonent  of  the 
elastic  energy  associated  with  the  mismatch  strain  A<J>(  reduces  to  zero 
(Fig.  2).  Whereas,  the  energy  that  derives  from  the  elastic  inhomogeneity 
changes  to  that  for  a  void  (obtained  from  eqn.  (8)  with  ED  =  0)  given  by; 


0.4  pApA  *  2.4  'p*. 


(9) 
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Th  .  change  in  elastic  energy  of  a  spherical  particle  following  separa¬ 
tion  is  thus; 

r 

2 


4tt  .3  J  /  T  \ 
^elastic - r  b  Ve  > 


0A  +  °-3En, 

£7)  ™ 


-8 


[o.4pApA  ♦  2.4  'pAj'pAJ 


V6+1> 


(10) 


2.2  The  Surface  Energy 

The  increase  in  surface  energy  that  accompanies  circumferential  frac¬ 
ture  of  an  ellipsoidal  oarticle  is 


where 


¥  = 


9 


F(k,4>)  and  E(k,4')  are  elliptical  integrals  of  the  first  and  second  kinds 
respectively,  and  is  the  pertinent  fracture  energy.  Simplification 

of  eqn.(13)  can  be  achieved  for  a  spherical  particle,  radius  b. 


A*s  *  4*b2yint  (14) 

2.3  The  Change  in  Thermodynamic  Potential  of  Separation 
Neglecting  compliance  effects,  in  accord  with  our  original  premise, 
the  change  in  the  thermodynamic  potential  of  the  system  following  circum¬ 
ferential  particle  separation  is; 


A*  =  ^elastic  +  A*s 


05a) 


and  hence,  for  a  spherical  particle. 


Ap 


4irb 

3 


3Yint  “  b(eT)2 


°-3Em 

(1+S) 


0.4Bb  pApA+e'p^j 


A  lower  bound  estimate  of  the  applied  stress  needed  to  induce  microfacture 
(obtained  by  setting  Ap^O)  is  thus  embodied  in  the  relation-. 


0.48 


Wl+B)(-e‘) 


r  A  A,  c 1  A 

[p  p  +6  p.j 


A  t 

Pij] 


E  ( ~eT ) 
m  ' 


3Y, 


int 


bEfeV 

.  m 


0.3 

T+F 


(16) 


The  above  result  pertains  to  complete  separation  of  the  particle  from 

the  matrix.  However,  the  change  in  thermodynamic?  potential  for  partial 

separation  is  expected  to  be  closely  similar  to  that  given  by  eqn.  (15). 

The  relative  invariance  of  4*  with  circumferential  crack  length  2a'  can  be 

deduced  from  the  recent  strain  energy  release  calculations  of  Ito  et  al.^6^, 

a 

which  yield  a  deviation  in  Ap,  at  p  =0  given  by; 


<5Ad>/APelastic  =  -0.56  [a'/b  -  sin  (a'/b)cos  (a'/b) -(l/3)cos3(a'/b)] 

-  1.06  cos  (a'/b)  +  0  .88  (17) 


which  differs  from  zero  by  <4%  for  all  a'/b.  This  small  deviation  may  account 
for  the  success  of  the  simple  thermodynamic  criterion  for  microcrack  formation. 
3.  THE  PROCESS  ZONE  SIZE 

Prediction  of  the  process  zone  size  requires  the  introduction  of  a 
criterion  for  microcracking  in  the  crack  tip  field.  An 


upper  bound  estimate 
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is  presented  in  this  section  by  sinply  invokinn  the  requirement  that  particle 
separation  occurs  when  at  the  particle  reduces  to  zero,  by  virtue  of 
the  influence  of  the  crack  tip  stress  field  on  the  elastic  energy  of  separa¬ 
tion.  The  local  stress  used  to  assess  the  process  zone  dimension  is  assumed 
to  be  the  linear  elastic  crack  tip  field  and  must,  therefore,  be  regarded  as 
an  approximation. 

The  zone  size  r  ,  for  small  zones  relative  to  the  crack  length,  can 
be  estimated  by  inserting  into  eqn.  (16)  the  crack  tip  stresses,  given  for 
plane  stress  conditions  by^^; 


A  K  /sin9\ 
1  2  ) 


where  K  is  the  stress  intensity  factor,  (r,e)  are  the  particle  coordi¬ 
nates  with  respect  to  the  crack  tip  and  a  is  the  crack  length.  The  resultant 
quadratic  equation  for  rc  is  unwieldly  and  is  thus  not  presented  in  detail. 
The  zone  width  over  the  crack  surface  h  (fig.  3)  is  of  principal  present 
interest.  This  dimension  coincides  with  the  angle  of  which  a  line  parallel 
to  the  crack  becomes  tangent  to  the  process  zone,  and  occurs  at  0  ^  0.42tt. 

For  this  angle  and  for  6=1,  the  zone  width  is  given  by  (h«a); 


Firstly,  a  solution  does  not  exist  for  b  <  7y,_  ,/E  (e  )  ,  because  the 

I  1 1  w  III 


elastic  energy  increase  from  the  inhomogeneity  effect  exceeds  the  elastic 
energy  decrease  associated  with  the  mismatch  strain,  in  regions  close  to 


the  crack  tip  (where  e^  e^").  A  minimum  particle  size  b  .  is  thus 

min 


needed  to  induce  a  process  zone; 

nt 

Ti2 - 


bmin 


^  Yint 


hb'V 


(20) 


Secondly,  the  zone  size  tends  to  infinity  as  the  denominator  in  eqn.  (19) 
tends  to  zero.  This  coincides  with  the  condition  for  spontaneous  (stress 
free)  microcracking; 


20 

k  _  mt 

bc~  7TJ7F 


Em(eT 

m 


(21) 


a  result  that  has  been  frequently  quoted  in  prior  studies  of  microcrack- 
ing(12,17,18,). 


'  For  most  useful  ranges  of  particle  size,  bc  > b  >  >  b^. n ,  eqn.  (19)  can 
be  expressed  by  the  simplified  result  (h«a); 


1 


h  %  13  (  K  _ 

71  W7  [br/b  -  l]2 


(22) 


nr 


This  relation  allows  h+®  as  b-*-bc;  a  limit  that  is  not  consistent  with 
the  requirement  that  the  zone  height  be  small  with  respect  to  the  crack 
length.  A  convenient  approach  for  establishing  a  limit  on  the  zone  height 
is  to  recognize  that  some  microcracking  will  be  induced  at  the  level  of  the 
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remotely  applied  stress,  p^.  These  transformed  particles  cannot  partici¬ 
pate  in  the  toughening,  because  they  do  not  reflect  a  change  in  state  during 

/v 

a  crack  increment.  An  approximate  upper  bound  h  for  the  process  zone 
height  that  participates  in  the  toughening  is  obtained  from  the  zone  dimen¬ 
sion  at  which  the  stress  given  by  eqn.  (19)  reduces  to  the  level  of  the 
remote  stress; 


Equating  eqn.  (23)  with  eqn.  (22)  then  permits  the  definition  of  an  upper 
bound  particle  size  b  engaged  in  the  toughening; 

b/bc  -  [1  +  5  eje7]  =  1  +  5  (24) 


where  e^  is  the  remotely  applied  strain. 


4.  THE  TOUGHENING  INCREMENT 


The  increment  in  toughening  due  to  the  formation  of  the  microcrack 
process  zone  is  governed  by  the  detailed  stress  changes  that  occur  in  this 
zone  and  hence,  the  increase  in  the  applied  stress  needed  to  induce  crack 
extension.  However,  an  approximate  measure  of  the  toughening,  can  be  de¬ 
duced  by  examining  the  total  change  in  energy  of  the  system  following  a 
1 3 

crack  increment  .  Neglecting  the  energy  contained  in  stress  waves,  the 

13 

toughening  T-p  is  given  by  ; 


rT*'(A4>s  +  A4>0)  h 


(4/3)irb 


T  +  T  0 


(25) 
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where  rQ  is  the  intrinsic  resistance  to  crack  extension,  A<{>0  is  ne 
strain  energy  change  after  microcracking  in  the  absence  of  an  applied  stress 
(the  increment  in  process  zone  can  be  considered  to  be  added  at  the  essen- 
tially  stress-free  central  portion  of  the  crack),  and  V^.  is  the  volume 
fraction  of  particles.  The  most  serious  deficiency  of  eqn. (23)  is  the 
neglect  of  interaction  effects  between  particles.  The  result  should  only 
be  considered  pertinent,  therefore,  for  dilute  particle  concentrations. 
Estimation  of  the  importance  of  interaction  effects  will  require  a  numerical 
treatment.  Relaxations  that  occur  at  the  crack  surface  are  also  neglected; 
but  this  effect  should  be  small  for  process  zones  h  >3b. 

The  microcrack  toughening  can  be  readily  deduced  from  eqn. (25)  for 
the  ideal  uniform  particle  size  case;  and  this  calculation  is  conducted 
first.  Thereafter,  an  analysis  of  toughening  in  the  presence  of  a  size 
distribution  of  particles  is  attempted. 

4.1  Uniform  Particle  Size  Toughening 

Combining  eqns.  (5),  (14),  (22)  and  (25),  we  obtain,  for  8  =  1, 

(b<b); 


2YfK2 

rT  *  7  £mLbc/b-1J  +  "o 


However,  K  must  be  at  the  critical  level  Kc  during  a  crack  growth  in¬ 
crement;  hence,  for  linear  behavior. 


E  r 


T 


(27) 
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The  toughening  then  becomes 


(28a) 


where 


2Vj 


[bc/b-l] 


(28b) 


It  should  be  noted  that  the  dependence  of  the  zone  width  on  K  has 
caused  IV  ■+  ®  as  n  ■*  1.  However,  this  condition  will  never  be  realised 
for  fully  confined  process  zones, characterised  by  the  requirements  (noted 
above)  that  the  particles  within  the  zone  be  smaller  than  b.  The  implica¬ 
tions  of  eqn.  (28)  in  the  range  n  1  should  thus  be  discounted.  An  ap¬ 
proximate  form  of  eqn.  (28)  suitable  for  subsequent  analysis,  that  pertains 
over  a  useful  range  of  n,  is; 


[bc/b-l] 


(29) 


The  significance  of  the  intrinsic  toughness  I"0  in  eqn.  (28)  requires 
amplification.  The  presence  of  detached  or  cracked  particles'adjacent  to 
the  crack  tip  will  evidently  reduce  rQ  from  that  for  the  dense  matrix. 

The  magnitude  of  the  reduction  will  depend  on  the  detailed  distribution  of 
the  particles  and  the  specific  nature  of  the  interactions  between  the  micro¬ 
cracks  and  the  primary  crack.  A  first  order  estimate  based  on  the- reduction 

★ 

in  load  bearing  area  would  be: 


This  estimate  is  probably  pertinent  to  irregularly  shaped  particles  larger 
than  the  matrix  grain  size.  However,  it  should  be  noted  that  voids  of 
regular  shape  in  a  homogeneous  matrix  can  impede  crack  extension,  even 
thouoh  the  effective  load  bearing  area  is  reduced.  Alternate  relations 
between  r  and  T  should  thus  be  sought,  if  the  microstructure  suggests 
void  impediment  effects. 
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r_  -  r_(l- 

o  nr 


Vf) 


(30) 


where  r  is  the  matrix  touqhness.  The  trends  with  volume  fraction  nre- 
m 

dieted  by  this  result  are  plotted  in  Fig.  4,  for  various  particle  radii 
(b  <  b  )  and  for  8  =  1.  It  is  evident  from  the  figure  that  the  appreciable 

v» 

increases  in  toughness  occur  for  volume  fractions  up  to  ^0.3.  More  signifi¬ 
cantly,  however,  it  becomes  apparent  that  the  particle  size  must  be  relatively 
close  to  the  critical  value,  in  order  to  achieve  appreciable  toughening. 

This  imposes  a  significant  practical  constraint  upon  the  consistent  attain¬ 
ment  of  high  toughness. 

Finally,  it  is  noted  that  the  toughening  exoressed  by  eqn.(28)  is  in¬ 
dependent  of  particle  shape  (however  the  shape  does  have  an  appreciable 
influence  on  the  critical  particle  size). 


4.2  Particle  Size  Distribution  Effects 

The  introduction  of  a  distribution  of  oarticle  sizes  evidently  results 
in  a  degradation  of  the  maximum  toughening  anticipated  by  the  uniform  particle 
size  analysis.  The  specific  influence  of  a  size  distribution  can  be  demon¬ 
strated  to  coincide  w.ith  the  replacement  of  Vf  in  eqn.  (26)  with  the  equivalent 


distribution  function,  yielding  the  relation 


b^4('b)db 


(31) 


where  $(b)db  is  the  particle  size  distribution  function  and  Vf  is  now 
the  total  volume  fraction  of  particles.  The  major  contribution  to  the 
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toughening  derives  from  particles  in  the  size  range  just  below  the  critical 
size  'B.  Hence,  since  high  toughness  requires  an  appreciable  fraction  of 

/  Av 

the  population  to  be  uncracked  (b  <  b)  in  the  absence  of  an  external 

stress,  the  large  extreme  of  the  particle  size  distribution  is  of  primary 

interest.  An  extreme  value  function  is  thus  selected  for  ascertaining 

1 9 

size  distribution  effects, 


$>  3  1  -  exp[-(bQ/b)k] 


4>(b)db  = 


k(bQ)k 

^T- 


exp[-(bQ/b)k]db 


(32) 


where  bQ  and  k  are  the  scale  and  shape  parameters  respectively.  Com 
bining  eqns.  (31)  and  (32)  gives 


k<b,A  f  b  e-(b°/b)kdb 

[r(l-3/k)]x3/\?  J  bk'3(b  -b) 
c  min  c 

vf*°-3/k)  f3k 

[r(i-3/k)]  J  k  (z1/k  -  l) 
Vf  A(Vbo* 


(33) 


k  k 

where  X  3  (b0/bc)  ,  z  =  (b  /b)  ,  r(k)  is  the  gamma  function  and 

A(bc/bQ,  k)  is  obtained  by  numerical  integration.  The  relative  tougnening 

given  by  eqn.  (33)  is  plotted  in  Fig.  5  for  two  choices  of  ;  (1  x  10"^  and 

.2 

1  x  10  ).  The  peak  value  of  the  toughness  Ty  varies  appreciably  with 
the  width  of  the  particle  size  distribution  (as  reflected  by  the  shape 


parameter  k),  such  that  Ty  increases  as  the  width  decreases  (i.e.  as 
k  increases).  There  are  also  important  effects  of  the  remote  strain  level 
(reflected  in  the  c  dependence)  which  imply  an  influence  of  crack  length 
(and  test  method)  on  the  fracture  toughness.  A  characterization  of  the 
particle  size  distribution  and  a  determination  of  geometry  effects  thus 
become  essential  prerequisites  for  the  interpretation  of  experimental  toughen¬ 
ing  results.  More  significantly,  the  results  imr.ly  that  a  narrow  particle 
size  distribution  is  needed  if  appreciable  toughening  is  to  be  realized, 
in  accord  with  this  mechanism. 

The  toughening  anticipated  by  eqn.  (33)  can  be  related  to  the  intrinsic 

toughness  of  the  matrix,  by  recognizing  that  separation  of  essentially  all 

particles  with  a  size  b>bm.jn  will  have  occurred  as  the  crack  tip  is 

a 

approached,  i.e.,  as  r-*-0  and  hence,  e  -*•«.  The  area  of  matrix  exposed 
to  crack  extension  is  thus  reduced  by  an  amount  related  to  the  total  volume 
fraction  of  particles  V^.  The  toughening  thus  becomes 


where  £  is  the  fraction  of  particles  with  size  >b  . 

min 

£  *  l/<b3>  /  b3ij>(b)db 

bmin 

Typical  values  of  5  are  plotted  in  fig.  6. 


(34b) 
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5.  IMPLICATIONS  AND  CONCLUSIONS 

An  analysis  of  toughening  caused  by  the  separation  of  seconc  phase 
particles  from  the  circumventing  matrix,  within  a  crack  tip  process  zone, 
has  been  presented.  The  most  critical  variables  that  influence  the  tough¬ 
ness  are  the  proximity  of  the  median  particle  size  distribution:  narrow 
distributions  close  to  the  critical  being  preferred.  The  magnitude  of 
the  critical  particle  size  is,  in  turn,  dependent  upon  several  material 
variables:  notably  the  crack  propagation  resistance  of  the  interface,  the 
elastic  moduli,  the  thermal  contraction  coefficient,  and  the  particle  shape. 
These  variables  offer  extensive  opportunities  for  designing  microstructures 
that  optimize  the  toughening  available  from  this  particular  microcrack 
toughening  mechanism. 

The  analysis  also  suggests  that  discrete  microcrack  process  zones  will 
only  be  observed  for  a  limited  range  of  particle  size  and  applied  strain 
levels,  as  characterised  by  a  minimum  size  for  particle  separation  and  a 
maximum  size  for  separation  at  the  level  of  the  applied  strain.  The  tough¬ 
ening  is  thus  anticipated  to  exhibit  trends  that  depend  sensitively  upon 
the  strain  induced  by  microcracking  relative  to  the  applied  strain. 
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FIGURE  CAPTIONS 


Fig.  1.  A  microcrack  process  zone  and  a  schematic  of  the  stress 

modifications  expected  within  the  vicinity  of  the  crack  tip. 

Fig.  2.  A  schematic  indicating  both  the  stresses  that  the  particle  is 
subject  to  in  a  crack  tip  field  prior  to  separation,  and  the 
stress  free  state  of  the  particle  following  separation. 

Fig.  3.  A  predicted  process  zone  contour,  indicating  the  transformation 
zone  width  h-p 

Fig.  4.  Trends  in  toughness  with  volume  fraction  for  given  values  of 
particle  size. 

Fig.  5.  The  predicted  toughening  plotted  for  several  shape  parameters 
k  as  a  function  of  the  ratio  of  the  scale  parameter  to  the 
critical  size  (a)  for  an  applied  strain  characterized  by  t,  -  10 
(b)  ?  *  10'2. 


Fig.  6.  The  variation  of  the  area  fraction  of  particles  separated  as 
a  function  of  the  ratio  of  the  scale  parameter  to  the  critical 
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ABSTRACT 

Brittle  materials  are  subject  to  microcrack  formation  at  grain  boundaries 
and  at  second  phase  particles.  These  cracks  are  induced  by  residual  stress  that 
results  from  Incompatibilities  in  thermal  contraction.  The  development  of 
residual  stress  and  its  partial  relaxation  by  diffusion  (at  elevated 
temperatures)  are  described.  The  evolution  of  microcracks  within  the  residual 
stress  fields  are  then  examined.  Particular  attention  is  devoted  to 
considerations  of  the  critical  microstructural  dimension  at  the  onset  of 
microcracking. 

INTRODUCTION 

Many  properties  of  ceramic  materials  depend  on  the  Incidence  of 
microcracking.  The  most  notable  physical  characteristics  that  exhibit  a  strong 
dependence  on  microcrack  formation  are  certain  mechanical  (fracture  toughness1 
and  fracture  strength2)  and  thermal  (thermal  diffusivity3)  properties.  The 
formation  of  stable  microcracks  is  primarily  related  to  localized  residual 
stresses  that  develop  because  of  thermal  contraction  mismatch  or  anisotropy  (the 
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former  In  multiphase  materials4  and  the  latter  in  single  phase  materials5). 
Significant  progress  has  recently  been  achieved  in  the  analysis  of  microcracking 
events  by  using  a  combination  of  stress  analysis  (based  on  the  Eshelby  concept) 
and  fracture  mechanics.5’6*7  The  intent  of  this  paper  is  to  examine  the 
microcracking  phenomenon  in  order  to  emphasize  both  the  progress  that  has  been 
achieved  and  the  limitations  of  the  available  analyses. 

One  of  the  dominant  characteristics  of  microcracking  is  its  dependence  on 
the  scale  of  the  microstructure.  Typically,  there  is  a  "critical"  microstruc- 
tural  dimension,  ic,  below  which  microcracking  is  not  generally  observed  and 
above  which  a  significant  density  of  microcracks  becomes  evident.4’ 8,9  The 
development  of  a  capability  for  predicting  ic  is  a  primary  objective  of  micro¬ 
cracking  analyses.  A  critical  comparison  with  measured  values  of  iz  is  also  a 
demanding  test  of  the  validity  of  such  analyses. 

The  amplitude  of  residual  stress  fields  produced  by  thermal  contraction 
mismatch  is  independent  of  the  scale  of  the  microstructure.  A  criterion  for 
microfracture  based  on  the  peak  tension  would  not,  therefore,  yield  a  size 
dependence.  This  dilemma  was  first  addressed  by  suggesting8  that  the  onset  of 
microfracture  be  dictated  by  an  equality  of  the  loss  of  strain  energy  and  the 
Increase  in  surface  energy  associated  with  the  microfracture  event.  The  former 
Is  a  volume  dependent  term,  and  the  latter  is  a  surface  area  term;  hence,  a 
critical  size  emerges  in  a  natural  way.  A  reasonable  correspondence  with 
experimental  observation  was  achieved  by  specifying  the  ratio  of  the  final  crack 
size  to  the  dominant  dimension  of  the  microstructure.  A  conceptual  difficulty 
with  the  approach  arises  because  only  the  thermodynamics  of  the  initial  and 
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final  stages  of  the  f.*acture  event  are  considered;  whereas,  fracture  Is  dictated 
by  the  rate  of  energy  change  at  the  critical  condition  for  unstable  crack 
extension.5*9 

Subsequently,  since  size  effects  in  brittle  fracture  often  derive  from 
statistical  considerations, 10  a  potential  role  of  flaw  statistics  was  sug¬ 
gested.11  Notably,  since  fracture  initiates  from  small  inhomogeneities  (pores, 
inclusions,  etc.),  the  spatial  and  size  distribution  of  these  fracture  initi¬ 
ating  sites  can  Influence  the  Incidence  of  fracture.  However,  if  the  size  dis¬ 
tribution  of  these  inhomogeneities  is  independent  of  the  scale  of  the  micro¬ 
structure,  the  fracture  probability  for  a  constant  volume  fraction  of  the 
responsible  microstructural  phase  would  either  be  Independent  of  size  (volume 
distributed  inhomogeneities) 12  or  decrease  with  increase  in  size  (interface 
distributed  inhomogeneities).13  A  statistically  based  argument  must,  therefore. 
Invoke  inhomogeneities  that  Increase  in  size  as  the  microstructure  enlarges. 

This  effect  Is  a  plausible  possibility,  because  inhomogeneities  (such  as  pores) 
tend  to  increase  in  size  during  sintering,  in  direct  proportion  to  the  size  of 
the  grains14  (or  other  microstructural  entities).  However,  in  the  absence  of 
well -substantiated  distribution  functions  to  describe  the  inhomogeneity  size, 
the  quantitative  utility  of  the  statistical  approach  is  limited. 

More  recently,  it  has  been  recognized  that  a  size  effect  can  stem  directly 
from  considerations  related  either  to  the  gradient  of  the  residual  stress 
field5*6  or  to  stress  relaxation.7  For  example,  if  the  fracture  initiating 
inhomogeneity  is  of  sufficient  size  that  it .experiences  an  appreciable  gradient 
of  stress,  then  dimensional  considerations  demand  that  the  fracture  be  size 
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dependent.  Specifically,  the  stress  intensity  factor,  K,  is  given  by5 


1 

K  -  /a*  /  o{x/t)  F(x/a)  d(x/a)  (1) 

0 


where  a  is  the  Inhomogeneity  size,  a  is  the  stress  and  F(x/a)  is  the  appropriate 
Green's  function.  Since  the  stress  can  always  be  expressed  in  the  form 

a(xh)  -  a  n{x/a,  a/t)  (2) 

a 

where  o  is  the  peak  residual  tension,  the  stress  intensity  factor  can  be  written 


a  /-* /  n(x/a,  a/t)  F(x/a)  d{x/a)  =  <{a /i)  (3) 

a  ft  0 


where  <{ a/z)  is  the  function  determined  by  integration.  Now,  if  the  stress 
Intensity  factor  is  equated  to  its  critical  value  for  crack  extension,  Kc, 
Eq.  (3)  yields  a  "critical  size”  given  by 


4 

J2294 A/ jbs 


An  additional  size  influence  derives  from  the  diffusive  stress  relaxation 
that  can  occur  at  elevated  temperatures. 15  The  rate  of  relaxation  will  be  more 
rapid  in  fine-grained  materials,  because  of  the  enhanced  diffusive  fluxes. 
Smaller  residual  stresses  will  thus  obtain  and  the  tendency  for  microcrack 
formation  will  be  reduced. 

The  considerations  of  microcracking  developed  in  this  paper  relate 
primarily  to  the  size  effect  that  derives  from  residual  stress  gradients  and 
relaxation  phenomena.  Beyond  the  scope  of  this  review  are  the  residual 
stresses6*17  produced  by  phase  transformations  during  cooling  (as  in  the  Zr02 
based  alloys)  and  the  effect  of  externally  applied  stress  fields  on  the  residual 
stresses  and  the  onset  of  microcracking. 


RESIDUAL  STRESSES 

Residual  stresses  typically  encountered  in  ceramic  materials  derive  from 
differences  in  thermal  contraction  (anisotropy  of  the  thermal  expansion 
coefficient,  a,  for  a  single  phase  material,  and  contraction  mismatch  for 
multiphase  systems).  Thermal  contraction  differences  are  important  because 
ceramics  are  fabricated  at  elevated  temperatures  (by  hot  pressing  or  sintering) 
and,  during  cooling,  stress  relaxation  (by  diffusion  or  viscous  flow)  becomes 
sufficiently  Inoperative  below  a  temperature  Tg  that  appreciable  local  stresses 
must  develop  from  the  contraction  mismatch.  The  elastic  stresses  that  evolve 
below  Tg  can  be  computed  by  using  adaptations  of  the  Eshelby  approach. 16 

I 

i 

I 
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Several  such  calcul utions  will  be  presented  below.  A  more  difficult  problem  to 
address  is  the  definition  of  Tg-,  this  is  also  examined  in  the  following  section. 

Elastic  Stresses 

The  stresses  that  develop  below  Tg,  calculated  by  using  the  Eshelby 
approach,  are  illustrated  for  the  anisotropic  contraction  of  a  hexagonal  grain 
in  Fig.  1.  This  method  of  calculation  first  extracts  the  microstructural  entity 
(or  entities)  subject  to  shape  deformation  and  its  shape  is  then  allowed  to 
change  (as  characterized  by  the  unconstrained  transformation  strain  e^). 
Subsequently,  its  shape  is  restored  to  the  shape  of  the  matrix  cavity  (by 
exerting  a  uniform  stress)  and  it  is  then  reinserted  into  the  cavity.  Finally, 
interface  forces  are  applied  (of  equal  magnitude,  but  opposite  sign,  to  the 
restoring  stress)  to  achieve  continuity  of  stress.  For  isolated  particles  of 
ellipsoidal  geometry  this  process  yields  a  uniform  stress  within  the  particle; 
hence,  stress  analysis  Is  relatively  straightforward.  More  complex  behavior  is 
expected  for  other  geometries,  such  as  individual  grains  within  a  poly- 
crystalline  aggregate. 

Multiphase  Materials 

The  stresses  within  a  spherical  particle  subject  to  transformation  strains 
eT  (hydrostatic)  and  1  (deviatoric)  are16*17 
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(5a) 


where  a*  is  the  hydrostatic  stress  and  is  the  deviatoric  stress.  The  same 
hydrostatic  stress  level  pertains  for  ellipsoidal  particles,  irrespective  of 
their  shape;  but,  the  deviatoric  stress  is  sensitive  to  the  particle  shape.16 
Two  extremes  are  of  interest.  First,  if  the  particles  and  matrix  have  isotropic 
thermal  contraction  coefficients,  the  resultant  stress  is  exclusively 
hydrostatic;  hence. 


where  -  ap  is  the  thermal  contraction  mismatch  between  matrix  and  particle, 
and  T  is  the  temperature.  For  a  particle  and  matrix  with  similar  elastic 
constants,  Eq.  (5)  becomes 


Second,  if  the  particle  exhibits  anisotropy  of  thermal  expansion  (e.g.,  01  and 
02),  such  that  the  average  expansion  matches  that  of  the  matrix,  then  the  stress 
within  the  particle  is  purely  deviatoric  and  is  given  by 
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which  for  uniform  elastic  properties  reduces  to 


(7  -  5^(0^  -  ag) (T( 

15(1  -  7) 


(9) 


In  general,  therefore,  the  resultant  stress,  for  uniform  elastic  constants 
is 


The  equivalent  result  for  ellipsoidal  geometry  is 


(10) 


(11) 


where  y  is  a  function  of  the  particle  shape.16  For  a  needle,  y  -  1/2;  for  a 
sphere,  7  s  (7  -  5v)/15(l  -  v);  and  for  a  disc,  y  *  ir(c/b)  (2  -  v)/4(l  -  v) , 
where  c  is  the  disc  thickness  and  b  is  the  disc  radius. 


The  stresses  within  the  matrix  are  more  difficult  to  analyze  and,  with  the 
exception  of  spherical18  and  cylindrical  particles, 19-21  have  not  been  computed 
exactly.  Eshelby  shows  that  if  the  harmonic  and  biharmonic  potentials  of  the 
particle  (however  arbitrary  its  shape)  are  known,  then  the  displacements  in  the 
matrix  are  related  to  the  transformation  strain  in  the  particle  by 


T  „  T 

u  *  Jk  Mjk  _  gik  q,jk 

1  ibiruU  -  v)  TTu 
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where  ol^  is  the  stress  derived  by  Hooke's  law  from  the  strain  e^j. 

The  matrix  stresses  can  then  be  obtained  from  the  displacement  derivatives.  In 
the  case  of  a  spherical  particle  subject  to  hydrostatic  strain,  the  matrix 
stresses  are  particularly  simple, 


V  ■  »(r0/r)3,  »99"-(a/2)(r0/r)3  (13) 

where  rQ  is  the  particle  radius  and  r  is  the  distance  from  the  particle  center. 

The  harmonic  and  biharmonic  potentials  for  a  very  long  cylindrical  particle 
(fiber)  have  recently  been  calculated20  and  may  be  used  to  compute  the  matrix 
stresses  from  Eq.  (12): 

4>  =  2irrQ2  in  (r/rQ)  (14a) 


=  irrQ2(r2  -  rQ2)  -  irrQ2  r2  +  rQ2/ 2  in  (r/rQ) 


(14b) 


For  the  case  of  thermal  contraction  mismatch  between  the  fiber  and  matrix, 
described  by  a  hydrostatic  strain  eT,  the  matrix  stresses  are 


c  _  c 
all  "  E 


T  (3  -  2v) 
e  2(1  -  vT 


r02(x2  -  y2) 


T  (3  -  2v) 
e  2(1"”- TF 


(15) 


A  case  of  particular  interest  is  that  of  the  cylindrical  particle  exhibiting  an 
anisotropy  of  thermal  contraction  in  the  plane  perpendicular  to  the  long  axis  of 
the  cylinder  such  that  the  average  contraction  matches  that  of  the  matrix.  The 
stresses  generated  within  the  matrix  are  of  the  form21 
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c  _  E  ell 

ffll  (1  +  V)(l  -  2v)" 


ro'  (xZ  -  yz)  .  V  (x4  -  y4) 
7"^  rz 


e  T 
*  11 , 

4rQ2  y2(3x2  -  y2) 

4(1  -  vZ) 

i* 

where 

•ll  ■  ,ai 

■  V(Tg  -  T> 

2  I A  .  .4  c  2.2, 


(16) 


However,  since  the  stresses  Immediately  adjacent  to  the  interface  are  of 
greatest  interest  for  microfracture  problems,  some  pertinent  information  can  be 
obtained  by  deriving  the  stresses  in  the  matrix  just  outside  the  inclusion.  For 
an  ellipsoidal  inclusion  subject  to  dilatation,  the  stresses  in  the  matrix  can 
be  written  quite  generally  as:16 


cij  =  °^n1nj  ”  l/3<$ij)  (17) 

where  are  the  normals  to  the  ellipsoid  surface.  Particular  values  for  the 
stresses  at  the  particle  matrix  interface  have  been  calculated  for  disc  and 
needle  shaped  particles.22  The  stresses  are  a  maximum  near  the  termination  of 
the  major  axis  of  the  ellipsoid.  However,  as  described  earlier,  the  gradient  in 
stress  (in  addition  to  the  stress  level  at  the  particle  interface)  is  of 
importance  in  dictating  the  size  effect.  As  far  as  the  authors  are  aware,  the 
stress  gradients  around  ellipsoidal  particles  subject  to  a  transformation  strain 
have  not  been  calculated  and  remain  a  subject  for  further  work. 


Single  Phase  Materials 

Grains  in  single  phase  materials  exhibit  relatively  complex  geometric 
configurations,  and  stress  analysis  is  more  complex  than  for  the  isolated 
ellipsoidal  particle.  However,  some  useful  approximations  can  be  obtained  quite 
straightforwardly.  The  general  level  of  residual  stress  within  the  grains  can 
be  obtained  by  simply  requiring  a  grain  to  be  contained  within  an  Isotropic  ma- 
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trix,  with  the  average  properties  of  the  polycrystilline  aggregate,  and 
inserting  the  anisotropic  contraction  coefficients  into  Eq.  (10).  However,  this 
simplification  neglects  the  stress  enhancing  influence  of  grain  junctions,  an 
effect  which  has  important  consequences  for  microfracture. 


Estimates  of  the  stresses  that  develop  in  the  vicinity  of  grain  junctions 
can  be  obtained  by  adopting  two-dimensional  analogues,  such  as  an  array  of 
hexagonal  grains.  The  stresses  that  develop  in  such  an  array  can  be  determined 
by  firstly  establishing  the  resultant  body  forces  at  each  grain  boundary  facet 
(see  Fig.  1).  These  body  forces,  p,  generate  non-uniform  stresses  that 
superimpose  upon  the  uniform  stresses  developed  during  shape  restoration.  Of 
principal  interest  are  the  stresses  at  the  grain  boundaries,  because  these  are 
the  dominant  sites  for  microfracture.  The  stresses  at  a  site  (x,z)  Inclined  at 
an  angle  s  to  the  boundary,  (Fig.  2)  are  of  the  form5 


4itff 


""wxx  _  S  _ (z  +  a  sin  8) _ 

p  cos  8  0  [a2  -  2a(z  sin  8  -  x  cos  g)  +  (x2  +  z2 


)3 


/(I  -  v) - - - °  C-°-S-,-?  ,xj_ - - - -S)  da  (18) 

L  [a  +  2a(z  sin  6  -  x  cos  6)  +  (x^  +  z*)]J 

Equation  (18)  should  be  used  to  obtain  the  stresses  on  boundary  AB  (Fig.  3) 
from  the  body  forces  of  the  four  adjacent  boundaries  (AA' ,  AA",  BB',  BB").  For 
more  remote  boundaries  sufficient  accuracy  can  be  achieved  by  placing  a  single 
force  at  the  grain  facet  center,23  which  represents  the  total  body  force  on  that 
boundary  (Fig.  3).  These  stresses  superimpose  on  a  uniform  stress  (equal  in 
magnitude  to  half  the  initial  stress)  which  derives  from  the  body  forces  on  AB 
coupled  with  the  initial  stress. 


The  component  of  the  stress  from  the  four  adjacent  boundaries  dominates  the 
behavior  in  the  vicinity  of  the  grain  junction.  This  stress  component  is  of  the 
form5 
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I 


pWT*  V3’v)  ln 


. 


2  2  2 
l  +  x  +  z  +  2i(z  sin  B  -  x  cos  6) 


1 - T 

x  +  z 


0 


a  (a  „  y  ^  tan"1  z  sin  B  -  X  COS  b! 

A2(8,v,x,z).  tan  j^li  cos  S  +  x  sin  3)  J 

■1  T  z  tan  8  -  x  ")  ,  ,  ,  „  „  „  _x 

J_2'(z  +  T TaTTsTj  A3(l»  8’  v’  x‘  z) 


-  tan 


(19) 


where  Ai,  A2,  and  A3  are  relatively  complex  functions  in  the  range  ±  2ir.  The 
logarithmic  term  is  singular  at  the  grain  junction  and  is  thus  the  most 
influencial  with  regard  to  microcrack  formation.5 


The  specific  stress  magnitudes  that  develop  depend  on  the  relative 
orientations  of  the  grains  circumventing  the  boundary  of  interest.  Preliminary 
calculations  have  been  conducted  for  the  orientation  that  has  been  assumed  to 
yield  the  maximum  stress;  this  depicted  in  Fig.  4a,  and  the  results  are  shown  in 
Fig.  4b.  Calculations  for  more  general  grain  orientation  relations  are  now  in 
progress.23  These  results  will  provide  a  full  perspective  of  residual  stress 
distributions  in  polycrystal  line  aggregates  in  which  there  is  a  random 
distribution  of  contraction  anisotropy  orientations. 


STRESS  RELAXATION  EFFECTS 

Stress  relaxation  in  ceramics  occurs  primarily  by  diffusion  (or  by  viscous 
flow  in  the  presence  of  an  amorphous  phase).  These  relaxation  processes  are 
usually  motivated  by  local  gradients  in  hydrostatic  stress,  and  thus  occur  in 
response  to  localized  thermal  contraction  stresses,  while  the  material  is  at 
elevated  temperatures. 


Multiphase  Materials 

For  isotropic  multiphase  materials,  there  is  no  gradient  of  hydrostatic 
stress  within  the  isolated  phase  (Eq.  5).  However,  large  shear  stresses  exist 
within  the  surrounding  matrix  (note  that  the  hydrostatic  stress  is  zero).  The 
shear  stresses  within  the  matrix  cause  grain  boundary  sliding,  and  diffusive 
deformation  will  occur  in  response  to  local  normal  stresses  induced  by 


12 

J2294A/ j  bs 


sliding.24  The  deformation  /i el d  will  be  similar  to  that  of  a  cavity  subject  to 
internal  pressure.  Initially,  racial  flow  in  the  matrix  will  redistribute  the 
residual  stress.  Reduction  of  the  stresses  will  begin  when  the  redistribution 
extends  across  the  sample,  following  the  onset  of  interaction  between  the  stress 
fields  around  adjacent  particles.  The  authors  are  not  aware  of  solutions  for 
this  problem,  although  the  analysis  Is  relatively  straightforward. 

Polycrystalline  Single  Phase  Aggregates 

It  has  already  been  demonstrated  that  anisotropic  thermal  contraction  in 
polycrystalline  aggregates  develops  tensile  or  compressive  stresses  on  grain 
boundaries.  A  gradient  of  chemical  potential  suitable  for  diffusive  relaxation 
(Fig.  5)  is  thus  established.  The  "initial"  stress  involves  singularities  near 
grain  junctions  (Fig.  4).  But  the  singularities  are  weak  (logarithmic)  and 
should  be  rapidly  dispersed  by  localized  diffusive  fluxes.  The  rate  controlling 
relaxation  process  involves  diffusion  between  adjacent  grain  facets  (Fig.  5), 
such  that  atoms  are  removed  from  the  boundaries  subject  to  compression  and  are 
deposited  on  boundaries  under  tension.  If  it  is  as*sumed  that  the  relaxation 
times  are  sufficiently  rapid  that  atom  deposition  and  removal  occurs  uniformly 
along  the  respective  grain  boundaries,  a  parabolic  "steady- state"  stress  dis¬ 
tribution  must  develop  along  the  boundaries  during  the  relaxation  process.  The 
extent  of  strain  relaxation  can  then  be  deduced  by  using  well-established  mathe¬ 
matical  procedures  for  diffusive  flow.  This  mode  of  analysis  is  only  permiss- 
able  when  the  diffusivities  are  large,  notably  at  the  highest  temperatures.  The 
stress  evolution  at  intermediate  temperatures  requires  "transient"  solutions 
involving  more  complex  formulations.  Such  analyses  have  not  yet  been  performed. 
Currently,  therefore,  it  Is  only  possible  to  obtain  approximate  solutions  by 
permitting  "steady-state"  relaxation  above  a  "freezing"  temperature  T$  and 
Invoking  fully  elastic  stress  development  below  Tg. 

The  stress  relaxation  problem  can  be  posed  by  first  establishing  the  normal 
elastic  displacement  5i  of  the  boundaries  during  cooling  (the  driving  force  for 
the  diffusive  flow.  Fig.  5)  and  the  displacement  relaxation  52  due  to  diffusion. 
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Then,  the  resultant  displacement  5  (=  Si  -  <$2),  which  determines  the  level  of 
the  relaxed  stress,  can  be  derived.  The  solution  for  a  constant  cooling 
rate  T  will  be  presented.25 

The  elastic  stress  level  on  a  grain  boundary  normal  to  the  direction  of 
maximum  contraction  is 

8  E(T  -  Tt)Aa 
ffxx  =  (1  +  v) 

where  TQ  is  the  initial  temperature,  Aa  is  the  deviation  of  the  contraction 
coefficient  from  the  average,  and  8  is  a  coefficient  that  depends  on  the 
orientations  of  the  adjacent  grains.  The  corresponding  elastic  displacement  is 
• 

f  _  /3  l  8  T  t  Aa  » 

51  -TTTTUJ -  <21> 


The  relaxations  of  these  displacements  by  diffusion  are  governed  by  the 
relation26 


d2q(x,t) 


kld2(t) 


where  Db5b  is  the  diffusion  parameter,  n  the  atomic  volume  and  5 

is  assumed  to  be  uniform  (as  noted  above).  Integration  of  Eq.  (22)  gives  the 

stress  distribution 

o(x,t)  =  -$£-  +  CjX  +  C2  (23) 

where  E  *  kT52(t)/5l  0^5^  and  C*  and  C2  are  the  integration  constants.  The 
positions  of  zero  flux  (do/dx  *  0)  In  the  system  are  at  the  grain  facet  center  0 
and  at  the  grain  junction  J  (Fig.  5).  Hence,  since  the  flux  must  be  continuous 
at  the  grain  corner  A,  the  constant  Ci  must  be  equal  to  £*/ 2.  If  a0  Is  the 
stress  at  A,  the  stress  distribution  becomes 

o(x,t)  *  -  1*1  +  +  <jo  (24) 
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The  equivalent  average  stress  is 

<0>(  t)  s  "yp  +  (25) 

Volume  conservation  requires  that  the  volume  of  material  deposited  on  the 
tensile  boundary  must  equal  the  volume  removed  from  the  compressed  boundary.  The 
stress  at  the  grain  corner  thus  becomes 

aQ  *  Cl2/ 12  (26) 

and  the  average  stress  reduces  to 
kTz2  5 ? 

‘mjr  (27> 

D  D 

The  average  stress  on  each  boundary  must  also  be  related  to  the  resultant 
displacement  of  the  grain:  . 

/5Z  d<0>  -  &  * 

TTST  ’  h  '  52  (28) 

*  • 

Substituting  5*  from  Eq.  (21)  and  S2  from  Eq.  (27),  and  noting  that  is 
temperature  dependent 


°b  ■  °o  e 


-Q/kT 


where  Q  is  the  activation  energy  for  boundary  diffusion,  the  following 
differential  equation  obtains 

_  12n«J)_E  -Q/kT 


The  solution  to  this  equation  must  be  conducted  numerically.  However,  an 
approximate  series  expansion  may  be  derived25  for  comparison  with  the  elastic 
result  expressed  in  terms  of  a  "freezing"  temperature  Tfl: 
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<0>{1  +  v) 


From  this,  Tg  is  given  approximately  by 

Tq - — - rr  (31] 

a  in[12fl  O05b£//3nlcr  T] 

The  trends  in  Tg  with  the  influential  variables  are  immediately  apparent. 
In  particular,  Tg  increases  as  the  grain  size  and  cooling  rate  increase  and  as 
the  diffusivity  and  modulus  decrease.  This  behavior  is  exemplified  for  A1 2O 3 
(D05b  *  10’9  m3  s’1,  Q  *  100  kcal/mole,  Q  -  10’29  m3,  n  =  30,  £  =  420  GN  m2), 
for  which: 

T  =  50,000 

9  29.5  -  in  i3  T 

where  1  is  in  microns,  T  in  K  s’1  and  Tg  is  in  K.  Specifically,  for  1  =  1  un 
and  T  »  IK  s"1,  Tg  *  1695K;  whereas  for  i  *'10  ym  and  T  =  IK  s’1,  Tg  =  2210K. 


MICROCRACK  FORMATION 


The  formation  of  microcracks  at  grain  boundaries  or  at  second  phase 
particles  has  been  considered  to  depend  on  the  existence  of  a  distribution  of 
small  inhomogeneities  that  pre-exist  at  the  boundaries  (especially  at  three 
grain  junctions)  or  interfaces. 5,6,7  These  inhomogeneities  have  been  proposed 
because  the  stress  intensification  levels  associated  with  the  residual  stress 
fields  do  not  appear  to  be  of  sufficient  magnitude  to  induce  fracture  in  defect 
free  material  (although  further  study  is  needed  to  establish  whether  this  possi¬ 
bility  can  be  discounted).  The  role  of  the  proposed  Inhomogeneities  is  to  fur¬ 
ther  enhance  the  stress  intensification  to  a  level  suitable  for  microcrack  evo¬ 
lution.  It  is  certainly  reasonable  to  presume  that  inhomogeneities  do  exist  at 
boundaries  or  Interfaces  in  ceramics,  e.g.,  small  pores  remaining  at  grain 
triple  points.  However,  little  effort  has  yet  been  devoted  to  the  elucidation 
of  the  inhomogeneities  that  Induce  microfracture  in  specific  miocrostructural 
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situations.  It  Is  thus  generally  assumed  that  the  inhomogeneities  exhibit  the 
stress  concentrating  properties  of  small  microcracks:  a  presumption  that  is 
evidently  an  over-simplification.  Thereafter,  stress  intensity  factors  can  be 
calculated  (from  the  residual  stress  levels)  and  compared  with  the  critical 
values  for  grain  boundary  separation.  Approximate  stress  intensity  factors  are 
conveniently  calculated  with  a  superposition  method,  based  on  the  prior  stress 
field. s»6  A  typical  example,  illustrated  in  Fig.  6,  is  where  a  microcrack 
develops  along  two  symmetrically  stressed  grain  facets,  initiating  at  the  common 
triple  junction.  The  stress  intensity  factor  for  such  a  crack  is  given  by 

-(-1  *  <(a /i)  (32) 

EAaAT/l 

where  «(a/jt)  is  the  function  plotted  in  Fig.  6.  It  is  noted  that  the  stress 

A 

intensity  factor  exhibits  a  maximum,  K.  This  is  typical  of  crack  extension  in 
residual  or  spatially  varying  stress  fields.  The  principal  maximum  in  the 
present  analysis  essentially  coincides  with  a  crack  front  located  at  the  first 
trip>e  junction,  where  the  residual  stress  changes  sign  (i.e.,  the  residual 
stress  becomes  compressive  along  the  impinging  boundaries).  Equating  K 
to  the  boundary  separation  resistance,!^,  yields  an  absolute  minimum  condition 
for  the  formation  of  microcracks.  This  corresponds  to  an  upper  bound  for  the 
critical  grain  facet  size; 


i 


c 


(1  +  v) 
<  EAaAT 


2 


(33) 


A 

where  *  is  the  magnitude  of  the  normalized  stress  intensity  factor  at  the 
maximum.  Estimates  of  specific  values  of  the  grain  facet  sizes  that  induce 
microcracking  involve  statistical  considerations  based  on  flaw  distributions. 

More  exact  calculations  of  the  stress  intensity  factor  can  be  obtained  by 
using  numerical  (finite  element  or  finite  difference)  methods.  A  recent 
example7  is  the  use  of  a  finite  difference  scheme  for  calculating  the  stress 
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intensity  factor  for  a  crack  at  tne  interface  of  a  spherical  second  phase 
particle.  A  convenient  use  of  the  finite  difference  method  involves  the 
calculation  of  the  strain  energy,  U,  as  a  function  of  crack  length,  a.  The 
stress  intensity  factor  is  then  deduced  from  the  crack  length  derivative  of  the 
strain  energy.  A  maximum  value,  K,  is  obtained.  The  corresponding  absolute 
minimum  requirement  for  microcrack  initiation  is 


1.2(K*/<c>)2 


(34) 


where  K*  is  the  resistance  of  the  interface  to  fracture  and  <a>  is  the  stress  in 
the  uncracked  particle. 


Further  progress  in  the  elucidation  of  microfracture  is  achieved  by 
incorporating  defects  that  reduce  the  critical  size  below  the  upper  bound 
values.  Little  progress  has  yet  been  achieved  in  selecting  appropriate  defects 
and  defect  size  distributions;  although  results  could  clearly  be  obtained  by 
selecting  arbitrary  distributions  (such  as  one  of  the  extreme  value  distribu¬ 
tions).  Careful  obsevations  coupled  with  pertinent  stress  intensity  factor 
calculations  are  needed  to  establish  a  more  fundamental  appreciation  of  micro¬ 
cracking.  Comparison  of  the  available  calculations  of  K  with  experimental 
observations  of  microcracking  suggest5*6'7  that  triple  point  defects  in  the  size 
range,  0.1  <  2a/z  <  0.3,  are  typically  involved  in  the  microcrack  initiation 
process.  However,  more  direct  observations  and  further  calculations  are  needed 
to  substantiate  and  refine  this  result. 


CONCLUSION 

This  paper  has  described  methods  for  calculating  the  residual  stresses  that 
develop  at  the  microstructural  level  because  of  thermal  contraction 
inhomogeneity.  The  stresses  are  typified  by  locally  large  amplitudes  (with  some 
singularities)  and  rapid  gradients.  These  characteristics  are  central  to  the 
onset  of  microcracking,  both  microcrack  Initiation  and  arrest. 

The  stress  level  is  also  shown  to  depend  on  the  rate  of  stress  relaxation 
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at  elevated  temperatures,  by  diffusion  or  viscous  flow.  The  relaxation  rate  is 
a  strong  function  of  microstructure:  rapid  relaxation  rates  obtain  in  fine 
grained  materials  or  in  materials  containing  an  amorphous  boundary  phase.  These 
relaxation  effects  have  been  demonstrated  to  be  manifest  in  an  effective  freez¬ 
ing  temperature,  a  temperature  at  which  elastic  residual  stresses  begin  to 
develop. 

The  onset  of  microcracking  within  the  residual  stress  field  has  been 
considered  to  depend  on  the  presence  of  small  microstructural  inhomogeneities 
(such  as  voids)  at  the  susceptible  interfaces.  These  features  certainly  exist, 
but  have  not  yet  been  uniquely  correlated  with  the  onset  of  microcracking.  By 
treating  these  pre-existent  inhomogeneities  as  crack-like  entities,  stress 
intensity  factors  have  been  calculated.  The  level  and  variation  in  stress 
intensity  factors  indicate  the  potential  for  microcrack  initiation  and  arrest  at 
grain  boundaries  and  interfaces. 

In  particular,  a  lower  bound  for  the  critical  microstructural  size  needed 
to  initiate  microcracks*has  been  identified  (no  microcracks  can  be  observed  at 
size  levels  below  this  bound).  The  actual  formation  of  microcracks  above  the 
lower  bound  depends  on  the  statistical  characteristics  (size  and  spatial)  of  the 
pre-existent  inhomogeneities.  This  issue  has  not  been  addressed. 
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FIGURE  CAPTIONS 


Fig.  1 

Fig.  2 

Fig.  3 

Fig.  4 

Fig.  5 

Fig.  6 


A  schematic  indicating  the  Eshelby  method  for  calculating  the  residual 
stresses  and  strains  generated  by  anisotropic  thermal  expansion  of  an 
hexagonal  grain. 


The  linear  boundary  segment  used  to  compute  the  relaxation  stresses, 
showing  the  coordinate  system  (x,  z). 


An  hexagonal  grain  array  showing  the  body  forces  used  to  calculate  the 
stress  at  the  central  facet  AB. 


The  grain  configuration  that  yields  large  values  of  the  residual  stress 
at  facet  AA'  and  the  stresses  calculated  to  exist  along  that  facet. 


The  elastic  and  diffusion  displacement  that  occurs  during  cooling, 
indicating  the  direction  of  the  diffusive  flux.  Also  shown  are  the 
stresses  that  develop  during  steady-state  diffusive  flow. 


The  variation  of  the  normalized  stress  intensity  factor  with  crack 
length  for  the  grain,  crack  configuration  indicated  on  the  inset. 
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ABSTRACT 

An-  indentation  method  for  determining  the  adhesion  of  inter¬ 
faces  between  thin  films  and  substrates  has  been  developed.  The 
method  provides  a  quantitative  measure  of  the  interface  fracture 
resistance  and  has  the  advantage  of  simplicity  and  reproducibility. 
The  method  has  been  demonstrated  for  a  range  of  ZnO/Si  systems  and 
the  adherence  has  been  correlated  with  acoustic  properties. 

I 

INTRODUCTION 


The  quantification  of  adhesion  between  films  (or  coatings) 
and  substrates  is  an  important  but  inadequately  resolved  problem. 

A  meaningful  method  for  measuring  adhesion  is  needed  both  for  qual¬ 
ity  control  and  for  the  development  of  improved  film  properties 
,(e.g.  monitoring  variations  of  adhesion  with  different  treatments 
and  fabrication  methods) .  This  need  has  become  increasingly  im¬ 
portant  with  the  recent  expansion  in  the  use  of  thin  film  substrate 
in  the  electronics  industry. 

Ideally,  the  requirements  of,  an  adhesion  test  are:  quantita¬ 
tive  interpretation,  applicability  to  a  wide  range  of  film  thick¬ 
ness  and  film/substrate  combinations,  economy  (minimal  material  use 
and  .machining  requirements),  reliability,  sensitivity,  and  adapta¬ 
bility  to  routine  testing.  Despite  many  attempts  to  develop  tech¬ 
niques  for  adhesion  measurement^-”®,  an  entirely  adequate  quantita¬ 
tive  test  has  not  yet  emerged.  The  objective  of  the  present  study 
is  to  develop  an  adhesion  test,  based  on  indentation  fracture  - 
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procedures,  tha'  satisfy  these  criteria.  _ _ 

Adhesion  tests  based  on  fracture  mechanics  principles  are  the 
most  likely  source  o;.  quantitative  procedures.  Fracture  mechanics 
methods  exhibit  several  attractions.  Firstly,  an  ahesion  parameter,^ 
such  as  the  interface  fracture  resistance  (or  toughness)  K^nter  __ 

can  be  defined.  Secondly,  the  potential  exists  for  developing  _ - 

sound,  theoretically,  based, .  measurements  of  interface,  toughness.  _ 
Finally,  fracture  mechanics  provides  a  mathematical-framework  for_  _V~ 
using  the  fracture  toughness  as  -a  design  parameter..  —  Rowever^_cur= 

rently -available- tests^j ^-0  are  of  restricted  application  and_are _ L. 

too  complex  for  most  purposes. 

The  test  method  proposed  in  this  paper— involves-the  -intro- — - 
duction  of  a  mechanically,  stable -crack  into -the— interface,— bv—em-— 
ploying  conventional  -indentation  procedures . — The  -resistance  -to  — 

propogation -of— ther^crack- along  the— interface— is— then-used— as— a - 

measure  of  -adheslon;^-^  Analysis- -of — the— r-esult-s—ia— based— on— recently 
developed  indentation  fracture  mechanics11*13.  The  test  is  dem¬ 
onstrated  to  yield— drsensitive  meas u r e,  of - in terf a c e— to ughne s  s. - 

Also,-by~virtue— of— Its -inherent-  simplic-ityy-and-STiiall— scale-,— it - 

appears  to  meet  most  of  the  requirement^  for  a  widely  applicable 
adhesion  measurement  technique. 

INDENTATION  FRACTURE 

The  "well  developed"  deformation/fracture  pattern  which  re¬ 
sults  from  indentation  of  the  surface  of  a  homogeneous  brittle 
solid  with  a  Vickers  pyramid  is  shown  schematically  in  figure  1. 

The  elastic/plastic  nature  of  the  deformation  response  is  charac¬ 
terised  by  the  plastically  deformed  zone  (with  a  residual  hardness 
I  impression)  surrounded  by  cracks.  Two  distinct  crack  systems  are  -t-r 
apparent;  the  "median/ radial"  system  comprising  two  orthogonal,  4-r 
semi-circular  cracks  parallel  to  both  the  load  axis  and  an  indenta — ra¬ 
tion,  diagonal,  and  the  penny-shaped  "lateral"  crack  parallel  to 
the  surface  and  centered  near  the  base  of  the  plastic  deformation 
zone.  Both  crack  systems  develop  mainly  during  the  unloadings  half  -=-?• 
cycle  and  are  driven  primarily  by  the  residual  elastic/plastic  ■di¬ 
stress  field..  1 

I  ^  ~ 

I  - 

Recent  analysis  of  these  crack  systems  has  provided  under-  -^=- 
standing  of  the  behavior  of  flaws  in  brittle  surfaces,  and  enabled  -=4 
development  of  predictive  models  if or  material  degradation  processes^— r 
The  radial/median  crack  system  is  pertinent  to  strength  degrada- 
tion  (static  contact13*14,  particle  impact15*16,  machining  dam-  — 
age17)  and  the  lateral  crack  system  to  material  removal  processes  —  - 
(erosion  and  wear  »  ).  In  addition,  the  derivation  of  the  func¬ 

tional  dependence  of  the  radial  crack  dimension  or  the  pertinent 
parameters  (contact  load  P  ,  material  hardness  H  .--elastic - 


INDENTATION 


Fig.  1.  A  schematic  showing  crack  patterns  generated  by  Vickers 
indentation  in  a  homogenous  material,  a)  Four  radial 
cracks  propagating  along  the  indentation  diagonal  direc¬ 
tion.  b)  Sectional  view  of  the  indentation;  the  dotted 
region  represents  the  plastic  zone  created  by  indenta¬ 
tion.  The  median  crack  coalesces  with  the  radial  crack 
to  form  a  half-penny  shaped  crack.  Lateral  cracks  propa¬ 
gate  parallel  to  the  surface. 

modulus  E  ,  and  toughness  Kc  )  has  led  to  a  simple  and  efficient 
technique  of  toughness  measurement^* 21.  The  method  has  been,  ex¬ 
tended  to  measure  the  interface  toughness  of  very  thick  coatings 
by  sectioning  normal  to  the  interface  and  aligning  the  indentation 
such  that  one  "radial  crack  propagates  along  the  interface^.  How¬ 
ever  the  method  cannot  be  applied  to  thin  coatings  (less  than  a>1 
mm)  because  of  free  surface  effects  on  the  crack  propagation. 

A  more  suitable  basis  for  the  present  film  adhesion  measure¬ 
ments  is  provided  by  the  lateral  crack  system.  The  film  surface 
parallel  to  the  interface  is  indented,  and  the  lateral  crack  is 
induced  to  propagate  along  the  interface  (figure  2).  Then  measure¬ 
ments  of  the  crack  radius  c  ,  the  indentation  load  P  ,  the  inden- 

tational  angle  2a  ,  and  the  film  thickness  h  .,  .will  he  demon= 

strated  to  provide  a  measure  of  the  toughness  K|nt.e-r  .  -.Lateral- 
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Fig.  2.  A  schematic  of  an "Indentation  made "on"  a 'film/substrate 
s ampTeT  The- later al  "c  r  a ck  pro pagates  prefe fen tially 
along  the  interface. 


cracks  in  a  homogeneous  material  usually,  but  not  always,  form 
~  near  the  base  of  the  plastic  zone.  However,  the  driving  force  for 
“  lateral  crack  propagation  exists  over  a  range  of  depths  and  the 
~  actual  depth  at  which  a  crack  propagates  is  influenced  by  the 
37  availability  of  crack  nucleation  sites ^3.  The  proposed  test  is 
~7  based  on  the  rationale  that  an  interface  (in  the  vicinity  of  the 
plastic  zone)  with  a  lower  toughness  than  that  of  either  the  film 
7" r  material  or  the  substrate  material  will  be  the  preferred  site  for 
77  lateral  crack  formation.  When  fracture  cannot  be  induced  at  the 
~ 7  interface  (but  occurs  in  the  film  or  substrate)  it  can  be  immedi- 
77  ately  concluded  that  the  interface  toughness  is  at  least  as  large 
T-  as  that  of  the  weaker  material  constituent. 

7“  The  present  study  is  confined  to  combinations  of  ZnO  on  Si 

~~r  and  Si  on  Si  because  the  materials  have  similar  values  of  both 
7'  hardness  and  elastic  modulus  (the  parameters  that  dictate  the  mag- 
77  nitude  of  the  residual  stress  which  provides  the  crack  driving 
force).  This  choice  simplifies  the  theoretical  analysis  and  per- 
7”  mits  direct  correlation  between  theory  and  experiment.  The  effects  77 
77  of  a  mismatch  in  hardness  and/or  elastic  modulus  will  be  addressed  777 
77  in  a  subsequent  analysis.  777 

“  EXPERIMENTAL  OBSERVATIONS  7r 

The  viability  of  the  indentation  technique  for  measuring  ad¬ 
hesion  was  firstly  established  by  examining  the  influence  of  a  sub¬ 
surface  interface  on  the  development  of  lateral  fractures.  This 
was  achieved  by  comparing  lateral  cracks  obtained  in  a  Si  single 


crystal  with  those  that  develop  at  the  in'  irface  between  a  2  ^m 
polycrystalline  Si  film  deposited  onto  si  igle  crystal  Si  (fig.  3). 

A  circular  sub-surface  lateral  crack,  located  at  the  interface, 
was  clearly  evident  in  the  latter,  while  a  sea.  ler,  less  distinct 
lateral  crack  developed  in  the  single  crystal.  These  observations 
indicate  that  interfaces  are  preferred  sites  for  the  lateral  frac¬ 
ture  process,  thereby  satisfying  the  basic  prerequisite  for  apply¬ 
ing  the  indentation  technique  to  the  adhesion  problem..  All. sub¬ 
sequent  experiments  were  conducted  on  ZnO/Si  systems. - 

19  18  9  7'  - - - 

It  can  be  anticipated  from  previous  work  ’  V  .in  indents-., 

tion  fracture  that  the  radius  of  the  lateral  crack  will  depend  on 

the  indentation  load  and  the  film  thickness- (i-e.  the. . crack— loca- _ 

tion)  in  addition- to  the  interface  toughness -(the  parameter-of - 

interest).  A  typical  inf luence  of  the  load  on  the  -crack  length. - 

is  illustrated  for  ZnO/Si  in  fig.  4. —  Quantitative  application  of 
the  indentation  method  will  require  determination  of-the  -functional 
dependence  ofuihe^crack-length  on  -each-of— these-variables, — This — 1 
issue  will  be  addressed  in  section  4.  At  this  jucture,  some  im¬ 
portant-trends— in  adhesion  are ; established ._-_This  is-achieved-by — 
comparing  -crack-lengths-for-the -same -load  and thickness - - - 


Fig.  3.  Optical  micrographs  of  indentation  on  Si.  A)  An  Indenta¬ 
tion  made  on  single  crystal  (111)  Si  surface,  indicating 
four  radial  cracks,  and  S)  an0 Indentation  made  on  thin 
film  specimen,  2  pro  poly  S1/500A  Si02/ (111) Si._  A  circular 
.  sub-surface  lateral  crack,  located  at  the  interface,  can 
be  seen  using  an  interference  contrast  method. 


Fig.  4.  Optical  micrographs  illustrating  the  effect  of  the  inden-  -f 
tation  load  on  the  lateral  crack  length  for  a  ZnO/Si 
system.  A)  5  N  and  B)  10  N. 

obtained  on  indented  samples  of  ZnO/Si,  prepared  under  different 
conditions.  In  some  instances,  the  adhesion  trends  are  correlated  -» 
with  other  property  measurements,  in  order  to  demonstrate  the 
general  utility  of  indentation  as  a  reproducible,  comparitive 
^aethod  for  determining  adhesion. 

The  indentation  results  obtained  at  300  g.  for  the  following 
two  systems, g6  ym  ZnO/lOOOA  SiO»/ (111) Si-N  type  and  6  ym  ZnO/ 

1500A  Pt/500A  Ti/(111)  Si-N  type,  are  illustrated  in  fig.  5.  It 
is  noted  that  faceted  lateral  cracks  develop  in  latter,  indicating 
the  existence  of  preferred  paths  of  fracture  resistance  (typical 
of  lattice  orientation  effects  in  single  crystals),  and  that  the 
average  crack  radius  for  this  system  is  smaller.  Superior  adhe¬ 
sion  thus  obtains  for  the  system  with  the  Pt/Ti  interlayer^. 

"^It  should  not  be  concluded  that  the  Pt/Ti  interlayer  invariably 
provides  superior  adhesion.  Other  film  formation  parameters--are— 
of  substantial  importance  with  regard  to  adhesion.-  -  -  - 


Fig.  5.  Optical  micrographs  of  300  g  Indentations  made  on  A) o 
6  ym  ZnO^lOOOA  Si02/(lll)Si-N  type,  B)  6  ym  Zn0/1500A 
Pt  -  500A  Ti/(lll)5i-N  type.  The  average  crack  radius 
is  smaller  in  the  latter,  indicating  a  superior  adhesion. 

24 

Independent  measurements  of  the  acoustic  properties  indicate 
appreciably  larger  acoustic  signals  for  the  latter.  A  correlation 
between  adhesion  and  acoustic  properties  is  thus  established.  The 
origin  of  the  correlation  is  not  evident  at  this  juncture,  but  it 
is  presumed  that  deposition  conditions  that  develop  a  good  quality 
interface  also  lead  to  films  of  high  quality  which,  in  turn,  pro¬ 
vide  good  acoustic  properties. 

0 

Another  series  of  tests  on  10  ym  ZnO/lOOOA  Si02/ (111)  Si  sam¬ 
ples  are  summarised  in  fig.  6.  One  sample  exhibits  superior  adhe¬ 
sion  to  the  other  two,  as  manifest  in  a  faceted  lateral  crack  con¬ 
figuration,  with  a  smaller  average  crack  radius.  Sectional  views 
of  these  samples  (fig.  7)  indicates  that  the  sample  with  the 
superior  adhesion  consists  of  a  ZnO  film  with  columnar  grains, 
oriented  with  (0001)  normal  to  the  substrate  surface.  This  film 
also  exhibits  greater  optical  transparency.  Additionally,  x-ray 
rocking  curve  experiments  show  that  the  full  width  at  half  maximum 
of  the  (0001)  diffraction  peak  is  smaller  for  the  sample  with  su¬ 
perior  adhesion  than  for  the  other  samples,  (1.9°  compared  to-2-.8r) 
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ther  studies  have  indicated  that  films  with"  a  highly  oriented 
structure,  good  optical  quality,  and  columnar  grain  structures 
typically  provide  optimum  acoustic  properties'-^ .  Again,  there¬ 
fore,  the  measure  of  adhesion  obtained  by  indentation  seems  to  re¬ 
flect  trends  in  other  important  physical  properties. 

The  sensitivity  of  the  indentation  test  to  inferior  adhesion 
is  demonstrated  in  fig.  8.  Very  large  lateral  cracks  develop  on 
poor  quality  films;  in  some  instances  causing  sections  of  the  film 
to  be  removed  (fig.  8B) .  The  variability  in  the  extent  of  .lateral 
fracture  obtained  at  various  locations  over  the  surface  indicates 
that  adhesion  in  such  films  exhibits  substantial  spatial  variation 
(good  quality  films  show  very  littlescatter- in .  the— Lateral.  crack _ 
dimension). —  —  -I;-- - - - - - - - 

Finally,.  the-application-of.  indentation— tests, -as-a  simple _ 

means  for  monitoring  trends  in  adhesion  is- emphasized. —  For— this _ 

purpose ,  10  ym  ZnO/Si-samples-have— been -annealed— in  -air— and— in - 

vacuum  and  then  the  lateral  fracture  extension  determined  in  each 


case  (fig.  9).  It  is  immediately  evident  that  vacuum  annealing  _ 

substantially  enhances  the  adhesio  ,  while  annealing  in  air  causes  _ 

a  corresponding  degradation. 

THE  INTERFACE  TOUGHNESS  _ 

A  quantitative  measure  of  the  interface  toughness  can  be  pro-  — 
vided.  by„  developing. an  adequate  model  of  lateral-  crack  extension.  --  — 
The  extension  oL  well  developed  indentation  fractures  has  recently-  -- 
been  subject  to  successful  analysis  by  treating  the  indentation  ...  ... 
plastic  zone  as  a  precompressed  spring  that  provides  the  driving 
force  for  crack  extension,  and  relaxes  as  the  crack  extendsll>l2.  i. 

This  same  approach  can  be  used  to  anticipate  the  crack  extension _ 

c  along  an„interface,. -.located  at_depth,  h..^_beneath..the  surface  1_ 
(fig.  10).  .  The  analysis,  in.  this  case  _dif f ers_from..that_in„a_homor.  _1 
geneous  material^2  only -in  the  . treatment  of  depth  h-(for  inter- 
face  cracks  h. ...  is._constant.and  equal_to  the  .film, thickness,,  where- 

as  cracks  in^a— homo geneous-material-are— located . at . the-base.  .of.,  the  2. 

deformation  zone  and  h  varies  with  indentation  load) .  The  result  2. 
obtained— for -the— in t erf ace  -latexa  1— c racks— is  j  -  2. 


«ti  -  p„/gi1/2  ;1/4 


a) 


where  p  is  the  peak  indentation  load,  PQ  is  a  threshold  load 
for  crack  formation  and  a  is  given  by; 


a2  -  O],  h3/2  H1/2/K 


IC 


(2) 


where  is  a  material  independent  coefficient  that  can  be  deter-  3 

mined  by  calibration,  on  a  homogeneous  material  with  known  H  and  j 
Kj£  ;  Such  calibration  experiments  have  been  performed  on  glass^2; 
experiments  which  yield  a  value  for  ai  =  6.4  x  10-8.  Hence,  an 
experimental  determination  of  a  ,  obtained  by  fitting  crack  lengths2 
at  several  load  levels  to  eqn.  (1),  permits  Kj£ter  to  be  obtained  J 
from  eqn.  (2);  provided  the  f ilm^thickness  h  is  also  measured 
(H  is  determined  directly  from  p  and  the  indentation  diagonal) . 

Experimental  results  obtained  for  a  ZnO/SiO^/Si  system  are 
shown  in  fig.  11.  A  fit  of  data  to  eqn.  (1)  indicates  that  the 
threshold  load  is  P0  =  3N  and  that  a  =  4.5  x  10'^.uN"1'  .  In¬ 
serting  h  ■  10  um  and  H  =  8.5GPa  into  eqn.  (2),  the  interface 
toughness  is  determined  to  be  0.1  MPaVm.  This  compares  with  ma¬ 
terial  toughness  of,  0.6  MPaVm  for  Si,  ~0.7  MPaVm  for  Si02  and 
•vlMParVm  for  ZnO.  The  interface  thus  exhibits  a  lower  fracture 
toughness  than  the  constituent  materials,  as  required  for  the  ob¬ 
servation  of  interface  fracture.  Preliminary  credence  in  the  model 
is  thus  provided.  Much  additional  experimentation  is  needed,  how¬ 
ever  to  determine  if  the  analysis  correctly  predicts  trends  in 
film  thickness,  load  and  interface  toughness.  . 


Tig.  9.  Optical  micrographs  of  1000  g  indentations  on  10  ym  ZnO/ 


(lll)Si  samples  annealed  at  1000°C  for  1  hr  in  different 
atmospheres.  (A)  as  deposited  sample  (B)  annealed  in 
air.  Film_  peeled  off  dramatically  after  indentation  in- 
dicating  adhesion  is  very  poor.  (C)  annealed  in_yacuum_ 
at  low  oxygen  partial  pressure,  adhesion  is  imporved  as 
cocparad  to  (A) .  (D)  annealed  in  vacuum  at  extra  low 

oxygen  partial  pressure.  Both  the  hardness  and  the  ad¬ 
hesion  of  the  film  exhibited  the  greatest  improvement. 


INDENTATION  MODEL 


Fig.  10.  A  schematic  of  the  Indentation  Model  indicating  the 
residual  force  associated  with  the  plastic  zone  that 
acts  as  the  driving  force  for  lateral  crack  extension. 


Fig.  11.  Dependence  of  the  lateral  crack  length  on  the.  indenta-_ 

.  tion  load.  The  specimen  is  10ym  ZnO/lOOOA  Si02/(lll)Si, 
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DISCUSSION  AND  CONCLUSIONS  '  ; 

A  simple  indentation  test  for  measuring  the  adherence  of  thin 
or  thick  films  to  substrates  has  been  developed.  The  general  util¬ 
ity  of  this  test  as  a  simple  way  for  establishing  trends  in  the  ad¬ 
herence  of  film  systems,  achieved  by  employing  different  deposition 
or  post-deposition  techniques,  has  been  demonstrated.  This  was 
achieved  by  conducting  experiments  on  various  ZnO/Si  systems.  Im¬ 
portant  influences  on  adhesion  of  the  interlayer- -material  (Pt/Ti  -- 
being  preferred -to- Si02)  end  of  the  deposition  conditions- (columnar 
ZnO  grain  morphologies  yielding  superior  adhesion)  were  discussed  - 
and  correlated  with  the  acoustic  properties.  Additionally,  sub¬ 
stantial  effects. of  heat ^treatment  procedures  on  the  adherence  were 
detected  in  this  manner.  The  simplicity  of-the  indentation  tech — 
nique  renders  -it-  ideally  _  suitable  -  for— trend-determinat-ions-of— this 
type.  The  technique^can-be-used  with-  small-specimens— <a  f-ew-mill-i- 
meters  in  length)-j-rno- special-specimen  preparation-  is-required ,~a— 
standard  hardness^testing-machine— and— optical.— micros  cope-  -are— the — 
only  requirements,  the  adherence  measurements  are  reproducible, 

and-results-can— berobtainedrroutinelyT^quicklyr.and-wlthr-minimal - 

material— damage? — — -:rr- "  ■ — - - - 

The  development  of  indentation  methods  for  the  absolute  mea¬ 
surement  of  interface  fracture  resistance  has  been  studied.  An 
analytic  solution  for  well  developed  lateral  cracks  has  been  shown 
to  conform  quite  satisfactorily  with  data  obtained  on  a  ZnO/Si 
system.  Further  experimental  studies  are  needed  to  substantiate 
the  fracture  model.  Thereupon,  the  method  should  be  capable  of 
providing  direct  information  concerning  the  fracture  toughness  of 
interfaces.  j 
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INTRODUCTION 


It  is  now  well  established  that  zirconia  is  amongst  the  most 
refractory,  thermal  shock  and  corrosion-resistant  oxides.  At  atmos¬ 
pheric  pressure  zirconia  exists  as  three  modifications:^ 

Melt  2-8-°  -  Cubic  (I)  -  2-00-£  Tetragonal  (II)  -°  C  Monoclinic  (III). 

~95Q°C 

The  tetragonal  to  monoclinic  transformation  is  martensitic  and 

accompanied  by  a  3%  volume  increase,  both  of  which  are  physically 

delecterious  to  the -ceramic.  The  transformation  may  be  suppressed 

however,  by  small  additions  of  other  oxides,  notably  CaO,  ^2^3  an<^ 

( 2 ) 

MgO  .  The  cubic,  fluorite-type  solid  solution  then  persists  to 
room  temperature,  as  fully  stabilized  zirconia,  FSZ.  The  more  useful 
mechanical  properties  of  zirconia  arise  from  composite  microstructure 
of  monoclinic ^  or  tetragonal ^ 4 ^  particles  of  pure  zirconia  dis¬ 
persed  in  a  matrix  of  cubic  zirconia,  constituting  partially 
stabilized  zirconia,  PSZ. 

The  crystal  structures  of  each  phase  are  summarized  in  Fig.1. 

Due  to  the  similarity  of  tetragonal  and  monoclinic  face  centred  cells, 

three  different  lattice  correspondences  may  arise ^  as  illustrated 

in  Fig. 2.  If  a  right  hand  screw  convention  is  adopted  for  unit  cell 

axes,  lattice  correspondences  A,  B  or  C  may  be  defined,  depending 

on  which  monoclinic  axis  a  ,  b  or  c  respectively  is  parallel  to  c. . 

m  m  m  t 

(5) 

Bailey'  examined  the  II  to  III  transformation  by  transmission 
electron  microscopy  and  found  twinning  on  (110),  ( 1  To )  and  (100) 
monoclinic  planes,  as  well  as  shape  changes  and  an  orientation 
relation,  LC  C.  He  found  indirect  evidence  for  orientation  relations 
derived  from  all  three  lattice  correspondences.  A  martensitic 
mechanism  for  the  tetragonal  to  monoclinic  transformation  was  first 


reporJ  ed  by  Wolten^  who  found  symmetry  options  and  variants  for 

the  jrientation  relation  derived  from  lattice  correspondence  C. 

He  concluded  that  the  interface  had  a  monoclinic  habit  plane  (101) 

which  becomes  (101),  (110)  and  (Oil)  of  the  tetragonal  cell.  His 

X-ray  orientation  relations  imply  all  three  lattice  correspondences. 

As  reviewed  by  Subbarao  et  al.^,  several  other  determinations  of 

orientation  relations  essentially  arise  from  LC  B  or  LC  C. 

Bansal  and  Heuer v  '  '  examined  single  crystal  II  and  III 

transformations  by  transmission  electron  microscopy  and  X-ray 

precession  experiments.  They  found  that  an  orientation  relation 

derived  from  LC  B  occured  for  an  Ag  temperature  above  1000°C. 

As  is  the  temperature  at  which  the  reverse  (monoclinic  to  tetragonal) 

temperature  starts.  An  LC  C  (i.e.,"type  3")  orientation  resulted 

from  an  M  below  lOOO°c.  Habit  planes  of  the  tvpe  (671)  or  (761)  , 
s  mm 

being  approximately  7°  off  (110)m,  constituted  type  A,  and  usually 

occured  inside  the  crystal  as  opposed  to  surface  regions  where  type 

B  habit  planes  were  found.  Type  B  plates  were  also  lenticular  shaped 

with  a  midrib  close  to  (100)  or  (010)  ,  corresponding  to  internal 

m  m 

twinning  on  these  planes.  Coarse  (100)  twins  were  also  observed 

m 

in  plate-free  regions. 

( o  \ 

Martensite  calculations  were  performed  by  Eansal  and  Heuer 

(  Q  \ 

utilizing  the  algebraic  method  of  Bowles  and  MacKenzie  '  ,  for  eight 
slip  lattice  invariant  shear  (LIS)  systems.  They  were  done  for 
lattice  correspondence  B  and  C  as  summarized  below: 

LC  B  (110)  [00 11  (111)  [1  Tol  (100)[00l]  (100)  [pi  1] 

lc  c  dTo)[ooi]  dTo)[no]  (ioT)[oiol  ( 1 1 1 )  [i  TcT 
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Martensite  analvses  on  LC  C  by  (iTo)  [001]  slip  svsten  predicted 

(571>  habit  nlanes, while  (10Q)  habit  Dianes  were  nredicted  bv 
t  j  m  m 

si  p  on  ( 1  To)  0 10]  .  No  martensite  analvses  using  LIS  twinning 

svstens  have  as  vet  been  reported. 

In  the  details  of  the  Bansal  and  Heuerv  '  calculations,  a 

numerical  error  was  noted.  A  check  on  the  volume  increase  in  the 

transformation  using  the  algebraicallv  calculated  magnitudes  of 

the  princinle  distortions  had  values  of  1.3360%  and  1.3356%  for 

the  two  lattice  corresnondences  B  and  C,  respective lv .  The  volume 

change  from  the  unit  cell  dimensions  at  10OO°C  as  used  bv  Bansal 
(8) 

and  Heuer'  showed  that  the  volume  increase  is  actuallv  1.9030%. 
Errors  in  these  initial  stages  of  the  work  could  be  misleading. 

In  this  work  we  present  calculations  of  possible  crystallo¬ 
graphic  transformation  mechanisms  in  bulk,  unconstrained  zirconia 
crystals.  For  all  three  lattice  correspondences  slip  as  well  as 
one  twin  LIS  svstems  are  analvsed .  Habit  planes,  shape  changes 
and  orientation  relations  for  each  variant  are  predicted  in  terms 
of  the  tetragonal  structure.  The  shape  changes  associated  with 
each  mechanism  and  its  variants  are  compared. 


ANALYSIS 

Lattice  parameters  of  both  tetragonal  and  monoclinic  phases 
corresponded  to  a  transformation  temperature  of  950°C.  They  were 
calculated  from  exDerimental  parameters  and  thermal  expansion 
coefficients  determined  by  Patil  and  Subbarao ^ 10 ' 1 1 ^ 


and  gave  a 
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II  to  IHvolume  expansion  of  approximately  3.0%.  Thermal  expansion 

was  shown  to  be  anisotropic.  The  current  analysis  at 

950°C  gave  a  value  of  3.047%  which  was  consistent  with  the  volume 

increase  as  calculated  from  lattice  parameters.  Buljan,  McKinstry 

and  Stubican  determined  lattice  parameters  at  1000°C  which  were 

/  q  \ 

used  by  Bansal  and  Heuer '  '  and  also  in  this  analysis  to  compare 
the  effect  of  lattice  parameter  differences. 

(?) 

The  martensite  analysis  as  formulated  by  Bowles  and  MacKenzie 

was  programmed  *into  Fortran  IV  computer  language  by  Ledbetter  and 
(12) 

Wayman  in  PR0GRAM  MRTNST.  Calculations  were  done  in  an  ortho¬ 

normal  frame  of  reference  in  which  input  -data  also  needed  to  be 
expressed.  In  zirconia,  part  of  the  lattice  deformation  may  be 
visualized  as  a  simple  contraction  or  expansion  in  the  direction 
corresponding  to  bm  which  was  here  designated  as  the  n3  direction. 

The  remainder  of  the  lattice  deformation  occured  in  the  plane  which 
is  normal  to  this  direction.  The  determination  of  the  two  principle 
axes  and  distortions  in  this  plane  was  readily  accomplished  by 

standard  geometrical  theory  for  homogeneous  strain  in  two  dimen- 
(13) 

sions  .  The  details  are  given  in  Appendix  1.  These  calculations 

of  principle  strain  magnitudes  and  directions  were  confirmed  by  the 

( o ) 

general  method  of  Bowles  and  MacKenzie  .  Table  1  presents  lattice 
parameters  of  the  two  phases  at  950°C,  while  calculated  values  of 
the  deformations  are  given  in  Table  2. 

Important  features  arising  from  the  above  are  firstly,  that 
two  of  the  principle  deformation  directions  are  irrational.  Secondly, 
the  magnitude  of  varies  only  by  approx.  1%  from  unity  for  all 
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three  lattice  correspondences.  Finally,  the  values  of  -  3  which 
are  usually  taken  as  measures  of  the  strain  energy  favour  lattice 
correspondence  C,  but  not  greatly. 

Physically,  lattice  invariant  shears  occur  in  the  product, 
monoclinic  phase,  but  LIS  systems  were  specif ield  in  the  tetragonal 
phase  to  which  they  were  referred  from  the  monoclinic  cell,  by  in¬ 
spection  of  lattice  correspondence.  The  lattice  invariant  shear 
systems  investigated  were  derived  mainly  from  dislocations, 
stacking  faults  and  twinning  shears  found  in  the  monoclinic  phase 
after  transformation,  by  previous  experimental  studies.  For  lattice 
correspondence  B,  the  twin  system  (010)  [JOCQ^  was  analysed,  where 
the  direction  of  twinning  shear  was  specified  in  the  input  data 
and  treated  as  a  slip  system.  Mathematically,  computations  are 
identical  for  slip  and  twinning.  This  approach  was  necessary  as* 
PR0GRAM  MRTNST  was  unable  to  calculate  twin  directions  for  non- 
orthogonal  systems.  A  simple  flow  chart  describing  the  method  of 
calculations  is  drawn  in  Fig. 3. 


RESULTS 

It  was  sometimes  difficult  to  asses  whether  or  not  a  predicted 
magnitude  of  lattice  invariant  shear  was  physically  reasonable. 
Lattice  invariant  shears  by  twinning  posed  no  problem  as  magnitudes 
determine  relative  proportions  of  each  twin  in  the  product.  The 
lattice  invariant  shear  angle  must  therefore  be  less  than  the  angle 
of  twinning  shear.  For  slip  or  stacking  fault  shears  the  shear  angle 
associated  with  an  individual  dislocation  or  stacking  fault  was 


-  t- 
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readily  determined  and  in  this  work  any  LIS  angle  greater  than  half 
the  maximum  value  was  taken  to  indicate  an  unreasonable  defect 
density.  Thus,  a  LIS  of  20°  C  was  taken  as  an  upper  limit  for 
"  reasonable  "  solutions.  Shape  strains  (  m1  )  greater  than  20% 
were  considered  unlikely  and  hence  not  reported. 

( 8 ) 

Table  3  compares  the  calculations  of  Bansal  and  Heuer  with 
those  computed  by  the  current  method,  using  the  same 
lattice  parameters  at  1000°  C.  It  is  seen  that  although  there  is 
general  agreement  there  is  a  scatter  of  habit  planes,  and  variation 
in  the  magnitude  of  shape  strain  by  approximately  10%.  Table  4 
summarizes  habit  planes  and  shape  strains  (  ,  d1  )  predicted 

at  950°  C  for  each  mechanism.  The  combination  of  all  three  lattice 
correspondences  with  LIS  systems  yielded  32  distinct  reasonable 
solutions.  Comparison  of  above  results  indicate  that  quantitative 
predictions  of  martensite  mechanisms  depend  on  lattice  parameters 
_of  both  phases  and  their  values  at  the  same  transformation 
temperature. 

Orientation  relations  were  determined  by  multiplying  vectors 

(12) 

and  plane  normals  by  the  total  strain  matrix  .  A  vector  V  was 

rotated  to  become  a  vector  V'  both  in  the  same  cubic  basis.  The 
mutual  angle  of  rotation  was  then  calculated  for  that  crystal  system. 
With  reference  to  the  operating  lattice  correspondence,  the  product 
vectors  and  plane  normals  were  then  re-labelled  according  to  the 
monoclinic  axes.  The  details  of  predicted  orientation  relations 
are  presented  in  Table  5. 
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Twinning  on  {010}<100>  as  lattice  invariant  shear  resulted  in 
four  symmetry-equivalent  solutions  for  each  twin  plane.  The  s  jear 
angle  associated  with  twinning  in  this  way  was  8.78°.  The  pro¬ 
portion  of  each  twin  in  the  product  was  predicted  to  be  0.409  : 

0.591  and  0.617  :  0.383  (for  LIS  angle  =  3.61°  and  5.44°,  respec¬ 
tively).  These  ratios  were  close  to  2:3  and  3:2.  Habit  planes 
approximated  to  {T21}^and  the  total  shape  change  was  significantly 
less  than  the  shape  change  of  the  lattice.  The  four  different  (but 
equivalent)  orientation  relations  for  each  twinning  LIS  system 
could  each  be  described  by: 

I  001 1 t  A|0l0|m  =  1.18°,  { 1 00} ^  A(l00)m  =  1.14°. 

DISCUSSION 

The  calculations  presented  above  are  quantitative  analyses 

of  possible  crystallographic  transformation  mechanism  operating  at 

a  temperature  of  950°C.  Comparison  of  results  based  on  different 

lattice  parameters  and  values  of  principle  strains  indicates  that 

such  predictions  are  dependent  on  lattice  parameters  of 

both  parent  and  product  phases,  and  that  their  values  at  the  same 

transformation  temperature  should  be  used.  Since  lattice  parameters 

(14) 

are  affected  by  solute  content  which  also  affects  the  M 

5 

temperature,  different  mechanisms  may  thus  be  favoured  or  altered 
under  different  conditions. 

For  all  three  lattice  correspondences  the  outstanding  feature 
of  the  results  obtained  with  stacking  fault  or  slip  modes  (apart 

from  {010}<101>)  is  the  small  (<i°)  Df  LIs  shear  required. 

The  predicted  crystallography  is  therefore  very  similar  for  a 
number  of  possible  LIS  systems.  Individual  orientation  relations 


therefore  night  only  be  resolved  with  precise  transmission  electron 
microscopy  since  particular  reflections  from  one  variant  would 
occur  within  2-3  degrees  of  eauivalent  reflection  from  ther 
variants . 

The  shape  strain  mav  be  resolved  into  components  parallel 
and  perpendicular  to  habit  nlanes  (Fig. 4).^^  Fron  Table  3  it  is 
seen  that  most  of  the  shape  strain  (m^ )  is  resolved  parallel 

to  the  habit  plane.  Table  4  shows  that  slip  LIS  mechansims  produce 

i 

shaoe  strains  of  11%  -  17%,  while  the  (010) twin  svstem 

tet. 

produced  strains  of  5%  per  variant. 

As  mentioned  in  the  introduction,  there  is  some  experimental 
evidence  for  the  occurrence  of  both  lattice  correspondence  S  and  C. 
There  is  no  apparent  reason  why  lattice  correspondence  *  should 
be  less  probable.  From  Table  5  it  is  seen  that  che  rotational 
component  of  the  total  mechanism  is  small  in  all  instances.  Hence 
Qucisj]  t  [axis]  m  does  implv  that  [«is]  t  becomes  JaxisJm  in 
the  transformation. 

The  (010) t  twin  system,  bv  lattice  correspondence  3  (Fig. 2) 

becomes  (001)  or  4100)  when  eouivalent  a.  and  b.  axes  are 
m  m  t  t 

interchanged.  (100)  and  iliol-  def ormation  twinnina  was  found,  bv 

m  l  j  m 

3ailey^  in  thin  TEM  specimens.  Kriven^1f^  also  found  the  same 

twin  systems  in  small  included  ZrO£  Darticles  in  AljO^-ZrO^  and 

pure  Zr©2  (42%  volume  fraction  tetragonal)  ceramics. 

The  calculations  presented,  here  however,  imply  that  the 

(100)  and  (001)  can  also  act  as  lattice  invariant  shear  svstems 
mm 

in  a  martensitic  mechanism,  giving  rise  to  a  nacrosconic  shape 


change  and  habit  plane.  This  finding  may  be  relevant  to  understanding 

how  ■  stress-induced  martensitic  transformation-,  in  included 

zirconia  particles  are  able  to  toughen  comrjsite  ceramics. 

The  fact  that  (100)  mav  act  as  an  LIS  tvin  system  further 

m 

suggests  that  martensite  calculations  be  performed  for  all  the  (100) r 

and  {110}  twin  systems  in  combination  with  the  three  lattice 
m 

correspondences  .  Such  calculations  would  require  modification  of 
the  method  here  presented. 

CONCLUSION 

In  conclusion,  therefore,  we  have  quantitatively  analysed 
some  martensitic  transformation  mechanisms  in  zirconia.  The  analyses 
presented  were  purely  geometrical  and  based  on  the  invariant  plane 
strain  criterion  of  a  martensitic  transformation  in  a  bulk  crystal. 
Lattice  parameters  were  corrected  for  both  tetragonal  and  monoclinic 
phases  to  the  same  trans formation  temperature  of  950°C.  Three 
principle  strain  axes  were  determined  in  the  monoclinic  unit  cell, 
and  the  different  tetragonal  LIS  systems  were  obtained  from  the 
monoclinic  cell  by  inspection,  with  reference  to  three  lattice 
correspondences.  The  computed  martensite  analyses  then  yielded 
solutions  which  predicted  habit  planes,  shape  strains  related  to 
the  habit  planes  and  orientation  relations  for  each  variant^ given 
with  respect  to  the  tetragonal  structure.  Most  of  the  LIS  slip 
systems  had  shape  strains  of  the  order  of  11%  -  17%,  while  one 
twin  LIS  system  had  5%  shape  strains  resolved  essentially  parallel 
to  the  habit  plane. 

Thus,  the  work  presented  here  has  shown  that  (1)  the  (100) 

m 

or  (OO’Ujn  twin  systems  may  act  as  lattice  invariant  shear  systems, 

giving  a  martensitic  solution  with  macroscopic  shape  chances 

and  habit  planes  etc.  (2)  It  illustrates  the  different  shape  changes 
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associated  with  each  mechanism  and  variant.  (3)  It  indicates  the 
importance  of  using  lattice  parameters  corrected  to  the  same 
transformation  temperature. 
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Table  1  .  Lattice  Parameters 
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Table  4. 


PREDICTIONS  OF  HABIT  PLANES  AND  SHAPE  STRAINS. 


Lattice 

INVARIANT 

Shear  System 


Habit  Plane 
indices  (PI) 


Lattice  Correspondence :  A 


Direction  of 
Shape  Strain 
D1 


(010)  [10J 


(110)  0 1 o) 


*?(161) 


~(212) 


(010) 


0.638 

-0.281 

0.701 

-[212] 

.0.126 

-0.967 

0.217 

^[l82] 

-0.002 

-0.266 

-0.942 

^\02T\ 


/-0 .003 
-0.346 
^-0.960^, 
~(013) 


-0.01 2 
-0.984 
0.171 

~'[06lJ 


(ioi)  [ToJ 


~(027) 


0.000 

0.995 

-0.094 

~[oiq] 


(010) 


-0.012 

-0.203 

0.957 

,v[pT5] 


Magnitude  of 
Shape  Strain 
Ml 


0.122 


0.122 


0.173 


0.173 


0.164 


0.164 


LATTICE 
INVARIANT 
SHEAR  SYSTEM 


-  15  - 


DIRECTION  OF 
SHAPE  STRAIN 
D1 


MAGNITUDE  OF 
SHAPE  STRAIN 
Ml 


HABIT  PLANE 
INDICES  (PI) 


(in)  [no] 


^-0.057^ 
-0.996 
^  0.075^ 

Moio) 

(  0.025N 
-0.324 
> -0.968, 
M013) 

/  0.661N 
0.354 
0.678, 
^(212) 


'0.839 
-0.515 
k  0.180 


M531) 


Lattice  Correspondence:  B 
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(Oil)  [Oil] 


0.615 
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0.730 


~(212) 

^-0.1 03N 
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Figure  Captions 

Figure  1.  Crystal  structures  of  monoclinic  and  tetragonal  zirconia 

with  lattice  parameters  at  950°C  for  both  phases.  Monoclin 
twin  planes(5)  are  (100), (110)  and  (110)  while  slip  occurs 
on (010) . 

Figure  2.  Lattice  correspondences  for  zirconia. 

Figure  3.  Flow  chart  illustrating  method  of  calculations. 

Figure  4.  The  relation  between  shape  deformation,  angle  of  shear 
and  habit  plane. 


Figure5.  The  determination  of  magnitudes  and  directions  of 
principle  strains  by  a  coordinate  geometry  method. 
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Directions  .nd  magnitudes  of  the  principal  distortions  associated 
with  the  three  possible  lattice  correspondences  between  tetragonal  and 
monoclinic  zircor.ia  have  been  determined  by  the  method  of  analytical 
geometry . 

The  problem  is  two  dimensional  as  one  principal  axis  is  obtained 
directly  by  inspection  f  n3//bm) .  The  remaining  two  principal 

distortions  must  lie  in  the  plane  of  which  the  direction  is  the  normal, 


The  theory  of  finite  homogeneous  strain  in  two  dimensions  (e.g.  Jaeger 
gives  a  simple  method  for  the  location  of  these  two 

vectors: 

The  changes  in  the  relevant  tetragonal  plane  can  be  factorised 
into  (a)  a  change  of  dimensions  (expansion/ contraction)  and  (b)  a 
change  of  shape  (simple  shear) .  This  is  illustrated  in  figure  5^ 
Consider  a  point  (x,y)  in  the  tetragonal  plane: 


03) 


(  i) 


x 

m 

x  ■>  —  x, 
Xt 


y  sin  S 

y  -►  - - y  (by  expansion) 


(x  -*•  ax)  Cy  -*•  dy) 

(ii)  ax  ax  +  tan  (90  -  g)dy,  dy  unchanged  (by  shear) 

(ax  -*■  ax  +  by)  (dy  -►  dy). 

Hence  x  -*■  x  +  tan  (90  -  S)  y  sin  6 

a  ® _ 

—  x  - -  Y 

xfc  y*. 


(x  *  ax  ♦  by) 


and 


'm 


sin  3 


:•/  *  dv) 


» W-. it  .,<V 


\  wi  .*i 
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From  two  dimensional  strain  theory,  for  x  -*■  ax  +  by  and 
y  cx  +  dy  directions  of  the  principal  distortion  axes  are  given  by  : 


tan  2  c 


2  (ah  +  cd)  2ab 


(since  c  =  o  in  this  case) . 

(-  is  the  angle  made  with  the  x  axis ,  and  the  required  angles  are 
a,  a  +  It/2) . 


Magnitudes 

(A, 

3)  are 

given  by: 

(A  +  B)2  = 

(a 

+  a)  + 

(b  - 

c)2  = 

(a  + 

d)2  +  b 

(A  -  B)2  - 

(a 

-  d)2  + 

(b  + 

c)2  = 

(a  - 

d)  2  ”+  b 

Directions  of  the  principal  distortion  axes  in  the  monoclinic 
phase  can  also  be  determined: 


>1 


A 


tan  2c' 


2  (ac  +  bd) 


2bd 


2  ^  w2  2  2 
a  +  b  -  c  -  d 


2  __  w2  ,2 
a  +  b  -  d 


(a',  a'  +  H/2  are  the  angles  formed  by  the  principal  distortion  axes  and 
the  x  axis) . 

Rotation  of  the  principal  axes  during  the  deformation  is  therefore 
(a'  -  a) . 

These  formulae  can  be  applied  to  tetragonal  £  monoclinic, 
orthorhombic  monoclinic  and  cubic  +-  monoclinic  changes. 
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ABSTRACT 

The  mechanisms  of  material  removal  during  wear  have  been  examined. 
Primary  emphasis  has  been  placed  on  the  lateral  fracture  mechanism,  in  an 
attempt  to  elucidate  both  its  realm  of  importance  and  the  concomitant  ma¬ 
terial  removal  rates.  Correlations  with  existent  data  are  shown  to  be 
encouraging.  The  enhancement  of  material  removal  by  heat  generation  in 
the  plastic  zone  has  been  examined  and  shown  to  be  generally  of  minor  sig 
nificance.  Some  considerations  of  material  removal  below  the  fracture 
threshold  have  also  been  presented. 
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I.  INTRODUCTION 


The  wear  of  ceramics  has  been  the  subject  of  extensive  empirical 


investigation  .  For  example,  material  removal  rates  have  been  measured 


on  various  materials  as  functions  of  the  normal  force,  the  horizontal 


(frictional)  force,  the  wheel  speed,  etc.-  However,  the  detailed  mechanisms 


of  material  removal  have  only  been  cursorily  explored.  Yet,  some  under¬ 


standing  of  mechanisms  is  an  essential  prerequisite  to  the  development  of 


optimum  wear  characteristics  (especially  for  complex  processes  such  as 


abrasive  wear).  The  intent  of  the  present  paper  is  to  examine  possible 


mechanisms  of  material  removal. 


Several  investigators'^  have  recently  recognized  the  important  role 
of  lateral  fracture  (fig.  la)  in  material  removal.  This  mechanism  of  wear 


is  explored  in  some  detail,  to  determine  both  its  realm  of  importance  and 


the  parameters  that  might  influence  the  removal  rate. 

5  6 

Lateral  cracking  occurs  above  a  threshold  force  ’.  This  threshold 


coincides  with  a  transition  from  a  relatively  severe  lateral  fracture  domi¬ 


nated  material  removal  process  at  higher  force  levels  to  other  less  damaging 


wear  mechanisms.  The  fracture  threshold  is  thus  an  important  concept  in  the 


study  of  wear,  and  the  nature  of  the  threshold  should  be  afforded  appreciable 


emphasis.  The  material  removal  mechanism  below  the  fracture  threshold  is 


presently  uncertain.  An  operative  (although  not  necessarily  exclusive) 
mechanism  involves  the  plastic  cutting  that  occurs  during  the  grooving 
process  (fig.  lb). 


The  material  removal  rate  above  the  fracture  threshold  evidently  de¬ 
pends  on  the  extent  of  the  lateral  fracture;  An  analysis  of  lateral  crack 


J 


extension  will  thus  be  used  to  predict  wear  rates  as  a  function  of  the 
dominant  material  variables.  Data  for  a  range  o-  materials  will  be  corre¬ 
lated  with  the  predictions. 

The  plastic  flow  associated  with  surface  penetration  must  lead  to  a 
local  temperature  increase.  This  temperature  rise  inevitably  results  in 
thermal  stresses,  that  enhance  the  extension  of  the  lateral  cracks.  The 
importance  of  thermal  effects  vis-a-vis  residual  stress  motivated  lateral 

i 

fracture  is  examined  in  the  third  section  of  the  paper.  Then,  some  consid¬ 
erations  of  material  removal  below  the  fracture  threshold  are  presented. 
Finally,  some  general  implications  for  the  development  of  wear  resistant 
materials  are  discussed  from  the  perspective  provided  by  the  present  analy¬ 
ses. 

2.  THE  LATERAL  FRACTURE  MODEL 

2.1  The  Threshold 

The  fracture  threshold  is  a  phenomenon  with  major  significance  for 
the  control  of  wear,  because  it  signifies  a  change  in  wear  mechanism  and 
a  concommitant  order-of-magnitude  increment  in  wear  rate.  The  general 
characteristics  of  fracture  thresholds,  at  sites  of  plastic  penetration, 
are  reasonably  well  comprehended5,6,7,  although  detailed  quantitative  pre¬ 
dictions  are  subject  to  some  uncertainty.  Of  particular  interest  to  the 
wear  mechanism  transition  is  the  minimum  threshold  load  P*  .  This  threshold 
load  pertains  in  the  presence  of  an  adequate  spatial  density  of  crack  pre¬ 
cursors. in  the  vicinity  of  the  plastic  zone6;  a  situation  that  is  expected 

s 

to  prevail  near  surfaces  subject  to  wear  damage.  The  minimum  threshold  is 
dictated  by  the  hardness  H  ,  toughness  Kc  and  elastic  modulus  E  of  the 
material,  in  accord  with  the  relation 
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(1) 


where  the  function  f(£/H)  and  the  constant  c  depend  upon  the  type  of 

crack.  For  lateral  cracks,  recent  analysis^  indicates  that  f  is  insensi- 

5 

tive  to  E/H  and  that  ?f  '5=  2  x  10  .  Threshold  predictions  for  lateral 
fracture  in  several  important  material  systems  are  summarised  in  Table  I. 
The  threshold  load  levels  are  relatively  low,  but  loads  of  this  magnitude 
can  exist  at  individual  wear  particles:  depending  upon  the  nature  of  the 
wear  medium,  the  contact  points  and  the  constraints  imposed  on  the  wear 
system.  There  have  been  few  experimental  studies  of  lateral  fracture  ini¬ 
tiation  and  validation  of  the  predicted  threshold  characteristics  is  not 
possible  at  this  juncture.  Such  studies  are  now  in  progress.  However,  it 
is  noted  that  predictions  of  the  radial  fracture  threshold  using  the  same 
analytic  procedures  are  in  good  accord  with  the  more  comprehensive  experi¬ 
mental  observations  that  exist  for  this  system6,7. 


2.2  The  Lateral  Crack  Extension 

Lateral  cracks  develop  exclusively  due  to  the  residual  stress  created 

g 

by  indentation0.  The  maximum  extension  of  the  crack  is  thus  realised  when 

the  penetrating  particle  is  removed.  The  cracking  behavior  can  be  analysed 

by  regarding  the  plastic  zone  as  a  precompressed  spring  which  motives  crack 
9  10 

extension  *  .  The  crack  advance  is  then  dictated  by  the  compliance  of  the 

material  above  the  crack,  treated  as  an  elastic  cantilever.  The  specific 
crack  extension  depends,  in  part,  upon  the  morphology  of  the  plastic 
penetration.  For  axisymmetric  indentation  (most  pertinent  to  intermittent 
particle  contact)  the  lateral  crack  length  c  is  given  by^; 
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(2) 


^  ★ 
where  P  is  the  peak  load  during  particle  penetration,  PQ  is  an  apparent 

threshold  load  for  fully  developed  cracking  and  c*-  is  given  by; 


where  is  a  material  independent  constant  that  depends  on  the  particle 
shape  and  can  be  determined  by  calibration  on  a  material  with  well  charac¬ 
terised  fracture  properties  (e.g.  glass).  For  a  linear  plastic  groove  (as 
might  pertain  for  continuous  particle  penetration)  the  equivalent  result 
(Appendix,  eqn  A-9)  is; 


where  is  another  material  independent  constant  The  only 

significant  difference  between  the  two  solutions  for  cL  occurs  in  the 
( E/H )  exponent.  But,  since  the  range  of  E/H  values  for  most  ceramics  is 
less  than  an  order  of  magnitude,  the  influence  of  the  grinding  morphology 
upon  crack  extension  is  considered  to  be  of  minor  importance. 

For  present  purposes,  the  indentation  forces  are  considered  to  be 
sufficiently  in  excess  of  the  apparent  threshold  to  permit  -the  direct  use 
of  eqn.  (3)  for  the  prediction  of  crack  extension.  The  wear  behavior  under 
these  conditions  will  exhibit  relatively  simple  functional  dependencies  on 

m 

the  dominant  material  variables  (as  described  in  the  subsequent  section).  The 
predictions  can  thereupon  be  straightforwardly  correlated  with  experimental 
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results.  The  cracking  behavior  at  intermediate  loads  is  less  extensive 
than  that  predicted  from  eqn.  (3)  and  hence,  the  wear  predictions  based 
on  eqn.  (3)  represent  an  upper  bound  wear  rate.  A  more  complete  analysis 
of  crack  extension  at  intermediate  loads  can  be  conducted  if  experimental 
correlations  indicate  that  some  important  trends  are  inadequately  predicted 
by  neglect  of  the  apparent  threshold. 

The  above  solutions  for  the  lateral  crack  length  pertain  to  conditions 
wherein  all  of  the  plastically  deformed  material  is  retained  within  the 
plastic  zone  and  thus  represent  an  upper  bound  crack  driving  force.  When 
some  of  the  plastically  deformed  material  is  removed  by  plastic  cutting, 
the  residual  forces  will  be  reduced  and  hence,  the  lateral  crack  extension 

I 

diminished.  It  is  not  known  at  this  juncture  whether  the  relative  allevia¬ 
tion  of  the  residual  forces  by  cutting  exhibits  important  material  depen¬ 
dencies.  Such  effects  would  emerge  from  inadequate  correlations  of  experi¬ 
mental  results  (between,  materials)  with  models  based  on  eqn.  (3). 

2.3  The  Material  Removal  Rate 

The  material  removal  analysis  considers  the  extension  of  the  lateral 
cracks  normal  to  the  motion  of  the  abrading  particle  (fig.  la),  since  this 
distance  defines  the  width  of  the  potential  material  removal  zone.  The 
situation  to  be  modelled  is  defined  in  fig.  la;  it  consists  of  an  abrasive 
particle  subject  to  a  normal  force  Pn  ,  moving  with  a  prescribed  lateral 
velocity  v^  . 

The  volume  Vi  of  material  removed  by  the  passage  of  each  abrasive 

★ 

particle  that  exceeds  the  fracture  threshold  force  P  is  determined  by 
the  extent  c^  and  the  depth  h^  of  the  lateral  fracture.  The  upper 
bound  volume  removal  is; 
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where  4^  is  the  distance  of  motion, 
scales  with  the  plastic  zone  radius  b 
plastic  penetration  radius,  a  ,  by^ 


(4) 

The  lateral  fracture  depth  typically 
which,  in  turn,  relates  to  the 


b  ^aCE/H)275 


(5) 


and  the  indentation  radius  is  related  to  the  hardness  and  the  peak  normal 
penetration  force  by 

a2  %  (Pn/H)  (g) 

Substituting  eqns.  (3),  (5)  and  (6)  into  eqn„(4),  the  volume  removed  be¬ 
comes  ; 

p9/8 

V°3  l"z  Hs/B  <E/H)VS  <7> 

where  is  a  material  independent  constant. 

Implicit  in  the  preceding  analysis  is  the  dominant  influence  of  the 
normal  force  Pn  upon  the  material  removal  process.  The  lateral  force 
P£  =  uPn  (presumed  to  be  dictated  by  the  resistance  of  the  material  to 
plastic  grooving,  i.e.  ductile  cutting  resistance)  is  considered  incidental 
to  material  removal  by  lateral  fracture.  Some  substantiation  of  this  premise 
will  emerge  later  from  experimental  correlations.  At  this  juncture,  it  is 
noted  that,  by  analogy  with  the  equivalent  elastic  problem^,  the  lateral 
force  may  tend  to  rotate  the  stress  field  in  the  direction  of  P£  ,  in  the 
sense  that  the  tension  is  enhanced  behind  the  contact  zone  and  suppressed 
ahead  of  the  contact.  However,  even  if  this  occurs,  there  should  be  no 


appreciable  influence  on  the  stresses  contained  in  the  plane  normal  to 
Pn  ;  stresses  which  have  the  primary  effect  on  lateral  crack  extension. 

Hence  under  quasi-static  conditions,  the  width  ofthe  lateral  fracture  zone 
should  berelatively  insensitive  to  the  magnitude  of  the  lateral  force. 

At  relatively  high  lateral  velocities,  local  instabilities  tend  to 
develop.  These  instabilities  presumably  occur  following  a  fracturing  event 
and  result  in  an  abrubt  increase  in  the  local  particle  velocity,  as  mani¬ 
fest  in  the  lateral  force  oscillations  typically  observed  in  studies  of 
the  motion  of  abrading  particles.  The  impact  events  that  result  from  the 
instabilities  are  likely  to  have  impulse  components  both  normal  and  parallel 
to  the  surface.  The  normal  component  affects  the  penetration,  and  thus 
has  the  more  significant  influence  on  the  material  removal  process.  The 
normal  velocity,  and  the  consequent  magnitude  of  the  lateral  fracture,  may 
depend  upon  the  amplitude  of  the  instability  and  hence,  indirectly  upon  the 
friction  coefficient.  These  complexities  may  be  of  secondary  significance, 
and  the  existence  of  dynamic  effects  should  be  sought  in  experimental  results 
before  embarking  on  an  analytic  study  of  the  role  of  friction  in  the  wear 
process. 


2,4  Correlations  With  Experiment 

Experimental  results  suitable  for  correlation  with  the  present  wear 
predictions  include  measurements  of  the  influence  of  the  normal  force  on 
the  wear  rate  (observed  for  a  specific  material)  and,  more  significantly, 
studies  of  the  trends  in  the  wear  rate  among  different  materials.  The  in¬ 
fluence  of  material  properties  on  wear  can  be  compared  with  the  lateral 
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fracture  model  either  from  measurements  of  the  wear  forces  at  a  specified 
(constant)  volume  rate  of  removal,  or  from  removal  rates  at  a  constant 
level  of  applied  force. 

Direct  measurements  of  the  lateral  crack  width  have  been  conducted  on 
MgZn  ferrite  as  a  function  of  the  normal  force  .  These  results  are  consis¬ 
tent  with  the  force  exponent  of  5/8  predicted  for  the  lateral  crack  exten¬ 
sion  (fig.  2).  The  influence  of  the  normal  force  on  the  volume  removal 

12 

rate  has  been  studied  for  several  materials  by  Wilshaw  and  Hartley  . 

These  results  (fig.  3)  indicate  a  force  exponent  that  exceeds  unity,  as 
anticipated  by  the  present  analysis  (eqn.  (7)).  The  experimentally  observed 
exponent  ranges  between  1.07  and  1.3:  values  which  compare  quite  favorably 
with  the  predicted  value  of  1.13. 

The  most  comprehensive  material  removal  data  suitable  for  comparison 
with  the  present  predictions  are  grinding  data  obtained  by  Koepke  and 
Stokes  .  These  data  comprise  force  measurements  determined  on  different 
materials  for  constant  volume  removal  conditions.  For  up-grinding  condi¬ 
tions,  which  best  simulate  the  behavior  analysed  in  the  preceding  section 
(because  fracture  chips  are  directly  ejected  from  the  system),  the  force 
data  for  a  wide  range  of  materials  (fig.  4)  are  in  remarkable  accord  with 
the  predicted  material  dependence,  K  H"5/^8  .  The  superior  correlation 
exists  for  the  dependence  on  the  normal  force;  while,  a  more  modest  corre¬ 
spondence  obtains  for  the  lateral  force.  These  relative  correlations  tend 
to  substantiate  that  the  development  of  a  lateral  force  is  incidental  to 
the  material  removal  process.  The  one  material  that  deviates  from  the 
correlation  is  lucalox  AlgO^.  This  is  a  relatively  coarse  grained  A^O^ 
that  also  exhibits  unusual  indentation  fracture  characteristics14.  Under 
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incintation  conditions  extensive  grain  boundary  separation  occurs,  rather 
than  the  formation  of  discrete  crack  systems.  Grinding  in  this  (and 
similar)  materials  could  thus  procede  by  grain  ejection  and  hence,  the 
material  removal  would  be  expected  to  exhibit  a  different  dependence  on 
material  properties  than  anticipated  by  the  lateral  fracture  model. 

Machining  data  at  constant  force  are  also  available  for  a  range  of 
1 5 

materials  .  However,  the  data  exhibit  an  appreciable  grain  size  depend¬ 
ence.  This  behavior  presumably  arises  because  the  lateral  fracture  in  the 
coarser  grained  materials  extends  over  a  relatively  small  number  of  grains 
and  hence,  the  Kc  and  H  values  required  for  the  volume  removal  predic¬ 
tions  are  smaller  than  the  isotropic  polycrystalline  values.  It  is  expe¬ 
dient,  therefore,  to  conduct  the  comparison  at  the  small  extreme  of  grain 
size  for  which  available  polycrystall ine  material  parameters  are  most  likely 
to  pertain.  The  resultant  correlation  is  summarised  in  fig.  5.  The  correct 
trends  are  predicted  by  the  material  parameter,  Kc“1//2  H"5//8  ,  although  the 
correlation  is  not  as  satisfying  as  that  evident  for  the  force  measurements. 

The  relatively  consistent  trends  in  both  the  material  removal  rate  and 
the  grinding  forces  with  the  material  parameters  suggested  by  the  simple 
lateral  crack  analysis  indicates  that  most  of  the  proposed  complexities,  if 
they  occur v.  do  not  impose  additional  material  dependent  influences  upon  the 
wear  process"*".  Dynamic  effects  and  residual  force  relaxation  by  plastic  re¬ 
moval  can  presumably,  therefore,  be  considered  to  introduce  only  small 


^Recent  abrasive  wear  studies  reported  by  Moore  and  King  8  indicate  a  poor 
experimental  correlation  based  on  K  5/8^1 /2#  However,  the  toughness  values 
determined  in  the  analysis  are  substantially  inaccurate  (e.g.  3.0  MPar^m  for  _ 
soda  lime  glass  compared  with  a  q_enerally  accepted  value  of  0.7  MPav^n:  9.1MPa/m 
for  SUN.  compared  with  4-5  MPa/mT  etc.)  Reanalysis  of  their  results  is  thus 
needed  before  embarking  on  an  assessment  of  the  inadequacies  of  the  simple 
lateral  fracture  model. 
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perturbations  upon  the  wear  characteristics  anticipated  from  the  quas"'- 
static  lateral  fracture  analysis,  at  least  for  the  test  conditions  used 
and  for  materials  within  the  range  of  properties  encompassed  by  the  existent 
experimental  results. 

3.  THERMAL  EFFECTS 

Heat  is  invariably  generated  during  a  grinding  operation.  This  heat 

derives  from  the  plastic  work  expended  by  the  penetration  of  the  abrading 

particle.  The  heat  is  generated  within  the  plastic  zone  and  leads  to  an 

expansion  of  this  zone  vis-a-vis  the  surrounding  elastic  material.  The 

zone  expansion  is  directly  analagous  to  the  effective  expansion  induced 

9  10 

by  plastic  penetration  ’  .  Hence,  heat  generation  can  be  considered  to 

yield  a  lateral  crack  driving  force  that  simply  superimposes  on  that  asso- 
ciated  with  the  residual  indentation  field.  The  ratio  of  crack  driving 
forces  attributed  to  these  two  sources  of  zone  expansion  will  determine 
their  relative  influence  upon  the  wear  rate,  by  lateral  fracture.  Addi¬ 
tionally,  it  is. noted  that  the  zone  expansion  induced  by  heating  creates 
shear  stresses  outside  the  isothermal  plastic  zone.  These  stresses  tend 
to  enlarge  the  plastic  zone,  increasing  its  volume  from  V  to  V*.  This 
zone  enlargement  also  contributes  to  the  crack  driving  force  ratio. 

The  work  done  in  creating  an  indentation  of  volume  AV  is  simply, 

HAV®.  If  this  work  is  dissipated  as  heat  within  the  plastic  zone,  the 
average  temperature  rise  0  for  adiabatic  conditions  is 

9  =  HAV/pcpV*  (8) 


-Mdb&i 
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where  p  is  the  density  -id  cp  is  the  specific  heat.  The  local  tempera¬ 
ture  rise  AT  at  a  radial  location  r  is  related  to  the  average  tempera¬ 
ture  6  by17; 

AT  =  30Jtn(b/r)  (9) 

The  adiabatic  unconstrained  volume  strain  is  thus+. 


where  a  is  the  linear  thermal  expansion  coefficient.  Since  AV/V  is  the 
unconstrained  volume  strain  of  isothermal  plastic  penetration,  the  adiabatic 
crack  driving  force  ratio  becomes; 

R  =  eT/(A V/V)  *  (3aH/pcp)(V*/V)  (11) 

Inserting  some  typical  values  for  the  material  parameters  (hMOGPa, 

a  ^5  x  10”®K"\  p  %3  x  103Kgm"3,  cp  ,v600JKg‘1;<“1 )  indicates  that  R  is 
★ 

M3.1  when  V  *5  V.  Usually,  therefore,  thermal  effects  would  be  expected  to 
exert  a  minor  influence  on  the  wear  rate.  A  significant  thermal  effect  may 
emerge,  however,  in  hard  materials  that  exhibit  a  larce  thermal  expansion 
coefficient  and  a  low  specific  heat. 

4.  SUB  THRESHOLD  HEAR  RATES 

Below  the  lateral  fracture  threshold,  material  is  removed  by  a  plastic 
cutting  process  (fig.  lb).  The  material  removal  rate  is  presumably  dictated 
by  the  plastic  penetration  and  should  thus  be  uniquely  related  to  the  material 

'It  is  noted  that  the  assumption  of  a  uniform  temperature  rise  6  within 
the  plastic  zone  yields  an  identical  result. 
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hardness,  the  normal  force  and  the  particle  geometr  ; 

A 

P  4. 

Vi  "  =ot(f/2) 


(1:2) 


where  4*  is  the  included  angle  of  the  penetrating  particle.  The  volume 
removal  rate  induced  by  plastic  cutting  is  about  an  order  of  magnitude 
smaller  than  that  associated  with  the  lateral  fracture  mechanism1®.  Appre¬ 
ciable  wear  resistance  is  thus  experienced  with  this  regime.  The  choice 
of  wear  conditions  that  exclude  lateral  fracture  is  evidently  a  highly  de¬ 
sirable  material/design  objective. 

5.  IMPLICATIONS  AND  CONCLUSIONS 

The  analyses  of  fracture  dominated  wear  processes  in  brittle  solids 
indicates* a  reasonably  good  correlation  with  available  experimental  data. 

This  correlation  substantiates  the  strong  importance  of  both  the  fracture 
toughness  and  the  hardness  upon  the  wear  rate:  high  toughness  and  high 
hardness  both  being  desirable  for  optimum  wear  resistance.  A  small  influ¬ 
ence  of  the  material's  thermal  parameters  is  also  anticipated,  but  only  for 
hard  materials  with  a  large  thermal  expansion  coefficient  or  a  small  specific 
heat. 

A  major  diminution  in  wear  rate  is  associated  with  force  levels  below 
the  fracture  threshold.  The  force  at  the  fracture  threshold  is  strongly 
influenced  by  both  the  hardness  and  the  toughness,  such  that  fracture  can 
be  suppressed  by  a  high  toughness  or  a  low  hardness.  Below  this  threshold, 
the  material  removal  rate  is  substantially  diminished  and  depends  principally 


upon  the  hardness. 

In  general,  materials  with  an  optimum  combine tion  of  high  toughness 
and  high  hardness  will  be  the  most  wear  resistant.  However,  there  is  some 
ambiguity  associated  with  the  role  of  the  hardness.  Although  a  high  hard¬ 
ness  minimizes  the  wear  rate  both  below  and  above  the  fracture  threshold, 
the  threshold  force  is  diminished  by  an  increase  in  hardness.  Situations 
can  arise,  therefore,  when  a  hardness  increase  in  the  vicinity  of  the 
threshold  leads  to  an  enhanced  wear  rate,  as  depicted  in  fig.  6. 

The  influence  of  the  fracture  toughness  also  requires  some  qualifica¬ 
tion.  The  toughness  needed  to  resist  lateral  fracture  must  be  developed 

on  the  appropriate  scale.  Specifically,  the  nicrostructural  scale  must 

•  ^ 

be  appreciably  smaller  than  the  length  of  the  lateral  fracture  in  order  to 
experience  the  macroscopic  fracture  toughness.  Some  important  effects  of 
the  dominant  microstructural  size  (e.g.  the  grain  size)  are  thus  to  be 
anticipated,  especially  at  force  levels  in  the  vicinity  of -the  threshold. 
These  effects  will  require  more  detailed  consideration  on  an  individual 
basis. 
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APPENDIX 

LATERAL  CP  >CK  EXTENSION 

Consider  a  two-dimensional  lateral  crack  configuration  created  by  a 
linear  groove  of  length  l  .  The  strain  energy  release  rate  is 

*  (Pp2l)  dX/dc  (A-l ) 

where  Pr  is  the  residual  normal  force  exerted  by  the  plastic  zone  on  the 

surrounding  elastic  material,  X  is  the  compliance  of  the  elastic  layer 

« 

above  the  crack,  and  c  is  lateral  crack  length.  Linear  elasticity  requires 
that^ 

Pr  ■  Pr°[l-ur/uJ]  (A-2) 

where  P°  is  the  initial  residual  force  (at  zero  crack  length),  ur  is 
the  opening  displacement  at  the  crack  center,  and  u°  is  the  opening  dis¬ 
placement  when  the  driving  force  relaxes  to  zero.  For  well  developed  cracks 
(c»h  where  h  is  the  depth  of  the  lateral  fracture)  the  compliance  may  be 
well  approximated  by+ 

X  -  4c3/EAh3  •  (A-3) 

Combining  eqns.  (A-l),  (A-2)  and  (A-3),  the  stress  intensity  factor  becomes; 


K  = 


cP°r 

£h3(l+6c3) 


(A-4) 


Eqn.  (A-3)  is  the  standard  solution  for  a  thin  elastic  cantilever.  The 
cantilever  in  the  present  model  is  disproportionately  wide  and  may  not 
extend  to  the  edges  of  the  material.  The  major  influences  of  these 
deviations  from  ideal  geometry  are  expected  to  be  reflected  in  the  numerical 
factor  in  eqn.  (A-3)  and  in  the  introduction  of  a  Poisson's  ratio  dependence. 


-16- 


where  8  =  (Pj/uj) (4,/E£h3) .  Equating  K  to  the  critical  value  K 

I 

lateral  crack  length  becomes; 

1  /2 

1  +  8c3  =  [G/(l-v2)]  P°c/tK  rh3/Z 

which,  to  the  same  level  of  approximation  used  in  eqn  (3a),  js 
c  *  [G/{1-v2)]1/4  (Eh3/2u°/4Kc)1/2 

Noting  that  (from  eqns  5  and  6) 
h-^b  *  (E/H)2/5(p/H)1/2 
and  that9’^ 

u°  -  AV/£b  *  a2/ Jib 

«  (H/E)2/S  (P/H)1/2  , 
the  lateral  crack  length  becomes 

e  .  ml!l 

c  Hl/s  K  1/2 


as  in  eqn.  (3b). 
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TABLE  I 

LATERAL  CRACK  NUCLEATION  EMMATESt 


Material 

Touohness  (MPa/m) 

Hardness  (GPa) 

I 

1 

Threshold  (N) 

1 

t 

Al2°3 

O  o 

t.  .  u 

20 

0.4  j 

Glass 

0.7 

G 

0.2  1 

MgO 

1.1 

8 

] 

0.9 

Si3N4 

3  (polycrystal) 

16 

31 

Si  3^4 
(estimate) 

• 

"“2  (single  crystal) 

16 

X).G  j 

•f* 

The  threshold  in  most  brittle  materials  occurs  at  load  levels  wherein  the 
scale  of  the  cracks  is  on  the  order  of  (or  smaller  than)  the  grain  size6.  Sinole 
crystal  values  of  toughness  are  thus  pertinent  unless  the  material  is  particularly 
fine  grained.  The  thresholds  are  thus  based  on  single  crystal  toughness  levels 
unless  otherwise  stated. 
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Fig,  1 


Fig.  2 


Fig.  3 

Fig.  4 


Fig,  5 


Fig.  6 


FIGURE  CAPTIONS 

A  schematic  indicating  the  mechanisms  of  material  removal  by 
(a)  lateral  fracture  and  (b)  plastic  cutting. 

The  influence  of  the  normal  force  on  the  extent  of  lateral 

•5 

fracture  , 

12 

The  effect  of  the  applied  force  on  the  material  removal 

1 3 

A  correlation  of  grinding  forces  at  constant  removal  rates 
with  the  toughness,  hardness  parameter  predicted  by  the  present 
analysis.  The  toughness  values  are  in  MPa vF  and  the  hardness 
values  are  the  the  Knoop  hardness, 

1 5 

A  correlation  of  material  removal  rates  at  constant  force 

with  the  toughness,  hardness  parameters  predicted  by  the  present 

analysis. 

A  schematic  indicating  the  wear  rate  as  a  function  of  normal 
force  in  two  materials  with  a  different  hardness  H  but  comparable 
levels  of  toughness  KQ  and  of  E/H. 
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SUMMARY 


The  effects  of  distributions  of  projectile  size  and  velocity  on  the 
erosion  of  brittle  materials  is  assessed.  Strong  dependencies  of  erosion 
rate  on  the  width  of  the  particle-size  distribution  and  on  the  velocity 
characteristics  of  the  erosion  equipment  are  predicted.  The  predictions 
are  confirmed  by  experiment.  The  results  indicate  that  careful  design  of 
controlled  erosion  experiments  is  essential,  in  order  to  avoid  misleading 
predictions  of  in-service  erosion  rates.  The  implications  of  the  analysis 
for  the  interpretation  of  threshold  effects  and  for  experiments  designed 
to  verify  erosion  theories  are  also  discussed. 
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1.  INTRODUCTION 

The  extent  of  solid  particle  impact  damage  in  brittle  materials  has 
been  demonstrated,  both  theoretically  [1-4]  and  experimentally  [1-8],  to 
depend  strongly  on  the  size  and  velocity  of  the  impacting  particles,  as 
well  as  the  relevent  material  properties  (hardness,  toughness,  elastic 
modulus,  density,  acoustic  impedance).  In  any  experimental,  or  in-service 
exposure  involving  multiple  impacts,  the  sizes  and  velocities  of  the  pro¬ 
jectiles  are  distributed  over  a  finite  range.  The  velocity  and  size  dis¬ 
tributions  must  be  incorporated  into  damage  analyses.  Strength 'degradation 
properties  are  simply  determined  by  the  maximum  projectile  size  and  velocity, 
which  dictate  the  most  severe  damage.  Material  removal  processes,  on  the 
other  hand,  are  cumulative  and  are  therefore  determined  by  the  entire 
particle  size  and  velocity  distributions.  Previous  erosion  analyses  have 
been  based  on  velocity  and  size  parameters  defined  simply  in  terms  of  the 
mean  projectile  diameter  and  velocity.  In  this  paper  the  influence  on 
erosion  of  an  additional  parameter,  the  particle  size  variance,  is  assessed. 

The  implications  of  the  analysis  for  the  controlled  experimental  measure¬ 
ment  of  in-service  erosion  rates,  for  experiments  designed  to  test  erosion 
theories,  and  for  the  interpretation  of  threshold  effects  are  discussed. 

The  analysis  is  confined  to  the  elastic/plastic  regime  relevent  to  angular, 
hard  particles  [1-4]. 

Two  fundamentally  different  erosion  test  configurations  are  in  current 
use.  In  one,  the  projectiles  are  accelerated  in  a  gas  stream  and  the  par¬ 
ticle  velocity,  for  a  given  gas  stream  velocity,  is  dependent  upon  the  par¬ 
ticle  size  [9].  In  the  other,  which  simulates  a  target  moving  through  stationar 
particles,  the  relative  velocity  is  independent  of  particle  size  [10].  The 
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influence  of  a  particle  size  distribut  on  would  be  expected  to  differ 
in  these  two  cases.  The  difference  is  assessed  and  the  predictions  con¬ 
firmed  by  experiment. 

2.  ANALYSIS 

The  general  features  of  single  particle  impact  damage  in  the  elastic/ 
plastic  regime  are  relatively  insensitive  to  projectile  shape  and  mechanical 
properties.  The  fracture  pattern  shows  strong  similarity  to  the  damage 
produced  by  quasi-static  indentation  [1,5].  Most  notably,  two  primary 
types  of  fracture  have  been  identified.  Half-penny  shaped  radial/median 
cracks  develop  normal  to  the  target  surface.  These  cracks  control  the 
strength  properties.  Lateral  cracks  form  parallel  to  the  target  surface 
and  are  centered  near  the  base  of  the  plastically  deformed  zone.  Pro¬ 
pagation  of  the  lateral  cracks  to  the  surface  produces  material  removal. 

Existing  erosion  models  are  based  on  analysis  of  the  volume  of 
material  removed  by  the  lateral  cracks  in  single  particle  impact.  Inter¬ 
action  effects  are  assumed  to  be  negligible,  so  that  the  cumulative  effect 
of  multiple  impacts  is  obtained  by  summing  the  volumes  removed  by  individ¬ 
ual  (independent)  impacts.  Two  quantitative  models  have  been  developed 
for  predicting  the  erosion  rate,  E,  (volume  loss  per  impact).  One  is 
based  upon  analysis  of  quosi-static  indentation  [1,2]  and  employs  an 
upper  bound  quasi-static  impulse  load  (determined  by  equating  the  projec¬ 
tile  kinetic  energy  to  the  work  required  to  produce  the  elastic/plastic 
indentation).  The  resulting  erosion  rate  is 

E  =  v2,4  R_3* 7  o  1,2  r1-3  h0*1  / 1  \ 


(1) 
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where  v  ,  Rp  ,  and  pp  are  the  projectile  velocity,  diameter  and  densi  .y, 
and  Kc  and  H  are  the  target  toughness  and  hardness.  The  other  model 
is  based  on  a  dynamic  analysis  of  the  elastic/plastic  stress  field,  and 
gives  the  result  [3,4]; 

E  =  v3*2  R  3,7  p  0,25  K"1’3  H-0*25  (2) 

P  c 

Despite  the  different  physical  assumptions  the  two  theories  predict  similar 

♦  • 

erosion  rates.  Here  we  confine  our  interest  to  particle  size  and  velocity 
effects  in  a  given  projectile  target  system  and  express  equations  1  and  2 
as; 

E  =  k  vn  Rm  (3) 

where  target  material  parameters  are  included  in  the  constant  k  .  The 
influences  of  projectile  size  and  velocity  distribution  on  the  erosion 
rate  follow  directly  from  the  explicit  size  and  velocity  dependence  of 
eqn.  3. 

Particle  sizes  generally  conform  to  a  logarithmic-normal  distribu¬ 
tion  [11]: 

4>(R)dR  =  [1/c(2tt)1/2]  exp  [-(£nR  -  £nRu)2/2o2]d£nR  (4) 

where  <j>(R)  JR  is  the  fraction  of  particles  with  diameter  between  R  and 

2 

R+dR  ,  Rp  is  the  median  particle  diameter,  and  a  is  the  variance  of 
£nR  (Fig.  1).  In  the  following  analysis  it  is  assumed  that  all  particles 
of  a  given  diameter  impact  the  target  at  the  same  velocity,  but  the 
impact  velocity  is  allowed  to  exhibit  some  dependence  upon  the  particle  diameter 
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(Appendix).  Hence,  by  adopting  a  relation  between  velocity  and  size 
(Appendix,  eqn.  Al),  the  erosion  rate  (eqn.  3)  can  be  expressed  exclusively 
in  terms  of  the  particle  diameter; 

E(R)  •  k  a"  Rnx+m  (5) 

where  x  and  a  are  characteristic  of  the  erosion  test  configuration. 

The  erosion  resulting  from  a  distribution  of  projectile  sizes  can  be 
directly  assessed  from  eqns.  4  and  5. 

The  erosion  rate  is  generally  measured  in  terms  of  the  total  volume 
of  target  removed,  normalized  by  the  calculated  number  of  impacts.  The 
calculated  number  of  impacts,  N  ,  (in  general  not  equal  to  the  actual 
number  NQ)  is  expressed  in  terms  of  the  measured  total  mass  of  projectiles 
divided  by  the  mass  of  the  projectile  with  the  median  diameter: 


at 

*  No  / 
r w 


R30(R)dR/R 


With  the  logarithmic-normal  size  distribution  (eqn.  4),  eqn.  (6)  becomes 

N  =  NQ  exp(9a2/2)  .  (7) 

» 

The  total  volume  of  target  removed  by  NQ  impacts  is 

oo 

V  »  f  E(R)  N  <J>  (R)  dR  ,  (8) 

0“' 

which,  in  conjunction  with  eqns.  (4)  and  (5)  becomes 

V  =  N  k  an  Rnx+17)  exp[02(nx+m)2/2] 


(9) 
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The  measured  erosion  rate  is  therefore, 

E  =  E(Ry)  exp  (a2/2)[(nx+m)2-9]  (10) 

where  E(Ry)  =  k  is  the  erosion  rate  corresponding  to  projectiles 
with  the  median  diameter  (i.e.  a  =0).  A  potentially  strong  dependence  of 
measured  erosion  rate  on  a  is  apparent. 

3.  EXPERIMENTAL 

In  order  to  test  the  predicted  dependence  of  erosion  rate  on  particle- 
size  distribution,  several  distributions  of  SiC  particles  were  prepared. 

The  particle  sizes  within  each  distribution  conformed  approximately  to  a 
log-normal  function,  with  a  median  diameter  of  38  ynr,  but  the  variances 
ranged  from  0.2  to  0.6  (Fig.  2).  The  distributions  were  prepared  by  mixing 
commercial  SiC  powders  (grit  sizes  120,  150,  180,  240,  280,  400,  600)  in 
the  appropriate  ratios,  as  calculated  from  the  measured*  size  distributions 
of  each  constituent. 

Erosion  experiments  were  performed  using  both  of  the  test  techniques 
described  in  the  Appendix.  Single  crystal  silicon7  was  chosen  as  a  target 
material  because  complications  associated  with  microstructural  variability 
and  subcritical  crack  growth  can  be  exclud-d,  and  because  impact  damage  in 
the  SiC/Si  projectile/ target  system  conforms  well  to  the  elastic/plastic 
damage  scheme  described  in  section  2[12].  Target  slabs  measuring  25  x  25 
x  5  mm,  with  (111)  planes  parallel  to  the  face,  were  eroded  at  normal 
incidence  and  at  a  projectile  velocity  of  100  m 

★ 

by  sedimentation;  sedigraph-L,  Micromeritics  Corp. 

^Monsanto  Commercial  Products,  St.  Peters,  Missouri. 
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The  erosion  rates  were  determined  from  measurements  of  the  target 

t 

mass  loss  and  the  mass  of  projectiles.  Several  successive  measurements 
were  performed  for  each  experiment  to  ensure  that  the  erosion  rate  was 
steady.  Typical  weight  loss  results  for  both  erosion  devices  are  shown 
in  Fig.  3, 

The  measured  erosion  rates  are  listed  in  Table  I  and  plotted  in  nor¬ 
malized  form,  as  a  function  of  a  ,  in  Fig.  4.  A  strong  dependence  of 
erosion  rate  on  a  is  evident.  The  solid  lines  in  Fig.  4  are  predictions 
based  on  equation  (10)  with  the  following  parameters;  n  =  3,  m  =  3.7  (ex¬ 
perimental  determination  for  Al203/Si  system  [8]),  x  -  0  for  the  slinger 
device,  and  x=  -0.16  for  the  gas  stream  apparatus  (Fig,  Al).  The  normalizing 
factor,  E(Ry)  ,  (the  erosion  rate  corresponding  toa=  0)  was  calculated 
from  equation  (10)  using  the  same  values  of  n,  m  and  x,  and  the  erosion 
rate  measured  for  batch  A  particles.  This  procedure  gave,  E(R^)  =  246  ym3/  impact 
for  the  gas  stream  apparatus,  and  E(Ry)  =  305  um3/impact  for  the  slinger 
device.  With  due  allowances  for  the  uncertainty  in  the  particle  size 
variances,  the  predictions  and  experiments  show  good  agreement. 

4.  DISCUSSION 

4.1  General  Discussion 

The  dependence  of  erosion  rate  on  particle  size  distribution 
is  sensitive  to  the  values  of  the  particle  velocity  and  size  exponents  n,  m 
and  x  (eqn.  10).  Consequently  the  absolute  erosion  rates  measured  by  the  two 
erosion  devices  (under  otherwise  identical  conditions)  are  significantly 
different  (Table  I).  It  is  of  interest  to  examine  this  difference  in  more 
detail.  The  errors  listed  in  Table  I  represent  the  scatter  in  measurements, 
and  are  therefore  appropriate  for  comparison  of  results  from  a  particular 
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appa' atus  under  constant  operating  conditions.  Comparison  of  results 
from  the  tvo  erosion  devices  requires  the  additional  systematic  error 
due  to  uncertainty  in  velocity  calibration  to  be  accounted.  Taking  E°=v^, 
along  with  5%  velocity  error  for  the  gas  stream  apparatus  and  2%  velocity 
error  for  the  slinger  device,  the  required  erosion  rate  errors  are  «*20% 
and  *10%  respectively.  With  these  errors  all  of  the  results  in  Table  I 
still  differ  significantly.  However,  the  calculated  erosion  rates  at  a=0 
(i.e.  E(Ry)  3  246  t 50  pm'Vimpact  for  the  gas  stream,  and  E(Ry)  =  305  ± 30  pm^/ 
impact  for  the  slinger  device)  are  consistent,  indicating  that  the  apparent 
inconsistency  in  Table  I  is  due  to  the  different  dependencies  of  E  on  a, 
arising  from  the  different  values  of  x. 

An  underlying  assumption  of  the  present  analysis  is  that  interaction 
affects  are  negligible,  so  that  the  cumulative  effect  of  multiple  impacts 
is  given  by  summing  the  volumes  removed  by  individual  (independent)  impacts. 
The  agreement  shown  between  prediction  and  experiment  would  appear  to  lend 
some  support  to  the  assumption.  The  subject  of  interaction  effects  has  been 
addressed  in  two  prior  studies.  In  the  first.  Gulden  [13]  showed  that 
interaction  effects  between  particles  of  uniform  size  are  negligible,  by 
demonstrating  the  absence  of  an  incubation  period  in  the  erosion  of  polished 
surfaces.  However,  an  enhanced  erosion  rate  was  observed  for  the  initial 
impacts  by  50  pm  particles  on  a  surface  which  had  been  previously  eroded 
by  270  pm  particles.  This  result  is  consistent  with  the  observation  by 
Routbort  et  al  [8]  of  an  enhanced  erosion  rate  from  a  mixture  consisting  of 
equal  weight  fractions  of  40  pm  and  270  pm  particles.  These  results  imply 
that  the  erosion  rate  for  small  particles  can  be  increased  by  the  prior 
damage  produced  by  much  larger  particles.  In  the  present  experiments  most 
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of  the  material  removal  is  produced  by  the  larger  particles  in  the  dis¬ 
tributions,  so  that  any  enhancement  of  the  contributions  of  the  smaller 
particles  is  not  expected  to  influence  the  total  erosion  rate  significantly. 
An  upper-bound  estimate  indicates  that  the  enhancement  is  less  than  2%. 


Imp! ications 


4.2.1  Simulation  of  In-Service  Erosion 


The  strong  dependence  of  E  on  a  has  important  implications 
for  the  use  of  controlled  experiments  to  simulate  in-service  erosion  rates. 
Firstly,  a  controlled  experiment  using  projectiles  matched  in  median  size 
to  the  in-service  projectiles,  but  with  a  narrower  distribution  of  sizes, 
would  result  in  a  serious  underestimate  of  erosion  rate.  The  underestimate 


is  largest  for  situations  most  closely  simulated  by  the  slinger  device 
(e.g.  target  moving  through  a  stationery  dust  cloud)  rather  than  the  gas 
stream  apparatus  (e.g.  stationary  target  subjected  to  dust  storm).  Secondly, 
use  of  the  gas  stream  apparatus  for  predict  erosion  rates  for  in-service 
conditions  approximating  the  slinger  device  will  also  result  in  an  under¬ 
estimate  (even  if  the  actual  in-service  projectiles  are  used).  Therefore 
careful  design  of  controlled  erosion  experiments  is  essential. 

4.2.2  Evaluation  of  Erosion  Parameters 

The  result  in  eqn.  (10)  has  potential  implications  for  experi¬ 
ments  designed  to  evaluate  the  parameters  m  and  n  in  eqn.  (3).  These 
parameters  are  usually  determined  from  the  gradients  of  logarithmic  plost  of  E  v 
at  constant  v,  and  E  vs.  v  at  constant  Ry.  However  it  is  imme¬ 
diately  evident  from  eqn.  (10),  that  if  x  or  a  vary  between  experiments, 
the  erosion  rate  will  not  show  the  functional  dependence  of  eqn.  (3).  Con¬ 
sider  first  experiments  involving  constant  particle  size  and  varying  veloc¬ 
ities.  Then  a  is  constant,  while  x  is  zero  for  the  slinger  device  and 
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non-zero  (but  constant)  or  the  gas  stream  apparatus  (Fig.  A1).  The  correct 
dependence  of  erosion  rate  on  velocity  is  thus  obtained  with  both  erosion 
devices.  This  result  contrasts  with  that  expected  for  experiments  in  which 
the  velocity  is  held  constant  but  the  particle  size  varies.  Then,  provided 
that  the  variance  remains  constant  for  each  particle  size,  the  correct 
dependence  of  erosion  rate  on  particle  size  is  obtained  for  the  slinger 
device,  but  not  for  the  gas  stream  apparatus.  (The  constant  variance  condi¬ 
tion  is  reasonably  well  satisfied,  for  example,  by  commercially  sized  Si C 
powders  (era*  0.2)  „)  For  the  gas  stream  apparatus,  where  x  varies  between 
-0.13  at  Ry  =  50  ym,  and  -0.32  at  Ry  =  1100  ym  (Fig.  Al),  the  size  exponent 
of  the  erosion  rate  deviates  from  the  required  value,  m  .  However  for  the 
typical  values  of  the  erosion  parameters  used  in  this  paper  (n  =  3,  m  =  3.7) 
an  insignificant  error  (*»1%)  in  the  experimental  determination  of  m  would 
result.  Therefore  it  appears  that,  provided  care  is  taken  to  use  constant 
width  particle  distributions,  the  conventional  experimental  determinations 
of  m  and  n  yield  accurate  results. 

4.2.3  Interpretation  of  Threshold  Effects 

The  existence  of  a  finite  distribution  of  projectile  size 

i 

bears  on  the  interpretation  of  threshold  effects  in  erosion.  Several 
studies  have  established  that  eqn.  (3)  is  well  satisfied  for  relatively 
large  particle  sizes  and  velocities.  However,  as  a  threshold  is  approached, 
the  erosion  rate  decreases  continuously  below  a  linear  extrapolation  from 
large  size-velocity  data  (Fig.  5)  [2,8].  This  behavior  can  be  described 
empirically  by  modifying  eqn.  (3)  to  incorporate  a  threshold  velocity  vQ 
and  size  RQ  [8]; 

E  ■  k(v-v0)"  (R-R0)"  <") 
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where  vQ  and  RQ  are  determined  by  fitting  eqn.  (11)  to  experimental 
data.  Equation  (11)  implies  that  the  extent  of  cracking  in  single  particle 
impacts  increases  continuously  from  zero,  as  the  excess  driving  force  above 
the  threshold  increases,  as  shown  schematically  by  curve  1  in  Fig.  6.  Single 
particle  impact  and  indentation  studies  have  verified  the  existence  of  a 
threshold  below  which  no  cracking  occurs  [6,12,13].  However,  at  the  thresh¬ 
old,  lateral  cracks  usually  initiate  discontinuous,^  as  shown  by  curve  2  in 
Fig.  6.  Erosion  with  uniformly  sized  particles  would  thus  be  expected  to 
exhibit  a  corresponding  discontinuous  increase  at  the  threshold.  However, 
a  distribution  of  particle  sizes  would  yield  a  continuous  increase,  because 
(near  the  threshold)  only  that  fraction  of  particles  with  radii  in  excess 
of  the  threshold  radius  would  remove  material.  A  discontinuous  cracking 
threshold  can  be  incorporated  into  the  analysis  of  section  2  by  replacing 
the  lower  integration  limit  in  eqn.  (8)  by  the  threshold  size  R0.  Comparison 
of  the  corresponding  prediction  with  data  from  Ref.  [8]  in  Fig.  5  shows 
good  agreement  (similar  to  the  agreement  obtained  from  eqn.  (11)).  Distinc¬ 
tion  of  the  two  threshold  behaviors  depicted  in  Fig.  6  is  not  possible 
from  the  erosion  data,  but  the  discontinuous  threshold  is  considered  to  be 
the  more  plausible  alternative. 

5.  CONCLUSIONS 

The  rate  of  material  removal  in  the  elastic/plastic  damage  regime  is 
sensitive  to  both  the  distribution  of  projectile  sizes  and  the  velocity 
characteristics  of  the  erosion  apparatus.  Consequently,  experimental 
erosion  measurements  can  appreciably  underestimate  in-service  erosion  rates 
if  the  experimental  projectile  size  distribution  is  narrower  than  that  of 
the  in-service  projectiles,  or  if  the  velocity  characteristics  of  the 
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experimental  and  in-service  erosion  configurati jns  differ. 

Results  of  experiments  designed  to  measure  the  velocity  and  particle 
size  exponents  in  the  single  particle  erosion  expression  can  also  be  in¬ 
fluenced  by  the  size  distribution.  However,  under  commonly  used  test 
conditions  the  influence  is  insignificant,  so  conventional  measurements, 
which  ignore  the  effect,  yield  accurate  results. 

A  straightforward  interpretation  of  erosion  threshold  behavior,  in 
terms  of  observed  single-particle  damage  initiation  characteristics  in 
conjunction  with  a  projectile  size  distribution,  is  suggested. 
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APPENDIX 


Velocity/particle-size  characteristics  of  erosion  experiments 

Two  different  experimental  erosion  facilities  are  used  to  investigate 
the  effects  of  particle  si2e  distribution  on  erosion.  In  one  the  projectiles 
are  introduced  into  a  uniformly  moving  gas  stream  in  a  0.95  cm  diameter 
tube,  and  accelerated  over  a  distance  of  3  m  before  impacting  the  target. + 

The  calculated  velocities  for  a  range  of  particle  diameters  and  gas  stream 
velocities  are  shown  in  Fig.  Al.  These  calculations  have  been  verified 
by  direct  experimental  measurement  using  the  double  rotating  disc  method 
[9].  For  a  given  gas  velocity,  the  particle  velocity  can  be  approximated 
by 

v  *  a  Rx  (Al ) 

where  a  and  x  can  be  taken  as  constants  for  a  given  experiment,  pro¬ 
vided  that  the  range  of  particle  diameters  is  restricted.  Over  the  full 
range  of  particle  sizes  shown  in  Fig.  Al  (a  factor  of  100)  x  varies 
between  -0.13  and  -0.32;  in  a  typical  experiment  the  range  of  particle 
diameters  is  about  a  factor  of  2. 

★ 

The  other  erosion  method  utilizes  a  "slinger"  device  [10],  in  which 
the  projectiles  are  introduced  into  the  center  of  a  rotating  tube  with 
open  slits  at  both  ends.  After  moving  slowly  along  the  tube  the  projectiles 
are  thrown  from  the  slits  in  a  tangential  direction,  with  the  speed  of  the 
tube  ends,  and  impact  stationary  targets.  The  experiment  is  performed  in 
vacuum.  Rotating  double  disc  measurements  have  verified  that  the  particle 
velocity  is  the  same  as  the  velocity  of  the  tube  ends,  and  independent  of 

+ Erosion  test  facility  at  Solar  Turbines .[14], 


*Erosion  facility  at  Argonne  National  Laboratory  [15]. 
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particle  size.  Therefore,  for  the  purposes  of  analysis,  the  partic’i 
velocity  can  be  expressed  by  eqn.  (A1)  with  x  =  0. 


i 


*  ***_  n 
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Table  1 


Erosion  of  single  crystal  silicon  by  the  silicon  carbide  particles 
of  Fig.  3.  Normal  incidence,  velocity  100  M  S~\ 


Erosion 

Erosion  Rate 

3 

urn  /impact 

Apparatus 

Batch  A 

Batch  B 

Batch  C 

Batch  D 

Gas  Stream 

254  ±10 

262  ±10 

325  ±10 

Slinger 

340  ±10 

492  ±10 

527  ±10 

646  ±10 

-  1.7  - 

FIGURE  CAPTIONS 

1.  Typical  probability  distribution  for  particle  distributions  (logarithmic- 
normal). 

2.  Size  distribution  for  three  batches  of  SiC  particles  used  in  erosion 
experiments. 

3.  Typical  measurements  of  target  weight  loss  as  a  function  of  mass  of 
projectiles.  Batch  A  SiC  particles,  v  =  100  ms"\  normal  incidence, 
silicon  target. 

4.  Erosion  rate  as  a  function  of  particle  distribution  width  for  two  erosion 
devices.  Data  obtained  from  the  SiC  particles  of  Fig.  3  impacting  a 
silicon  target  at  normal  incidence  and  v  =  100  ms~^ .  Solid  curves  are 
predictions  from  eqn.  (11). 

5.  Erosion  threshold.  Data  from  Ref.  [4];  A1203  projectiles  impacting 
silicon  target,  normal  incidence,  v  =  107  ms’1.  Straight  line  fitted 

to  data.  Solid  line  for  R£40  ym  predicted  from  the  analysis  of  section 
2.2,  with  lower  integration  limit  in  eqn.  (9)  replaced  by  RQ  =  27  ym. 

6.  Schematic  representation  of  single  particle  contact  damage  near  the 
threshold. 

Al.  Velocity/particle-size  relation  for  gas-stream  erosion  apparatus. 
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